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Abstract 



The influence of mesonic fluctuations on quantities in the Nambu-Jona-Lasinio model is 
examined. To that end different approximation schemes are introduced which guarantee 
the consistency with several relations following from chiral symmetry. In particular this 
refers to the Goldstone theorem, which states the existence of massless bosons, in our 
case pions, if the symmetry of the Lagrangian is spontaneously broken. The Gell-Mann- 
Oakes-Renner relation describes the behavior of the pion mass for small current quark 
masses, which explicitly break chiral symmetry. Three different schemes are presented. 
These are the inclusion of the ring sum in a so-called "$-derivable-method", an expansion 
in powers of 1/7V C and an expansion up to one- meson loop in the effective action formalism. 
Since the two latter schemes are used for explicit calculations, it is explicitely proved that 
the Goldstone theorem as well as the Gell-Mann-Oakes-Renner relation hold within those 
schemes. 

The influence of meson- loop effects on the quark condensate {4>ip), the pion mass, the 
pion decay constant f n and properties of p- and cr-meson are investigated. First we focus 
on the determination of a consistent set of parameters. In the 1 /A^-expansion scheme 
it is possible to find a set of parameters which allows to simultaneously describe the 
quantities in the pion sector and those related to the p-meson, whereas this turns out to 
be not possible within the expansion of the effective action. Results for the a-meson are 
also discussed. Here, similarly to the p-meson, mesonic intermediate states are essential. 
Besides, the relation of our model to hadronic models is discussed. 

In the last part of this thesis the behavior of the quark condensate at nonzero temper- 
ature is studied. In the low-temperature region agreement with the model independent 
chiral perturbation theory result in lowest order can be obtained. The perturbative 1/N C - 
expansion scheme does not allow for an examination of the chiral phase transition at higher 
temperatures, whereas this is possible within the meson-loop approximation scheme. A 
first order phase transition is found. 
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Introduction 



Since the substructure of the nucleus was discovered at the beginning of the century one 
tries to understand in which way its constituents, the nucleons, interact. Already in the 
thirties Yukawa suggested U that this interaction could be described by an exchange of 
a massive meson which was later identified with the pion. Presently this idea is still the 
basis for most phenomenological models of the interaction between nucleons which mainly 
rely on the exchange of mesons with various quantum numbers |2], |||] . 

Not only our present notion of the interaction between nucleons, but almost all phe- 
nomenological models which try to describe strongly interacting particles, mesons and 
baryons, in the low-energy region, are essentially influenced by the idea of Yukawa. These 
models incorporate mesons and baryons as the principal degrees of freedom in the low- 
energy region. In many cases reasonable descriptions of hadronic spectra, decays and 
scattering processes have been obtained within these phenomenological models. For in- 
stance, the pion electromagnetic form factor in the time-like region can be reproduced 
rather well within a simple vector dominance model with a dressed p-meson which is 
constructed by coupling a bare p-meson to a two-pion intermediate state j|, ^[ . 

Recently much experimental and theoretical effort has been undertaken to understand 
the properties of hadrons not only in vacuum but also in a strongly interacting medium 
and, in in the same context, the phase structure of strongly interacting matter. Besides for 
nuclear physics, these questions are of importance mainly for astrophysical applications: 
In the evolution of the early universe one presumes hot strongly interacting matter with 
small baryon densities having existed (temporarily) shortly after the "big bang" , and in 
the interior of neutron stars one expects to find densities of a few times nuclear matter 
density at very low temperatures. Experimentally one tries to gain information about the 
state of matter under various conditions with the help of heavy-ion experiments. Since 
experimentally only secondary hadrons are measurable, it is extremely difficult to obtain 
direct information about the state of matter in the early stages of the collision. Dilepton 
spectra (e + e~ or /i + p~) offer a possibility to learn, at least indirectly, something about the 
early stages, because leptons do not interact strongly. In this context, the vector mesons, 
especially the p-meson, play an important role. Here phenomenological hadronic models 
have been successfully extended to investigate medium modifications of vector mesons 
and to calculate dilepton production rates in hot and dense hadronic matter |J. 

However, hadrons do not represent the fundamental degrees of freedom. Nowadays 
quantum chromodynamics (QCD) is generally accepted as the theory of strong interaction, 
which contains quarks and gluons as fundamental degrees of freedom. QCD is constructed 
in the sense of a Yang-Mills gauge theory 0]. It is based on the color S'f/(3)-gauge group 
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with the gluons representing the gauge bosons. From the non-abelian structure of the 
gauge group it follows that the gluons themselves carry color charge and therefore can 
interact with each other. Quarks exist in six different flavors, "up", "down", "strange", 
"charm", "bottom" and "top", whereof the two first are almost massless compared with 
the other ones and typical hadronic mass scales. 

Two prominent features of QCD are "confinement", which describes the fact that 
in nature only colorless objects are observed, and the so-called "asymptotic freedom". 
The former means in particular that no free quarks and gluons can be observed, they 
are "confined" to hadrons. Up to know the origin of this phenomenon in QCD has not 
been understood. This is in contrast to asymptotic freedom, which reflects the fact that 
the coupling constant decreases with increasing energies. As a consequence QCD can 
be treated perturbatively at sufficiently high energies, whereas in the low-energy region, 
which is the relevant one for our hadronic world, perturbative methods are not applicable. 

In principle lattice calculations (see e.g. Ref. ||) provide us with the possibility to 
obtain results directly from QCD also at lower energies. Among the successes of these 
ab-initio lattice calculations are the simulation of the linear rising (confining) potential 
between two heavy quarks as well as the description of the mass spectrum of hadrons, and 
in particular calculations at nonzero temperatures. However, up to now these calculations 
mainly have dealt with thermodynamic observables. In the sector of light quarks, up 
and down, on the contrary, one is less successful, because of fundamental problems with 
treating light fermions on the lattice. In addition, till now no reliable results have been 
obtained at nonzero density. 

According to the Goldstone theorem |J a symmetry of the Lagrangian which is spon- 
taneously broken in the ground state enforces the existence of massless bosons. In the 
present case this symmetry is the chiral SU (2) l ® SU (2) R-symmetry in flavor space which 
is exhibited by QCD in the sector of up- and down-quarks. The small mass of the pions 
as compared with other hadronic masses is thought to reflect their nature as Goldstone 
bosons. The nonzero mass is thereby attributed to the fact that this symmetry is explicitly 
broken because of the nonzero masses of up- and down-quarks. The spontaneous breaking 
of chiral symmetry in the physical vacuum is usually attributed to instantons QlO | , semi- 



classical gluon field configurations. A method which takes advantage of the small mass of 
the pions and the fact that Goldstone bosons interact only weakly is "Chiral perturbation 



theory" |TT| , which corresponds to an expansion in powers of momenta and pion masses. 
However, this method also has its shortcomings: It is e.g. not suited to treat resonances. 
Thus, besides the small numbers of observables which can be studied on the lattice one 
has to rely here on calculations in effective models. 

Though hadronic calculations prove very successful in describing the properties of light 
hadrons, one might certainly ask in which way they can be understood from the underlying 
quark substructure. Since this question cannot be answered from first principles it has to 
be addressed within quark models. The main objection one can raise against most of the 
known quark models is the lack of confinement. However, at least for light hadrons chiral 
symmetry and its spontaneous breaking in the physical vacuum seems to play the decisive 



role in describing their properties with confinement being much less important [10|. Two 



of the most prominent examples of quark models incorporating chiral symmetry were 



already invented in the sixties: The Gell-Mann-Levy model [12|, dealing originally with 
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nucleons, pions and cr-mesons, and the Nambu-Jona-Lasinio (NJL) model [13], which 



originally was a model which exclusively contained nucleons. Nowadays these models are 
reinterpreted with quarks instead of nucleons. 

The NJL model has been applied by many authors to study properties of hadrons. 
Baryons can be constructed either by directly building a bound state of three quarks 



by solving Fadeev equations [14] or by considering chiral solitons [15]. Mesons of vari- 



ous quantum numbers have already been investigated by Nambu and Jona-Lasinio JTB 



and by many authors thereafter (for reviews see |T6|, [17|, ^8[). In most of these works 
one starts by calculating quarks in mean-field approximation (Hartree or Hartree-Fock) 
and then constructs mesons as correlated quark-antiquark states in "Random Phase Ap- 
proximation" (RPA). With a suitable choice of the model parameters, chiral symmetry is 
spontaneously broken in the vacuum and pions emerge as massless Goldstone bosons. The 
spontaneously broken symmetry is also reflected by the nonzero value of the quark con- 
densate, being closely related to the nonzero mass the constituent quarks acquire within 
this approximation. 

The breaking of chiral symmetry in the vacuum via dynamical effects and the consis- 
tent description of the pion is certainly one of the successes of the model. In contrast, 
the description of other mesons is more problematic, mainly because of the missing con- 
finement mechanism in the NJL model. This means that a meson can decay into free 
constituent quarks if its energy exceeds twice the constituent quark mass m. A typical 
parameter set, which is chosen in such a way that the values for the pion decay constant 
/„■, the pion mass m^, and the quark condensate are in agreement with the empirical 
values, leads to a quark mass of about m ~ 300 MeV. Thus, for instance the p-meson 
with a mass of 770 MeV would be unstable against decay into quarks. This is of course an 
unphysical feature. Besides, the physically most important decay channel of the p-meson 
into two pions is not included in those calculations. The same argument could be put for- 
ward for other mesons like e.g. the cr-meson which emerges in the standard approximation 
scheme as (almost) sharp particle just above the threshold for the decay into a quark- 
antiquark pair, although in nature it -if it exists at all- has a very large width which can 
be mainly attributed to the decay into two pions. The appearance of unphysical effects 
on the one hand and the missing of physical ones on the other hand obviously results in 
a poor description of the properties of those mesons and related quantities, e.g. the pion 
electromagnetic form factor which is mainly determined by the p-meson. 

The phase structure of strongly interacting matter in the NJL model within a mean- 
field calculation has also been investigated by many authors (see e.g. |L7|, pl| [19|, pOp. In 
this way many of the features induced by chiral symmetry, among others the expected 
restoration of this symmetry at nonzero temperatures and densities, can be illustrated 
nicely. 

At this point we should mention that an examination of the chiral phase transition 
is not the only thermodynamical application which has been considered up to know in 



the NJL model within the standard approximation scheme. For instance, in Ref. |21j a 
hypothesis raised by Farhi and Jaffe [22| that strange quark matter could be the absolute 
ground state of matter is critically examined with the use of the SU(3) version of the NJL 
model, which includes in addition to up- and down quarks also strange quarks. In the same 
context one could mention that a few years ago the existence of a "new" state of strongly 
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interacting matter was proposed [23, M: In cold dense matter a "color superconducting" 
phase is likely to exist. In that phase Cooper pairs of quarks and antiquarks are built 
similarly to electronic Cooper pairs in metallic superconductors (see e.g. Ref. p3| and 
references therein). Most of the corresponding calculations have so far been performed 
either in instanton models or in models of the same type as the NJL model. 

However, the same objections which can be raised against the description of mesons 
concern also the investigations of properties below the phase transition. These calculations 
suffer from the fact that the thermodynamics is entirely driven by unphysical unconfmed 
quarks even at low temperatures and densities. On the contrary the physical degrees of 
freedom, in particular the pions, are missing. 

Though there are some applications where the presence of unphysical free quarks 
and the missing of relevant degrees of freedom seem to be less important, e.g. the study 
of bulk quark matter, for most calculations the standard approximation scheme proves 
insufficient. Therefore several authors attempted to extend the standard scheme. For 
instance, in Ref. |26] the coupling of a quark- ant iquark p-meson to a two-pion state 
via a quark triangle was considered. Also higher-order corrections to the quark self- 
energy [^7j and to the quark condensate [^8] were investigated. However, these attempts 
are insufficient in another sense: Since one of the most important features of the NJL 
model is its chiral symmetry, any approximation scheme which is applied should conserve 
properties connected to that symmetry. Especially the existence of Goldstone bosons 
should be ensured if the symmetry is spontaneously broken. 

To our knowledge Dmitrasinovic et al. |29| for the first time discussed a symmetry 
conserving approximation scheme in the NJL model which goes beyond the standard 
Hartree (-Fock) + RPA scheme. The authors included a local correction term to the 
quark self-energy and, by an appropriate choice of diagrams, succeeded to construct meson 
propagators consistently. Among other things this guarantees that the Goldstone theorem 
holds. In addition to the explicit proof of the existence of massless Goldstone bosons, the 
authors of Ref. [29] showed the validity of the Goldberger-Treiman relation. A more 



systematic derivation was presented by Nikolov et al. |3"D|] , who used a one-meson-loop 
approximation to the effective action in a bosonized NJL model. 

An appealing feature of that scheme is that it incorporates mesonic fluctuations. An 
examination of the effect of mesonic fluctuations on various quantities, for instance on the 



pion electromagnetic form factor [Bl[ or on 7r-7r-scattering in the vector [32] and scalar 



channel ||33|| , has been performed, based on that approximation scheme. However, because 
of the difficulties which arise in connection with the numerical evaluation of the occurring 
multi-loop integrals, the authors of these references approximated the exact expressions 
by low-momentum expansions. 

Are there other possibilities to construct symmetry conserving approximation 
schemes? We could show that the most direct extension of the usual RPA in terms of 
techniques used in many-body-theory, the so-called "Second RPA" (SRPA) pi| , does not 
preserve the necessary symmetries |35| . On the other hand, in the literature other symme- 
try conserving approximation schemes are known. In addition to the one-loop expansion of 



the effective action [|3(| mentioned above, these are the "<3>-derivable" -method |37|, and 



an expansion in powers of 1/N C , the inverse number of colors. The mean-field (Hartree) 
approximation in combination with the RPA corresponds to the leading order in such an 
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expansion. In Ref. [29] it was shown that in the NJL model the Goldstone theorem also re- 
mains valid in an expansion up to next-to-leading order in 1/N C . A fundamental difference 
of this approximation scheme to that mentioned above is the perturbative character of 
the correction terms taken into account. More detailed investigations of quark and meson 
properties within the l/7V c -expansion scheme were performed in Refs. p9| , |4*0f . Recently 
such an expansion has been discussed also in the framework of a non-local generalization 
of the NJL model Q. 

In this paper we will calculate quark and meson properties as well in the non- 
perturbative scheme introduced in Refs. as in the l/iV c -expansion scheme, in- 

cluding the full momentum dependence of all expressions. A comparison of the results 
obtained in both schemes, focusing on the pion and the p-meson in vacuum, can also be 



found in Ref. [f42|. The investigation of pion properties was mainly motivated by recent 
papers by Kleinert and van den Bossche |43| who claim that in the NJL model chiral 
symmetry is restored in vacuum due to strong mesonic fluctuations. One argument which 
can be raised against this conjecture |39[] is the non-renormalizability of the NJL model 
which causes new divergencies to emerge if further loops, in this case meson loops, are 
included. And, new divergencies require new cutoff parameters. Following Refs. p9 



and [0 we introduce a cutoff Am for the meson-loops which is independent of the reg- 
ularization of the quark loops. For small values of that cutoff the results in the pion 
sector change only quantitatively as compared with those obtained in the Hartree + RPA 
scheme, whereas for large values of the cutoff within both extended schemes instabilities 
in the pion propagator are encountered which might be a hint for an unstable ground 
state [39, 42]. However, in the non-perturbative scheme a closer analysis revealed that 



these instabilities are not the consequence of restoration of chiral symmetry. 

Encountering those instabilities one might fear that thereby a description of meson 
properties will be made impossible. But since this is an effect which strongly depends on 
the choice of the model parameters, especially the meson- loop cutoff Am, it is not ruled 
out that mesons can be reasonably described with a parameter set far away from the 
region where instabilities appear. In fact, in the l/A c -expansion scheme we determined 
such a parameter set [|0| by fitting the values of m^, f^, the quark condensate (0'0)> an d 
the pion electromagnetic form factor. The inclusion of meson loops, in particular pion 
loops, is absolutely crucial to obtain a realistic description of the latter, which is therefore 
well suited to fix our parameters. 

We mentioned above that one disadvantage of hitherto performed calculations in the 
NJL model is the existence of quark-antiquark decay channels for all mesons if their energy 
is larger than twice the constituent quark mass. Of course, our calculations cannot cure 
this problem, but it can be by-passed: The analysis in Ref. f40| shows that a realistic 
description of the above listed observables is possible with a relatively large constituent 
quark mass of about 450 MeV, such that the threshold for the decay into a gg-pair lies 
above the peak of the p-meson spectral function. This is an important result since the 
constituent quark mass is not an independent input parameter. The same analysis was 



performed in Ref. H2] for the selfconsistent scheme. It turned out that no fit can be 
achieved with a constituent quark mass large enough to shift the threshold above the 
p-meson peak. 

At low temperatures and low densities pionic degrees of freedom are expected to be the 
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dominant ones. Thus not only the description of mesons but also the thermodynamics of 
the model should be considerably improved by including meson loop effects. This becomes, 
for instance, obvious from the temperature dependence of the quark condensate. For its 
low-temperature behavior chiral perturbation theory provides us with model independent 
results relying mainly on thermally excited pions P3| . In the standard approximation to 
the NJL model any attempt to describe this behavior fails whereas it could be shown that 
as well in the l/iV c -expansion scheme as in the non-perturbative scheme HS] it can be 
reproduced. In the latter scheme also an examination of the chiral phase phase transition 
is possible @j| g|. 

In the first chapter we establish the different approximation schemes to the NJL model. 
To that end we begin by briefly reviewing the standard approximation, Hartree + RPA, 
and in which way quark and meson properties are described within that scheme. Af- 
terwards we present three different possibilities to extend the standard scheme, first an 
expansion up to next-to- leading order in 1/N C , then a "<I>-derivable method", and finally 
an one-meson-loop expansion of the effective action. The two latter schemes have non- 
perturbative character in contrast to the former scheme. For the l/iV c -expansion scheme 
and the one-meson-loop approximation, the schemes which we will later apply to study 
properties of quarks and mesons, we will explicitly show the consistency with the Gold- 
stone theorem and the Gell-Mann-Oakes-Renner relation as well as the transversality of 
the polarization function in the vector channel. The chapter will close with a discussion 
of a particular approximation to the two extended schemes which points out the relation 
of our calculations to hadronic models. 

The numerical results at zero temperature will be presented in Chapter ^[ We begin 
by explaining the regularization procedure before we come to investigate the influence 
of mesonic fluctuations on the quark condensate and pion properties. Then we try to 
determine a set of model parameters which is consistent with quantities as well in the 
pion sector as in the p-meson sector. An examination of further quantities related to the 
p-meson and a comparison with hadronic calculations follows. Finally we briefly discuss 
properties of the a-meson. 

Chapter |3| will be devoted to a discussion of results at nonzero temperature. We will 
begin by establishing the formalism for calculations at nonzero temperature (and density) 
which will first be applied to illustrate some basic results in the standard approxima- 
tion scheme. Then we will discuss results for the quark condensate obtained in the two 
extended approximation schemes at nonzero temperature. 

Finally we will summarize our main results and present some possible objectives. 
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Chapter 1 

The Nambu— Jona-Lasinio model 



1.1 Basics of the model 

The Nambu- Jona-Lasinio model was originally introduced by Nambu and Jona- 



Lasinio [|13[] in 1961 to describe an effective nucleon-nucleon interaction. In its original 
form it contained an isospin-doublet built of proton and neutron. Nowadays it is rein- 
terpreted, incorporating quarks instead of nucleons. In analogy to the nucleonic isospin- 
doublet one now deals with up- and down quarks. If one extends the model to flavor- SU (3) 



also strange quarks can be described [461 H7L Hq, 49 . We will restrict our investigations 



to the model with two flavors. Thus, we consider the following Lagrangian: 

L = - m )if) + g s [(^) 2 + (^7s^) 2 ] - 9v [(^rif)) 2 + (^7 M 75^) 2 ] ■ (1.1) 

if) is a quark field with Nf = 2 flavors and N c = 3 colors. g s and g v are coupling constants 
of the dimension energy" 2 . This Lagrangian exhibits the same global symmetries as QCD 
if in flavor space only the S*?7(2)-isospin-sector is taken into account. Via the Noether 
theorem these global symmetries are related to conserved currents. The global symmetries 
of the above Lagrangian are: 

• Invariance under Uy(l) transformations 

if) — ► exp[— ia)if> . 

The corresponding current is: 

This symmetry is related to conservation of baryon number. 

• Invariance under S'L r y(2) transformations 

if) — ► exp[— ir ■ 0/2]ip . 

This can be seen by recalling that ipi^Tif) is a vector in isospin space and the 
transformation corresponds to a rotation in that space. The scalar product of a 
vector with itself then obviously transforms like a scalar. The same argumentation 
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leads to the conclusion that as well {ip^rij)) 2 and {ip^^^^rip) 2 transform like scalars. 
It is obvious that also (ipip) 2 transforms like a scalar. The corresponding current is: 

Ji = h»T k i> . 

• In the chiral limit, i.e. for mo = 0, in addition invariance under SUa(2) transforma- 
tions 

if) — > exp[— 275T • 9/2]if> 

is found. One can most easily convince oneself of this invariance by considering the 
following transformation properties: 

(ipip) — > (4>ip) cos 9 — (ipij^T ■ 9ip) sin 9 , 
(ipijsTip) — ► (ipij^rip) + (ipip)9 sm 9 — (ipijsT ■ 9ip)9(l — cos 9) . 



Here 9 denotes a unit vector in ^-direction, |. The corresponding current is: 

Jl = H^l5T k 1p. 

We can conclude that in the chiral limit the Lagrangian, Eq. (|1 . 1|) , is invariant under 
global SUv(2) < £>SUa('2) < £}Uv (1) transformations. The invariance under a SUy (2)® SUa{2) 
transformation is called "chiral symmetry". The missing of degenerate "chiral partners" 
in the hadronic spectrum suggests that chiral symmetry is spontaneously broken in the 
QCD vacuum. This has important consequences, among others this requires the existence 
of massless Goldstone bosons. Because of the small mass of the pions compared with 
other hadronic scales these are interpreted as Goldstone bosons. The nonzero mass is 
thereby attributed to the explicit breaking of chiral symmetry due to the nonzero current 
masses of up- and down-quarks. The NJL model has been used by many authors (see e.g. 



the reviews [[[(], [L7|, |18|| ) to study the spontaneous breakdown of chiral symmetry in the 
vacuum and its restoration at nonzero temperature and density. 

However, a principal difference to QCD is the treatment of color degrees of freedom: 
In QCD those are related to an invariance under local transformations of the gauge group 
SU(N C ), whereas in the present model iV c only counts the number of identical copies 
of quark fields. Usually one chooses the coupling constants g s and g v to be of order 
1/N C 3D | to reproduce the large- iV c behavior of QCD. Treating iV c as an expansion 



parameter will enable us to construct a scheme to describe mesons which is consistent 
with the requirements of chiral symmetry. All calculations will, however, be performed 
with the physical value iV c = 3. 



1.2 Standard approximation: Hartree-(Fock)+RPA 



Detailed analyses of the results one obtains in the NJL model using the standard approx- 
imation scheme can be found among others in the reviews [16[ 117, Before we will 



come to discuss possible extensions of the standard approximation scheme in the following 
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sections we will briefly review the approximation usually applied to the NJL model. On 
the quark level this corresponds to a (Bogoliubov-) Hartree approximation. The local 
four-fermion interaction within the NJL model allows to cast exchange terms, i.e. Fock 
terms, in the form of direct terms with the help of a Fierz transformation. Therfore a 
Hartree- Fock approximation is similar to a Hartree one with a redefined coupling constant. 
Since the contributions from Fock terms are suppressed by one order in 1/N C as compared 
with the Hartree terms, we will restrict ourselves here to the Hartree approximation. 
The Dyson-equation for the quark-propagator in Hartree approximation is shown di- 



agrammatically in Fig. |1 . 1| . Note that we have renounced here to show the arrows on the 
quark propagators which would be standard convention for fermion lines. However, since 
later on most of the diagrams exist in two versions which only differ by the orientation 
of the quark loop, it seems convenient to display only one version which has then to be 
understood to contain all orientations which lead to topological different diagrams. 

Via selfconsistently solving the Dyson-equation the quarks acquire a momentum in- 
dependent self-energy which leads to a nonzero "constituent" quark mass m 



m = itlq 

M 



Z H = m + J2 rV2i<7M / -0^ Tr [ T M S(p) }. (1.2) 



S(p) = (p' — m + ie)' 1 denotes here a quark propagator in Hartree approximation and the 
symbol "Tr" stands for a trace over color-, flavor- and Dirac-indices. In principle the sum 
over M contains all interaction channels, i.e. M = a, it, p, a% with 

Y a = 1, r^=ry 5 T a , 
r;° = 7 'V a and = ^ 5 r a . (1.3) 



V 



The corresponding coupling constants are gu = g s for M = a or M = n and qm = 9 
for M = p or M = a\. It can easily be seen, however, that in vacuum only the scalar 
channel, M — a, contributes. Consequently the quark self-energy is proportional to the 
identity matrix. One obtains 

/(It) ffi 
(27r) 4 p 2 — m l + ie 

Since g s is of order 1/N C , the constituent quark mass m and consequently the quark- 
propagator are of 0(1) in Hartree approximation. This corresponds to the leading order 
in a l/iV c -expansion for the quark self-energy. 

In the chiral limit, i.e. for vanishing current quark mass mo, it is obvious that always a 
"trivial" solution of Eq. (|1.4|) with m = exists. If the scalar coupling constant g s exceeds 




Figure 1.1: The Dyson- equation for the quark-propagator in Hartree- approximation 
(solid). Dashed lines denote propagators of bare quarks. 
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a certain critical value, in addition there exists a solution of Eq. (O) with a nonzero 
constituent quark mass m. Because of the resulting "gap" in the spectrum, Eq. (|1.4j) 
is called, in imitation of BCS theory for superconductors, "gap equation". A nonzero 
constituent quark mass reflects the spontaneously broken symmetry of the underlying 
ground state. This can e.g. be seen if one considers the relation between the constituent 
quark mass and the order parameter of chiral symmetry, the quark condensate. Generally 
the quark condensate is given by 

m = -i J -0-,Tr S(p). (1.5) 

In the Hartree approximation it is directly related to the constituent quark mass, 

= -rn^a . (i.6) 

We have introduced here the superscript (0) to indicate that we deal with a quantity in 
Hartree approximation. 

Mesons can be described via a Bethe-Salpeter equation. The leading order in 1/N C is 
here given by a so-called "Random Phase approximation" (RPA) without exchange terms. 
Diagrammatically this Bethe-Salpeter equation is displayed in Fig. |1.2| . 

The propagators of the mesons can be extracted from the quark-antiquark scattering 
matrix T, which satisfies the following equation 



■ M,ijkl 



(q) — K M)i j k [ + K Mi j ab J bcda (q)T Mcdk i(q) , (1.7) 



with Jbcda(q) = — i I -^iSi, c {p + q/2)Sd a (p — q/2). The indices a to d and i to I are 
multi-indices, indicating the components of the various quantities in color-, flavor- and 
Dirac-space. Km denotes the scattering kernel, which can be written in the following way: 

KM,ijkl = 2gM^M,ij^M,kl ■ (1-8) 



The vertices Tm, M = a, tt, p, a± have already been defined in Eq. (|1.3|) together with the 
corresponding coupling constants gM- If we make the ansatz 



■ M,ijkl 



(q) = -D M (q)T M>ij T M ,ki , (1-9) 



for the T-matrix, Eq. ( |1.7| ) can be transformed into an equation for the meson propagators 
Dm, 

D M (q) = -2g M + 2g M U M (q)D M (q) , (1.10) 



with the polarization functions I1m- These consist of a quark- antiquark-loop 

U M {q) = F M ,abJdabc(q)T M ,cd = -i J T^jTi Tr [ Tm iS (p + |) V M iS{p - |) ] . (1.11) 



Here again Tr denotes a trace over color-, flavor- and Dirac- indices. In the scalar and the 
pseudoscalar channel, i.e. for the cr-meson and the pion, we obtain 



12 






+ 



Figure 1.2: The Bethe-Salpeter equation for the meson-propagators in RPA (double line). 
The solid lines indicate quark propagators. 

D ' (q) = 1 - 2 gM ) ' 

D-to) . 0,W^ = 1 _- 2 ^ fa) ^. (1.12) 

Here a and 6 are isospin indices. We have used the following notation: n^f(g) = n^g) 5 a6 . 

In the vector channel we proceed in a similar way. The only difference is the the more 
complicated Lorentz structure of the polarization function and the resulting propagator. 
Since the polarization function can depend only on one four-momentum, in our notation 
q, its tensor structure is generally given by 

n^ afe (g) = W7i(<? 2 ) + ^r/ 2 (g 2 )) , (1.13) 

with two scalar functions f\ and f 2 . This can be separated into a tranverse and a longi- 
tudinal part using the corresponding projectors, 

T l»> = + JJW = . (1.14) 

qi qi 

However, because of vector current conservation, the polarization function in the vector 
channel has to be transverse, i.e. 

q,U^ a \q)=q u U^ a \q) = 0. (1.15) 

The tranverse structure, together with the assumption of Lorentz invariance, not only 
simplifies the determination of the tensor structure of the corresponding propagator but 
also the evaluation of the polarization function itself: With the ansatz 

U^ ab (q) = T^(q)5 ab U p (q) (1.16) 

we only need to evaluate the (Lorentz) scalar function 

nf(?) = -\g^ a \ q ) . (1.17) 

Finally one arrives at the following expression for the propagator of the p-Meson 

Dp ab (q) = (D p (q) T^ + 2g v IT) 5 ab = ( ^"^ (g) ^ + 2g v IT) S ab . (1.18) 
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In section L5 we will explicitly show that the polarization function in the vector channel is 
indeed transverse, provided that a suitable regularization procedure is used. In section ^J] 
we will discuss a further consequence of vector current conservation, namely that H p {q) 
should vanish for q 2 = 0, in the context of the regularization procedure we will apply. 

In the same way the propagator of the ai-meson can be obtained from the transverse 
part of the axial polarization function 11^ (q)Sab = — \T pu Tl l ^' ab (q). As discussed e.g. 



in Ref. |49]], 11^ in addition contains a longitudinal part which contributes to the pion 
propagator together with the pseudoscalar polarization function U^ b and the mixed ones, 
which contain one pseudoscalar and one axial vertex. There is no conceptual difficulty in 
dealing with this so-called "7r-ai-mixing" . We will neglect it in the following discussion, 
and only consider the pseudoscalar part of the pion propagator, in order to keep the 
formalism as simple as possible. The pion propagator including 7r-ai-mixing is discussed 
in App. [IT]. 

It follows from Eqs. ( |1.11|) to ( P-.18| ) that the functions £)jv/(g) are of order 1/N C . The 
explicit form of these functions can be found in App. [B]. We will call them "propagators" 
although strictly speaking they have to be interpreted as the product of a meson propa- 
gator with a meson-quark coupling constant squared. The latter is given by the inverse 
of the residue of the function Z?m(?) at the pole, whereas the mass of the corresponding 
meson is determined by the location of the pole, 

n-lf \ I n -2(°) dIl M (q) 

m w\ q 2 =m m = u , 9 Mqq = dq2 \ q 2= m \f ■ C 1 - 19 ; 

The superscript (0) here serves to indicate that rri^p and gff qq are RPA quantities. In 
a 1 /iV c -expansion these quantities are of leading order. One can easily convince oneself 
that they are of order 0(1) and l/y/N~ c , respectively. Dealing with 7r-ai-mixing, we have 
only to keep in mind that in addition to a pseudoscalar pion-quark coupling constant an 
axial one exists. This will be discussed in App. |C.1| . 



1.3 Extensions of the standard approximation 

The Hartree + RPA (only direct contributions) scheme, which was discussed in the pre- 
vious section is consistent with chiral symmetry. This means, for instance, that a spon- 
taneously broken symmetry enforces the existence of massless Goldstone bosons. Within 
that approximation scheme also the validity of the Goldberger-Treiman and the Gell- 
Mann-Oakes-Renner relation, which determines the behavior of the pion mass if the sym- 



metry is explicitly broken due to a small current quark mass mo, can be shown [110, O, PJI] . 



We will present a proof for these relations in section |1.4.1| . Since chiral symmetry is one of 



the main features of the NJL model, any approximation should be consistent with chiral 
symmetry. Within this section we will discuss different approximation schemes which 
have in common that on the one hand they fulfill the requirements of chiral symmetry 
and on the other hand they go beyond the standard Hartree + RPA scheme. We will 
begin with a strict l/iV c -expansion up to next-to-leading order, proceed by discussing the 
"^-derivable method" and then consider a one-loop expansion of the effective action. The 
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Figure 1.3: The 1/N C - correction terms 5£W (left) and (right) to the quark self- 

energy. 

latter corresponds to a selfconsistent extension of the gap equation including a meson-loop 
term. These methods enable us to include mesonic fluctuations in our investigations. 

1.3.1 l/iV c -expansion 

Within this section we will consider the quark self-energy and the polarization functions 
of the mesons in a l/A^-expansion up to next-to- leading order. To determine the corre- 
sponding correction terms let us first recall the l/A^-counting rules for the leading-order 
quantities: Quark propagators are of order unity whereas the meson propagators are of 
order 1/N C . In addition one has to keep in mind that, because of the trace in color space, 
every quark loop contributes a factor N c and every four-fermion vertex a factor 1/N C due 
to the order of the coupling constant. 

For the quark self-energy we find two correction terms, 

6E(p) = 5£ (a) +5E (b) (p) . (1.20) 



These terms are graphically shown in Fig. |1.3| . A solid line corresponds here to a quark 
propagator in Hartree approximation (0(1)) and a double line to an RPA meson prop- 
agator (0(1/N C )). With the help of the above mentioned counting rules it is then easy 
to convince oneself that these are indeed the only correction terms to the quark self- 
energy in next-to-leading order. Any further contribution inevitably contains in addition 
at least one quark loop and two coupling constants or one meson propagator, which leads 
altogether to an additional factor 1/N C . 

The correction terms to the quark self-energy determine those to the quark propagator 
to order 1/N C , 

6S(p) = S(p)SE(p)S(p). (1.21) 

S(p) here denotes the quark propagator in Hartree approximation. In leading order the 
self-energy has been iterated in order to obtain the corresponding propagator. This has 
been done via the selfconsistent solution of the gap equation, Eq. ( |1.2|) . However, in next- 
to-leading order an iteration of the self-energy would lead to terms of arbitrary order in 
1/N C in the propagator. Following Eq. (|1.5| ) the correction to the quark condensate is 
given by 

5(^) = -i J^TrSS(p). (1.22) 
We find four correction terms, 5U$ to SU.^ , to the meson polarization functions. 



These are displayed in Fig. |1.4j together with the leading order contribution. Solid lines 
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denote again quark propagators, double lines RPA meson propagators. This leads to the 
following form of the complete polarization function: 



U M (q) = U M (q) + SU^(q). (1.23) 



k=a,b,c,d 



The correction terms, which contain either one RPA propagator and one quark loop or 
two RPA propagators and two quark loops, are of order unity, whereas the leading-order 
term is of order iV c . 

The meson propagators are again, similar to the leading order, given as solution of a 
Bethe-Salpeter equation. The only difference to the propagators in RPA consists of taking 
into account the polarization functions up to next-to-leading order. Thus we arrive at the 
following expressions for the meson propagators, analogously to Eqs. Q1.12Q and ( |1.18|) , 



1 - 2g M U M (q) 

We should remark that these meson propagators contain arbitrary orders in 1/N C al- 
though we have determined the polarization functions by a strict expansion in powers of 
1/N C . This is due to the iteration of the product of the polarization functions with the 
corresponding coupling constant. 

This remark similarly concerns the meson masses. These are defined, in analogy to 
Eq. ( |1.19|) , as the location of the pole of the propagator, 



D- A l(q)\ q , =mli = . (1.25) 

Because of the implicit definition, tum also contains terms of arbritrary orders in 1/N C . 
This definition is consistent with the Goldstone theorem, but in the context of the Gell- 
Mann-Oakes-Renner relation we will encounter difficulties caused by higher-order (be- 
yond next-to-leading order) contributions to the pion mass. This point will be discussed 
in detail in section |1.4.3| . 

For the evaluation of the various contributions it is convenient to introduce effective 
meson-meson vertices which consist of quark loops. We need two different types of meson- 
meson vertices, a three-meson vertex, containing a quark triangle, shown on the l.h.s. of 



Fig. 1.5 and a four-meson vertex, containing a quark box, shown on the r.h.s. of Fig. 1.5 



For external mesons Mi, M2 and M3 the quark triangle has the following form: 

f d 4 k 

-«r Ml ,M 2 ,M 3 (g,p) = -J j^—^{Tr[r Ml iS(k)r M2 iS(k-p)r M3 iS(k + q)] 

+Tr [T Ml iS(k - q)r M3 iS(k+p)T M2 iS(k)]} (1.26) 



The operators T M have already been defined below Eq. (|1.2j ). The above expression 
already contains a sum over both possible orientations of the quark loop. 
The four-meson vertices can be written in the following way: 
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Figure 1.4: Contributions to the meson polarization functions in leading (RPA) and next- 
to-leading order in 1/N C . 



-i^Mi,M2,M3, M 4 iPU Vli P3 

d A k 



(2tt)< 



Tr [T Ml iS(k)T M2 iS(k - p 2 )T M3 iS(k - p 2 - p 3 )T M JS(k + p x )} 



+ Tr[T Ml iS(k - Pi)T M JS{k+p 2 + p 3 )T M3 iS(k + p 2 )T M jS(k)]) . (1.27) 



Here again we have summed over both possible orientations of the quark loop. 

With the help of the above definitions the l/iV c -correction terms to the quark self- 
energy as well as to the meson polarization functions can be written in a relatively compact 
form. For the momentum independent contribution ST,^ to the quark self-energy we 
obtain 



M - (2^) 

Here we have introduced the constant A 



1.28) 



2J (2tt) 4 m 



;i.29) 



which we will need later on for instance for the evaluation of diagram £11^. The factor 
1/2 in the definition of A is a symmetry factor, which is necessary in order to avoid double 
counting which otherwise would arise from the sum over both possible orientations of the 
quark loop in the definition of quark triangle, Eq. ( |1.26|) . A consists of a quark triangle 
coupling to a meson loop and an external scalar coupling. In principle we should have 
allowed also for other external couplings, for instance a pseudoscalar one, and sum over 
all these couplings and the corresponding meson propagators in Eq. (1.28). It turns out, 




P,M 2 \p3,m\- Ps,M 3 
pf = -p- q, M 3 /■ p u Mi i I \ Pi, M 4 



Figure 1.5: Quark triangle 

-«r Ml )Ma Mi Mi (pi ,p 2 ,p 3 ) ( right) . 



-^rM 1 ,M 2 ,M 3 (Q , ,p) (left) and quark box 
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however, that only the scalar coupling leads to a nonvanishing contribution to A because 
FmmmM -P) = for M x ^ a (cf. App. 0). 

For the momentum dependent correction term 5Y>^ to the quark self-energy we arrive 
at the following expression: 

5^ b \k)=iY: j ^r i D M {j>)V M S{k-p)Y M . (1.30) 



M 



The l/iV c -correction to the quark condensate is closely related to the quark self-energy 
(cf. Eqs. ( |1 . 5p and (|1.22|) ). Inserting the expressions for <5£( a ) and 5£< 6 ) into Eq. flPID , 
and performing the trace in Eq. ( p. .22 ) we obtain 



6(W>) = -AU a (0)D a (0) + A . (1.31) 

This expression can be further simplified with the help of the definition of the RPA 
propagator D a in terms of the corresponding polarization function ILj-, cf. Eq. ( |1.12| ). We 
finally obtain for the l/iV c -correction to the quark condensate 

= -m± . (1 . 32 ) 

The evaluation of the various contributions to the polarization functions give 

23 ( ) M\M% 

^ n S(g) = ~i I -ttAt J2 V m,Mi,Mi,m(^P^-p) d m 1 (p), 

J y A7r ) Mi 

^M(Q) = ~\ [ T^U Z) r M,Mi,M,Mi(liP,-q) D M 1 (p) , 

2J ( 2n r Ml 

Su m(q) = 77 r M,M, CT (?, -q) D*(0) I 777-T7 X) V M 1 ,Mx^-P) D Mi{p) , 

2 J ( 2n r Ml 

= -iT M)Mja (q,-q)D a (0)A. (1.33) 

The symmetry factor 1/2 in 511^ and 811$ has been introduced for the same reason as 
the factor 1/2 in Eq. flOg) . In <ffl£ } it is necessary in order to account for the fact that 



the exchange of Mi and M2 leads to the same diagram. 

For the further evaluation of Eqs. ( |1.28| ) to ( |1.33| ) we will proceed in two steps. First 



the quark loops for the RPA meson propagators and the effective three- and four-meson 
vertices will be calculated and then we will evaluate the remaining meson loops. 

In principle the sum in Eqs. ( [L.28} ) to (|1.33| ) is over all possible intermediate states, i.e. 



cr-,7r-,p— and ai-mesons. However, pions, the lightest particles in the system, are expected 
to yield the dominant contribution to most applications. For instance, the behavior of the 
quark condensate at low temperatures is mainly driven by thermally excited pions. The 
p-meson spectral function is also primarily determined by a two-pion intermediate state. 
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Other possible contributions to SU^ a \ like nai-, pa-, pp- or a±ai- intermediate states are 
expected to be much less important since the corresponding decay channels open far above 
the p-meson mass. In principle we have to keep in mind that the physical a± as well as 
the p are broad resonances which principally couple to three- and two-pion intermediate 
states, respectively. Therefore the physical threshold for the decay of the p into Tca\ or 
pp lies at Am-x ~ 560 MeV, i.e. below the p-meson mass of about 770 MeV. But of course 
considerable contributions will arise only above 2m p ~ 1.6 GeV and m ai -{-m^ ~ 1.4 GeV. 
Besides, experimentally the contribution of these intermediate states is negligible anyway, 
one finds that the branching ratio of p — > mi is about 100% |50] . 



From a phenomenological point of view it therefore seems well justified to consider only 
pionic intermediate states. However, consistency with chiral symmetry requires that also 
scalar states are considered, whereas vector and axial intermediate states can be neglected 
without any difficulties. Since the evaluation of the various contributions can be simplified 
considerably by neglecting vector and axial intermediate states we will exclusively take 
scalar and pseudoscalar intermediate states into account for most of the calculations. 
To describe a p-meson we have certainly to allow for an external vector coupling in the 
diagrams shown in Fig. |1 .4j . 

Another point which complicates the treatment of vector and axial intermediate states 
is that they cause strong divergencies in the correction terms to the polarization functions. 
On the one hand this has certainly to be attributed to the non-renormalizability of the 
NJL model, on the other hand it is related to a rather fundamental problem concerning the 
regularization of the RPA vector and axial polarization function. This will be discussed 
in more detail in section 2A and section 2.4.3| , respectively. 



1.3.2 ^-derivable theory 

The main disadvantage of the approximation discussed in the previous section is that 
the 1 /iV c -correction terms are treated perturbatively. For instance, though we have con- 
sidered corrections to the quark self-energy, these have not been taken into account in 
the gap equation. Consequently all quark propagators contributing to the l/iV c -corrected 
polarization functions are determined in Hartree approximation. As long as the 1/N C - 
correction terms are small compared with the leading-order contributions this approach 
seems reasonable. But a description of the chiral phase transition will certainly not be 
possible within that scheme. 

An approach which treats the mesons selfconsistently is highly desirable, since e.g. the 
thresholds for the decay p — > nai would then be located at Am^, but this is a very difficult 
task. For this reason we will be content with a scheme which enables us at least to treat 
the quarks selfconsistently. In the literature two functional methods can be found which, 
in principle, allow to construct a self consistent scheme which is consistent with chiral 
symmetry. One is the so-called "^-derivable" -method E7\, PSI, the other an expansion of 



the "effective action" j51], [52| . The latter will be discussed in the next section. We will 
begin here by illustrating the ^-derivable method on the Hartree + RPA scheme which 
emerges as the first approximation. Subsequently we will go one step further which will 
generate a non-local contribution to the gap equation. 

The starting point of a ^-derivable theory is a functional $, which in principle con- 
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tains all two-particle-irreducible skeleton diagrams which can be constructed from the 
interactions and fields of the underlying Lagrangian. Two-particle-irreducible diagrams 
are those which cannot be split into parts by cutting two lines. In the diagrams all 
Greens functions have to be full ones. The corresponding self-energies can be obtained 
by functional derivatives of $ with respect to the (full) Greens functions. Generally a 
symmetry conserving approximation is generated by taking any subset of diagrams [ j38fl . 
The thermodynamic potential per volume, Q, which is identical to the energy density of 
the system in vacuum, can be expressed as follows with the help of the functional $, 

n = -iTr ln(i S' 1 ) + Tr(S iS) + , (1.34) 

where S^ 1 = Sq 1 + X represents the full Greens function, i.e. in our case the full quark 
propagator. The symbol Tr denotes here a trace over internal degrees of freedom, such 
as color-, spin- or isospin degrees of freedom, as well as an integral over momentum or 
coordinate space. For instance, in vacuum we have 

d 4 p 



Tr(S5) = / * Tt(E(p)£(p)) 



(27T) 

An expression for the self-energy can be derived by requiring stationarity of Q with respect 
to variations of the full propagator. One finds 

(1.35) 



5 (iS) 

Evidently this means that the self-energy is given by the functional derivative of Q with 
respect to the full propagator. The scattering kernel for the Bethe-Salpeter equation 
is in turn given by the functional derivative of the self-energy with respect to the full 
propagator, 

5iS~ 6 (iS)8 (iS) • 1 6b) 

Let us now begin considering the simplest possible subset of two-particle-irreducible 
diagrams. In the NJL model this is given by the "glasses" , displayed in Fig. |1.6| . Again, 
solid lines represent quark propagators. The wavy lines have been introduced to visualize 
the direction of the interaction at the four-point vertices. To that approximation $ is 




Figure 1.6: The "glasses", simplest contribution to the functional $. 
explicitly given by 

0(0)^(0] = -^ 9 m(-J ^Mr Mi S^\p))) 2 . (1.37) 



M 
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The superscript (0) again denotes quantities to first approximation. Below we will demon- 
strate that the above expression for the functional $ generates the Hartree + RPA scheme 
for quarks and mesons, respectively. 

According to Eq. ( |1.35|) the quark self-energy is a priori momentum dependent and a 
matrix in color-, flavor- and Dirac-space. Performing the functional derivative (i and j 
are multi-indices, referring to color-, flavor- and Dirac-space), 



J 1 6 (tS> j '(x 1 ,x 2 )) 



= %i + g 2 ) £ 2^r Mji / A^Tr(T M S^(p)) , (1.38) 

M V / 

we find, however, that the self-energy is at least momentum independent. Remembering 
that only the scalar channel leads to a nonvanishing contribution we conclude that it 
is in addition proportional to unity in all spaces. After a comparison of the resulting 
expression for the self-energy with the gap equation, Eq. (|TT2"|), it becomes obvious that 
is identical to the Hartree self-energy derived in section T72. 
For the scattering kernel we obtain 



K ij?ki(Qii ?2, 33,24) = / d A x x J d 4 x 2 J d A x 3 J d 4 x 4 e 



5 {iSf? (x 3 ,x A )) 

£ 29M^M,ji^M,ikS(qi + <? 2 + <?3 + Qi) ■ (1-39) 



M 



As can be seen from a comparison with Eq. ( |1.8|) this expression meets our expectations: 
It is indeed the scattering kernel for the T-matrix in RPA. We therefore conclude that the 
Hartree + RPA approximation scheme can be derived within the ^-derivable method, if 
the functional <E> is restricted to the "glasses" . 







Figure 1.7: The " ring sum", contribution to the functional $. 



A suggestive extension for <3> is to take also the "ring sum", visualized in Fig. |IT7 , 
into account. Evaluating the diagrams, "glasses" and ring sum, we arrive at the following 
expression for <3>: 



m = -T,9M(-J-^-Jr(T M tS(p)f 



M 

74 



II 7^{ ln (l-2^H CT (p)) + 31n(l-2^(p))} . (1.40) 
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The functions Um{p) were defined in Eq. fll.ll ). We have to keep in mind that S here 



denotes a full propagator and not a propagator in Hartree approximation as it originally 
appeared in the definition of 11^ in Eq. ( |1.11| ). Proceeding in the same way as above (cf. 



Eq. ( |1.38|) ) we obtain for the quark self-energy, q 2 ) = <K<7i + 92)^^(5 



Z ij (q 1 ) = X H)ij + 6Xlf(q 1 ) . (1.41) 

S'E^(q) was defined in Eq. ( |1.30p . The same remark of caution as above concerning 
the quark propagator S is in order here: The expressions on the r.h.s. of Eq. (|L4l|) 
have to be evaluated with the full propagator. The upper panel of Fig. displays the 
corresponding gap equation. Here solid lines stand for quark propagators as they emerge 
from the selfconsistent solution of Eq. ( |1.41| ). 

To determine the scattering kernel for the Bethe-Salpeter equation which describes 
the mesons, the analogous procedure leading to Eq. ( |1.39| ) is applied. One finds 



r d 4 k 

+2 i J T^Ti D M{k + qs)D M (k - q A ){T M iS{k - Ql - q±)T m)^ M iS{k)T M ) kl ) . (1.42) 



Kij,ki(qi, <?2, <?3, 94) = $(qi + <?2 + <?3 + qA){Yl ^9M^M^M,ik + D M (q 1 + qi)T M ^T M ,jk 

M 

d A k 

The corresponding Bethe-Salpeter equation for the quark-antiquark T-matrix is displayed 
in the lower part of Fig. |1.8| . In Ref. [|^] this representation of the T-matrix was generated 
graphically by coupling an external current to each quark line in the corresponding gap 
equation, which is shown in the upper part of Fig. |1.8| . Moreover, with the help of axial 
Ward identities one can show generally the existence of a pole in T for (qi + q 2 ) 2 = 0. 

Because of the non-trivial momentum dependence of the self-energy it is a very difficult 
task to obtain an explicit solution of the gap equation, Eq. (|1.41| ). Therefore we will put 



this scheme aside and concentrate in explicit calculations on the expansion of the effective 
action which leads to a local contribution to the gap equation and is much easier to handle. 
This will be discussed in the next section. 



1.3.3 One- meson-loop expansion of the effective action 

The scheme we will derive within this section was first discussed by Dmitrasinovic et 



al. P9[ . The authors started from a gap equation which in addition to the Hartree self- 
energy contains a local meson-loop contribution. They found a consistent scheme to 
describe mesons and proved various relations following from chiral symmetry explicitly , 
like e.g. the Goldstone theorem. The same scheme was later derived more systematically 
by Nikolov et al. |5D[ using the effective action formalism. Here we will follow the method 



of that reference. Let us begin with a derivation of an expression for the effective action. 
We will mainly follow the presentation of Ref. [p6[| , Chapter 16. The interested reader is 
referred to that reference for further details. The starting point is a theory with an action 
I[4>]- Throughout the derivation we will, for simplicity, only regard a theory with a scalar 
field 0, which couples to an external classical source j. The generating functional W[j] is 
then given in path integral representation by 

e iwy\ = f V(<P) exp (il[<j)\ + % f d A x<j){x)]{x)) (1.43) 
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Figure 1.8: (upper panel) Gap equation as it results from $ if the ring sum is included, 
(lower panel) corresponding quark- antiquark scattering matrix. 



The functional W\j] contains all connected vacuum diagrams, which can be constructed 
from the underlying theory in the presence of the source j. The vacuum expectation value 
of the fields <fi can be obtained from a derivative of W [j] with respect to the external source 

5W[j] 



<0(*)> 



5j(x) 



j=0 



;i.44) 



This expectation value can in general also be defined in the presence of the source, we 
will call it <£>j, 

8W[j] 



5j(x) 



;i.45) 



For a given ipj this equation can in principle be solved for J. Throughout the following 
derivation we will denote the current, which is related to ifj via Eq. ( |1.45| ), by J. A 
functional r[<^j], which no longer depends on j but on (pj, can be obtained from W[j] by 
performing a Legendre transform, 

T&j] = W[J] - J d 4 x Vj (x)J(x) . (1.46) 

This functional ^[fj] is called "effective action". It contains all one-particle-irreducible 
diagrams (in the presence of the source) of the theory |3~6|| . Similarly to two-particle- 
irreducible diagrams, one-particle-irreducible ones are those which cannot be split into 
parts by cutting one line. One important property of T can be derived from the derivative 
with respect to (pj, 



d 4 x 



5W[j] 8j(x) 



Sj(x) 8(pj{y) 



d 4 xipj(x) 



5J{x) 
6<p(y) 



-J(y) 



(1.47) 
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Hence, if the source J vanishes, only the fields ips at the stationary points of T[(p d 



:i.48) 



tp=(p*- 



are possible fields ipj. This means in particular that the vacuum expectation values of the 
field <p are given as the values of (fj at the stationary points of T, i.e. ip s = (0). From this 
point of view the field <pj behaves similarly to a classical field whose equations of motion 
can be generated from the stationarity condition of the action I[<p]. Thus the effective 
action can be regarded as some sort of "action with quantum fluctuations" . 

Propagators can be obtained from the second derivative of W[j] with respect to the 
source, 

6 2 W[j] 5<pj(x) 



A(x,y) = ($(x)$(y)) 



This can also be derived from the effective action, 



3=0 



5j(y) 



3=0 



5 2 r\ 



8<p(x)5ip(y) 



6j(x) 



3=0 



A(x,y)~ 



1.49) 



;i.5o) 



It is convenient to introduce another quantity, the effective potential. Let us as- 
sume that we wish to calculate the effective action T[tp] with constant, i.e. space-time- 
independent, fields (p(x) = (fQ. Then V[tp] will be proportional to the volume of space-time 
V4 = / d 4 x. Therefore we can write the effective action in the form 



r>] = -v 4 v(^ ) 



;i.5i) 



where V((po) is not a functional of ip but an ordinary function. This function is known 
as the effective potential jp3"| . The effective potential can be interpreted in terms of the 
energy density |3£|. V(ipo) can be interpreted as the minimum of the energy density for 
all states with expectation values (p for the fields 0. 

After having explained some general aspects of the effective action let us now come 
to the determination of an explicit expression for T[(p]. Certainly it will be impossible to 
obtain an exact result for the effective action in an interacting theory because the path 
integral in Eq. ( |1.46|) can in general not be carried out. This is only possible for constant, 
linear or quadratic terms. Thus, expanding the argument of the exponential in Eq. ( 1.43 ) 
around some fixed field (po(x) which satisfies 



61[<f>) 
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<P=<P0 



-J 



1.52) 



up to quadratic terms in the fields, we can evaluate the path integral and arrive at an 
expression for the effective action T[ipo\ up to "one-loop" (see [[36), Chapter 16), 



I*> M =iT M ^7Vln( i *g). 



11-53) 



The first term in Eq. (|1.53| ) is the "zero-loop" or "tree-level", the second one the "one- 
loop" contribution. Strictly speaking the field ip which arises from the Legendre transform 
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in Eq. ( 1.46Q not necessarily needs to coincide with ipo which has been introduced for the 



evaluation of the path integral. However, at tree-level it is obvious that both fields are 
identical. Besides, it can be shown that in a perturbative expansion this is true up to 



one- loop order [54 



It can be shown (see e.g. |f56fl ) that in most cases the effective action T[(p] has the 
same symmetries as the action I [(f)]. Especially symmetries which are related to linear 
transformations of the fields are conserved. This has an important consequence (see 
e.g. |36|, |55|): One can show that a spontaneously broken continuous symmetry then 



generates a pole at q 2 = in the two-point function A(q). This corresponds to the 
existence of a Goldstone boson. 

We are now in a position to apply the effective action formalism to the NJL model. 
A detailed discussion can be found in Ref. p5] . We will begin our presentation with 



the lowest-order contribution, i.e. the tree-level approximation which will turn out to be 
equivalent to Hartree approximation + RPA. 

In the remaining part of this section we will drop the vector and axial vector interaction 
and start from a Lagrangian which contains only scalar and pseudoscalar interaction 
terms, i. e. 

£ = $(0 - m )^ + g s [(^) 2 + (ipiKTip) 2 ] . (1.54) 
The generating functional of the system can be expressed in terms of the path integral 

Z[ v ,f]] = e iff M = J D(^)D(^)e i;[ ^* ,+! / rf4l( '* (l) ' ,w ^ w, '' W) , (1.55) 

with the action 

I$,ij>] = J d 4 xC(x) . (1.56) 

To further proceed it is convenient to "bosonize" the action. This is achieved by intro- 
ducing auxiliary hermitian fields <&' a ,a = {0,1,2,3}. These collective bosonic fields are 
chosen in such a way that the action becomes bilinear in the quark fields, i. e. 



Z[r], fj] = J V(i[))V(if;)V(& a ) exp {il(ifj, rf 4 x(^(x)^(x) + fj(x)if)(x)) 

~ J d 4 x(K + 2^r a ^) 2 } , (1.57) 

with T a = (1, 1757*). The quark fields can now be integrated out. After performing a shift 
of the auxiliary fields, $ a = $' a + (m , 0), one finally arrives at the bosonized form of the 
generating functional 

Z\^f]\ = J P($ a )exp{i/[$ ] +i J d 4 x J d A yf}(y)B-\y,x)rj(x)} (1.58) 

with the action 

/[$] = -i Tr lm'D - — j rf 4 x($ 2 - 2m $ + mjj) . (1.59) 
D is the Dirac operator 

D(x,y) = (i$-T a $ a )6(x-y) . (1.60) 
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This notation should not be confused with the meson propagator -Dm(<?) we introduced 
earlier. The first term in Eq. ( |1.59|) arises from the "fermion determinant" . The symbol 



Tr is therefore to be understood as a functional trace and a trace over internal degrees of 
freedom such as flavor, color and spin. This is the quark-loop contribution. Expanding 
the fields around some constant field $ a = (M, 0) + A$ a we can rewrite the Dirac operator 
as follows 

D(x,y) = W-M)6(x-y)-r a AQ tt 6{x-y) 

= D (x,y)-V, (1.61) 

where V is of first order in the fluctuations A$ a . Do corresponds to a quark propagator 
with "mass" M, it can be written as 

Do(x,v) = J -^j^Sip) ■ (1.62) 

The logarithm in the quark-loop term to the action can then be expressed in powers of 
the fluctuations A$ a , 

ln(iD) = ln(iDo) + ln(l - D ~ V) 

= ln(iD ) - £ . (L63) 

n 

This expression will be helpful when we determine the derivatives of the action / with 
respect to the fields $. In which way are the expectation values of the auxiliary fields $ 
related to physical quantities we are interested in? This can most directly be seen for the 
quark condensate: It can be expressed via the expectation value of $o as 

dW 1 

<« = S^ = -2^«*°>- m °>' (L64) 

From Eq. ( |1.50| ) we can infer that the meson propagators can be obtained as second-order 
derivatives from the effective action. In "zero-loop" approximation the effective action is 
simply given by, 

r(°)[$] = /[$] , (i.65) 

where we have dropped the fermionic source terms r\ and fj for simplicity. To determine 
the expectation values of the fields <3> we first have to build the functional derivative with 



respect to the fluctuations A$. To that end only the first order term in Eq. ( |1.63|) has to 
be considered. We arrive at the following expression 

Since the trace in the first term on the r.h.s. vanishes for b ^ there is no contribution 
with b ^ 0. Comparing the first term on the r.h.s. of the above expression with the 
definition of the Hartree self-energy we conclude that the expectation value of the fields 



26 



$ in the present approximation can be written in the form (M, 0), where M is given by 
the solution of the following equation: 



M-mo-E H (M) = , (1.67) 



with Y,h as defined in Eq. ( |1.2|) . If we identify the expectation value of the zeroth compo- 
nent of the field <3>, i.e. M, with the constituent quark mass in Hartree approximation, 
this equation is identical with the gap equation (Eq. ( |1.2| )). We have to emphasize here 
that the interpretation of ($o) = M as a constituent quark mass is not clear a priori. 
In principle the constituent quark mass should be determined from the pole of the quark 
propagator. However, in Hartree approximation this pole coincides with ($o)- 

If we evaluate the second-order derivative of the effective action 1^°) at the stationary 
point, we obtain the following result 



(1.68) 

where we have already exploited the fact that in our case the polarization function IIm 
is diagonal in flavor space. Thus, comparing this result with the definition of the inverse 
meson propagators in Eq. ( |1.50| ), we conclude that we in this way exactly recover the 
inverse meson propagators in RPA, cf. Eq. ( |1.12| ). This allows us to draw the conclusion 
that the effective action formalism in "zero-loop" approximation yields the same results 
as the Hartree approximation + RPA. 

Extending the effective action to "one-loop", starting from a bosonized version of 
the NJL model, means that we take into account mesonic fluctuations. The effective 
action is then given by Eq. ( |1.53| ). The second term in Eq. ( |1.53| ) contains the mesonic 
fluctuations. As discussed in Ref. |5l| the method is only meaningful if the second-order 
functional derivative which enters into this term is positive definite. Otherwise severe 
problems arise due to an ill-defined logarithm, which would then be complex. We will 
come back to this point in section |2.3|. 



Although we introduced this approximation rather generally as "one-loop" approxi- 



mation we will throughout this work adopt the notation of Ref. []3D]| and call this approx- 
imation scheme the "one-meson-loop approximation" (MLA). This name is motivated by 
the fact that one-loop contribution in fact contains mesonic fluctuations due to the pre- 
ceding bosonization procedure. In addition we have to mention that the "quark-loop" 
approximation described above is not really a tree-level approximation which would be 
suggested by the introduction of r*- -* in Eq. ( |1.53| ). In order to arrive at the bosonized 
form of the action, we have integrated out the quark fields, i.e. the bosonized action in 
principle contains all quark loops. 

For the following derivation of explicit expressions in one-meson-loop approximation 
it is convenient to first establish relations between the effective meson-vertices introduced 



in section 1.3.1 and the third- and fourth-order derivative of the action /[$], 
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J d i x 1 d i x 2 d i x^- i{qi - Xl+q2X2+qTX3) 



<5 3 /[$] 



5($ a (xi))5($ 6 (x2))5($ c (x 3 )) 



A$=0 



-i(qi-x\+q2-X2+q3-X3+qi-X4.) _ 



= S(qi + q 2 + g3)r o ,6,c(gi, ^2) 
5 4 /[$] 



A$=0 



5($ a (x 1 ))5($ b (x 2 ))5($ c (x3))5($ d (x 4 )) 
= -iS{q t + q 2 + qs + qi) a ,b,c,d{qi, £2, 93) + r a , MiC (gi, g 2 , + r 0)CiM (gi, g 3) g 2 ))(l-69) 

These relations enable us to straightforwardly derive an equation for the stationary point 
of the effective action in Eq. ( |1.53| ). The following "gap equation" holds |30f : 

($0) -mo- £«$„» = ($0) -mo- - <*£«$ » = . (1.70) 

Here is the Hartree contribution to ($0) as defined in Eq. (|1.2|). In the upper part of 



Fig. |1.9| this gap equation is visualized graphically. The cross indicates the constant term 
($0) — m 0i the second term corresponds to the Hartree contribution and the correction 
term <5E corresponds to the third diagram. It consists of a quark loop dressed by an 
RPA-meson loop. It can be shown, that at the external vertex only the scalar interaction 
contributes. Hence ST, is given by 

5E«$o» = -2<7,A(($ », (1.71) 



where A is the constant defined in Eq. ((L29). We should emphasize that these diagrams 



have to be evaluated self consistently at ($0), which is a solution of Eq. ( |1.70| ). Thus 
all "quark propagators", also those which enter into the calculation of the RPA meson 
propagators, have the form S(p) — (J) — ($ )) _1 - Because of the new diagram 8T,, these 
are in general different from the Hartree quark propagators. 

We have to comment on the notations we use here. The utilization of the symbol E 
for the different contributions and the name "quark propagator" for S(p) suggests that 
we again, as in the "zero-loop" approximation, deal here with quark self-energies. But 



x 



o 




O <Q> © <X> 

Figure 1.9: The "one-meson-loop approximation scheme": gap-equation (upper part), and 
the corresponding mesonic polarization functions (lower part). The double line denotes 
an RPA meson propagator (see Fig. which is self consistently constructed from the 

quark propagators dressed with the local self-energy of the present equation (solid line). 
The cross stands for the constant term ($ ) — m o i> n the gap equation (Eq. ( \1.7Q) ). 
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as already pointed out, the interpretation of ($o) as a constituent quark mass should be 
considered with great care. In the Hartree approximation the expectation value of the 
field $o is equal to the pole of the quark propagator. Therefore within this approximation, 
the interpretation of ($o) as a constituent quark mass is correct. The form of the quark 
lines which enter the calculations in the MLA seems to point in the same direction. But 
this is not the case. In Ref. |Kj it has already been emphasized that ($o) is not equal to 
the pole of the "real" quark propagator. This will become clearer if we look at the quark 
condensate. If we take S(p) being a quark propagator literally, we expect to obtain the 
quark condensate by evaluating Eq. ( |1.5|) with the "quark propagator" S(p), 

m = _M<M> = _<M^£ _ A(( * o)) . (1 . 72) 

2g s 2g s 

This result does not agree with the result we obtain following the prescription in Eq. (|1.64j ) 
to determine the quark condensate as 

This expression is reduced to the perturbative result in a strict 1/N C expansion, whereas 
by expanding Eq. (|1.72j ) one only recovers the contribution of 5T,( a \ Moreover, as we will 
discuss in section |1.4.2| , (iptp) according to Eq. ( |1 . 73| ) is consistent with the Gell-Mann- 
Oakes-Renner relation. Hence, within the MLA, Eq. ( |1 . 73|) is the correct expression for 
the quark condensate, whereas the gap equation (Eq. ( |1.70| )) should not be interpreted 
as an equation for the corresponding inverse quark propagator. For simplicity, however, 
we will still call ($o) = m' & "constituent quark mass" and S(p) = (J) — m') -1 a "quark 
propagator", although this is not entirely correct. 

Because of the additional diagram in the gap equation, the RPA is not the consistent 
scheme to describe mesons: In the chiral limit, RPA pions are no longer massless. Con- 
sistent meson propagators can be obtained in the same way from the effective action as 
the RPA meson propagators resulted from the effective action in zero-loop approxima- 
tion. According to Eq. ( |1.50|) , these are given as second-order derivatives of the effective 



action, evaluated at the stationary point. With the help of Eq. ( pL -69| ) this derivative can 
easily be evaluated. This leads to three additional mesonic polarization diagrams which 



are displayed in the lower part of Fig. L~9 together with the RPA polarization loop. Obvi- 
ously these diagrams are identical to SW^ , 8U$ and 811$, which we defined in Sec. |1.3.1| 
(Fig. [D|, Eq. ( |1.33| )), i.e. the new meson propagators are given by 

D,M) = - ' 2 % - - , (1-74) 
1 - 2g M U M (q) 

with 

fl M ( q )=U M (q)+ £ Sn%\q). (1.75) 

k=a,b,c 



Here we have dropped the momentum conserving 5-function (see Eq. ( [1.681 )). As men 



tioned above, all diagrams have to be evaluated at the stationary point, i.e. with a con- 



stituent quark mass m' which is a solution of Eq. ( 1.7C ). For simplicity we renounced to 
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indicate this point by adding an additional mass argument in the above expressions. Since 
the general structure of the expressions we have to evaluate is the same in both schemes, 
we will introduce the convention that unless otherwise stated all quantities within the 
Hartree + RPA and the 1 / iV c -expansion scheme will be calculated with the Hartree quark 
mass rriH and within the MLA with the selfconsistently determined quark mass m' . 

The above discussion clearly reveals that the 1/N C - arguments, put forward by the 



authors of Ref. [29| in order to motivate their choice of diagrams for the extended gap 



equation, Eq. (|1.70|) are questionable. In addition to the fact that in a perturbative 1/N C - 
expansion there would be two corrections to the quark self-energy in next-to-leading order, 
a momentum independent term and a momentum dependent term (cf. Fig. |1.3|) , whereas 
the last one is missing within this non-perturbative scheme, one should stress again, that 
the selfconsistent solution of the gap equation mixes all orders in 1/N C anyway. In addition 
to the systematic approach we presented here following Ref. |5D| , a possibility to engender 
this scheme diagrammatically was discussed in Ref. [@. 

In fact, our discussion of the one-loop approximation to the effective action shows that 
no reference to l/iV c -counting is needed, neither within the diagrammatic derivation in 
Ref. Besides, if one performs a strict 1/N C expansion of the mesonic polarization 



diagrams up to next-to-leading order one exactly recovers the diagrams shown in Fig. |1.4| 
p9| . This is quite obvious for the diagrams 511$ to 5H^, which are explicitly contained 
in Eq. ( |1.75| ). Diagram SH$ , which seems to be missing, is implicitly contained in the 
quark-antiquark loop via the next-to-leading order terms in the quark propagator, which 
arise from the extended gap equation. 

Finally, we would like to stress that this scheme is selfconsistent with respect to the 
quarks, but not with respect to the mesons. This can be seen easily: The quark propaga- 
tor, obtained from the selfconsistent solution of Eq. (|1.70|) , is used in all loops, whereas the 
intermediate meson propagators, i.e. RPA propagators calculated with the selfconsistent 



quark propagators, are not identical to the improved meson propagators in Eq. (|1.74j) . 
We will come back to this point in section 0. 



1.4 Consistency with chiral symmetry 

It was mentioned in the previous section that the effective action formalism provides 
approximations which are consistent with the Goldstone theorem. This can be shown 



on a rather general level |36], pql - The only condition is that the approximation for the 
effective action conserves the symmetries of the action and the Lagrangian respectively, 
which is fulfilled by the loop expansion we performed in the previous section ||36|| . Hence, 
the Hartree approximation + RPA as well as the MLA yield massless Goldstone bosons, 
the pions. Since the mesonic polarization functions of the MLA contain all diagrams up 
to next-to-leading order of the i/iV c -expansion scheme and the various contributions to 
the pion mass have to cancel order by order in the chiral limit, this implies that the 1/N C 



scheme discussed in section |1.3.1| is also consistent with the Goldstone theorem. 



Nevertheless, for the numerical implementation it is instructive to show the consistency 
of the different schemes with chiral symmetry on a less formal level. Since most of the 
integrals which have to be evaluated are divergent and must be regularized one has to 
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ensure that the various symmetry relations are not destroyed by the regularization. To 
this end, it is important to know how these relations emerge in detail. This will also 
enable us to perform further approximations without violating chiral symmetry. As we 
will see in section |2.3| , this is very important for practical calculations within the MLA, 
which cannot be applied as it stands. 

For instructive reasons we begin by an explicit proof of the Goldstone theorem and 
the Gell-Mann-Oakes-Renner (GOR) relation within the Hartree approximation + RPA 



although this is completely standard and can be found in many references (see e.g. [13 



lq , [17], Afterwards this will be discussed within the MLA and the 1/A^-expansion 



scheme. An explicit proof of the Goldstone theorem in both schemes was given first by 



Dmitrasinovic et al. fl29fl . The GOR relation is of particular interest in the context of the 
proper definition of the quark condensate in the MLA. (cf. Eqs. ( |1 . 72| ) and ( |1.73| )). 

To keep the formalism as simple as possible, we restrict ourselves in this section to 
scalar and pseudoscalar interactions. The explicit proof of the Goldstone theorem includ- 
ing vector and axial vector interactions can be found in App. |C.2 . 



1.4.1 Hartree approximation + RPA 

The Goldstone theorem states the existence of massless pions in the chiral limit. Therefore 
one has to show that the inverse pion propagator vanishes in the chiral limit for zero 
momentum, 

2g s U 7T (0) = l for m = 0. (1.76) 
As before we use the notation U^ b = 5^11-^. An evaluation of the RPA loop gives 

2^1^(0) = ^. (1.77) 
m 

We designedly called the constituent quark mass in Eq. ( |1.77| ) m and not m# because 
the validity of this relation between the RPA polarization loop and the Hartree quark 
self-energy is not restricted to the Hartree approximation + RPA itself. If we take this 
relation at the Hartree mass tjih, together with the Hartree gap equation m# = mo + X/f, 
the validity of Eq. (|1.76|) in the chiral limit is ensured. 

Going away from the chiral limit the pion receives a finite mass. To lowest order in 
the current quark mass tjiq it is given by the Gell-Mann-Oakes-Renner (GOR) relation, 

m^/^-moW). (1.78) 

To check the validity of this relation we have to expand m 2 to linear order in m whereas 
f% and the quark condensate have to be evaluated in the chiral limit. The expression 
for the latter is taken from Eq. (|1.6|) . The pion decay constant /„- is calculated from the 
one-pion-to- vacuum axial matrix element. In Hartree approximation + RPA we obtain 

(°) r 

*S ab fl 0) q» = -™ J tJ^I Tr [^r a zS(p + q)z l5 T b zS(p)} . (1.79) 

This has to be taken for on-shell pions, i.e. q 2 = m 2 . The superscript (0) again denotes 
that we deal with a quantity in Hartree approximation + RPA. The above expression 
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resembles the definition of the quark-antiquark loop for the mesonic polarization functions 
in RPA in Eq. ( |1 . 1 1|) . Basically we only replace one vertex by 7Al75 T a /2 and the other by 
i7 5 r b times the pion-quark coupling constant which has been defined in Eq. ( |1.19| ). 

We have to emphasize here that most of the following formal derivations are inde- 
pendent of the actual value of the quark mass. One should keep in mind, however, that 
observables, like /„■, m„ or the pion-quark coupling constant g nqq are only meaningful when 
they are calculated with a quark mass consistent with the scheme used. 

The pion decay constant can most easily be evaluated with the help of the axial Ward 
identity 

7 5^ = 2m 75 +l5 S- 1 (q + p) +S~ l (p)^ . (1.80) 



Contracting Eq. ( 1.79 ) with g M , which means that we take the divergence of the axial 



current, and inserting Eq. ( |1.80| ) on the right hand side of Eq. ( |1.79[ ) we obtain 



✓ Il^g)- 11,(0) 



q 2 



where we have already divided out the common factors i8 a b and have used Eq. (|1 . 77| ) 
to cast the last term in the above form. In the chiral limit, q 2 = — > 0, we can 
use Eq. ( |1.19| ) to replace the difference ratio on the right hand side of Eq. ( |1.81| ) by a 
pion-quark coupling constant. We then arrive at the Goldberger-Treiman relation^, 

fi%& = m H • (1.82) 

Strictly speaking, at this point we encounter an inconsistency if we use another quark mass 
than the Hartree mass m#: We assumed that the pion is massless in the chiral limit, which 
is not the case for the RPA pion calculated with a mass m ^ m#. Therefore, although 
we till now at no stage of the derivation explicitly needed the RPA pion propagator, here 
and in the following part of this section we will use m# instead of a general m. 

For the pion mass we start from Eqs. ( |1 . 1 2| ) and ( |1.19| ) and expand the inverse pion 
propagator around q 2 = 0, i. e. 

1 - 2^(0) - 2 9s (An T ( 9 ))| 52=o mf + 0(mf 4 )=O. (1.83) 

To find m 2 ^ ' in lowest non- vanishing order in m , we have to look at the uiq dependence 
of the different quantities entering Eq. (|1.83|) . Using the relation given in Eq. ( |1.77| ) and 
the Hartree gap equation (Eq. ([TT2"|)) we see that 1 — 2^11^ (0) contains a term linear in 
rriQ. The derivative only implicitly depends on m via the constituent quark mass m# 
and contains no term of linear order. We have to take it therefore in the chiral limit and 
can then replace it by a pion-quark coupling constant, Eq. ( |1.19| ). Thus finally we obtain 
for m 2 }- ' in lowest (linear) order in mo 

m 2 ( o) = ™o (Lg4) 

m H 2g s 

lr The Goldberger-Treiman relation [j56| was originally introduced for nucleons instead of quarks in the 
form f^g-nNN — QAfriN- If we had considered vector and axial vector states, we would have to deal with 
a gA 7^ 1 also in the case of quarks, see e.g. Ref. |16|. 
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Combining Eq. ( |1.6| ) for the quark condensate, Eq. ( |1.82|) , and Eq. ( |1.84j ) we find that the 
Hartree and RPA quantities fulfill the GOR relation (Eq. ( 1.78Q ). 

As already mentioned above, this proof in the Hartree approximation + RPA has 
been listed for purely illustrative reasons. In the following section we will proceed with 
the MLA because in this scheme the proof of the Goldstone theorem as well as the GOR 
relation works in a similar way as in the Hartree + RPA scheme. The main relations are 
formally the same except that we have to replace the Hartree + RPA quantities by their 
MLA analogues. In the l/iV c -expansion scheme, in contrast, especially the GOR relation 
is more difficult to show as we have to carefully expand the various quantities in orders 
of 1/N C . 



1.4.2 One- meson-loop approximation 

Being consistent with the Goldstone theorem requires that, in the chiral limit, the inverse 
pion propagator vanishes at zero momentum, 

2# s fL T (0) = l for m = 0. (1.85) 



The function tl^ b was defined in Eq. ( 1.75 ). It consists of the RPA polarization loop U^ b 
and the three additional diagrams = a, b, c. Restricting the calculation to the 

chiral limit and to zero momentum simplifies the expressions considerably and Eq. (|1.85| ) 
can be proved analytically. 

We should keep in mind that the constituent quark mass is now given by the extended 
gap equation, Eq. (|1.70|) . Therefore, the r.h.s. of Eq. (|1.77|) is different from unity in 
the chiral limit and RPA pions are not massless. This has important consequences for 
the practical calculations within this scheme, which will be discussed in greater detail in 
section |273| . 

The correction terms 8U^ are defined in Eq. ( |1.33|) . Let us begin by diagram 6H^\ 
As mentioned above, we neglect the p and a\ subspace for intermediate mesons. Then 
one can easily see that the external pion can only couple to a ira intermediate state, since 
a oa intermediate state would violate isospin conservation and a intermediate state 
is not possible for parity reasons. Evaluating the trace in Eq. (|1.26| ) for zero external 
momentum one gets for the corresponding 7T7rcr-triangle diagram 

T%^p) = -5 ab AN C N } 2m I(p) , (1.86) 

with a and b being isospin indices and the elementary integral 

, v f d k 1 /x 

^ = J (2tt) 4 (fc 2 - m 2 + ie)((k + p) 2 - m 2 + ie) ' ^ > 



Inserting this into Eq. ( |1.33|) we find 

,ffl(^(0) = i8 ab J -0- 4 (W c N f I(p))Hm 2 DM DM . (U 
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Now the essential step is to realize that the product of the RPA sigma- and pion propa- 
gators can be converted into a difference [p9j1 , 



DM dm 



DM - DM 

AN c Nj 2m 2 7(p) 



to finally obtain 



;i.89) 



5n (a)a 6(0) = ^ 4AyV/ J _g_ 2 /(p){P CT (p) - £>„(?)} . 

The next two diagrams can be evaluated straightforwardly. One finds 

fllW-^O) = -5 ab 4tf c ty / -0u{DM (HP) + HO) - (P 2 - 4™ 2 ) K(p)) 



;i.9o) 



<5IT; c)a6 (0) 



+DM (3/(p) + 3/(0) - 3p 2 tf(p)) } , 
-5 ab 4iV c ty / jAiM{ ~ DM) ~ DM) • (1.91) 



The elementary integral i^(p) is of the same type as the integral I(p) and is defined in 
App. 0. 

Recalling that in the Hartree + RPA scheme Eq. ( 1 . 77| ) played an essential role in the 
proof of the Goldstone theorem we presume that a similar relation is valid for and 511. 
Therefore we look at the correction to the quark self-energy which is given by Eq. ( |1.71| ). 
It is proportional to the constant A, defined in Eq. ( |3.11| ), which is explicitly given by 



A = 4AyV7 m J DM (2 Hp) + 1(0) ~ (p 2 ~ 4m 2 ) K(jp)) 



+DM (3/(0) - 3p 2 K(p))} . 



[1.92) 



Comparing Eqs. ( |1.90| ) and ( |1.91|) with Eq. (|1.92|) we obtain for the sum of the three 
correction terms to the pion polarization function 



E *n«(o) = --. 

, , m 

k=a,b,c 



;i.93) 



Now, Eq. (|1.71 ) can be employed to arrive at an expression similar to the relation be- 
tween the Hartree quark self energy and the RPA pion polarization function for vanishing 
momentum as given in Eq. ( |1.77 ): 



E *n( fe >(o) 



k=a,b,c 



2g s m 



1.94) 



It has already been emphasized in section 1.4.1 that Eq. ( T.77|) is valid for any quark 
mass. Thus we can combine it with Eq. ( |1.94| ) to finally arrive at the following result: 



20 A (o) 



+ 



2g s m 2g s m 
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1.95) 



Hence, together with the modified gap equation (Eq. ( |1.70| )) we obtain 



20.^(0) = 1 - ^ (1.96) 

in agreement with Eq. ( |1.85| ). 

As already pointed out, most of the integrals we have to deal with are divergent and 
therefore have to be regularized. Therefore one has to make sure that all steps which lead 
to Eq. (|1.94j ) remain valid in the regularized model. One important observation is that the 
desired result can be obtained on the level of the p-integrand, i.e. before performing the 
meson-loop integral. This means that the regularization of this loop is not restricted. We 
did not need either to perform the various quark loop integrals explicitly but we had to 
make use of several relations between them. For instance, in order to arrive at Eq. ( 1.90| ) 



we needed the similar structure of the quark triangle r 7r 7r (T (0,p) and the RPA propagators 
D a {j>) and D^p). Therefore all quark loops, i.e. RPA polarizations, triangles and box 
diagrams should consistently be regularized within the same scheme, whereas the meson 
loops can be regularized independently. 

To linear order in the current quark mass mo the pion mass is given by the GOR 
relation (Eq. ( |1.78|) ). This relation holds exactly in the MLA, if we choose the appropriate 
definition of the quark condensate. This will be demonstrated below. 

For the pion decay constant /„- we follow the same steps which led to Eq. (|1.81 



in 



the previous section. It is calculated from the one-pion to vacuum axial vector matrix 
element, which basically corresponds to evaluating the mesonic polarization diagrams, 
Fig. |L4] , coupled to an external axial current and to a pion (cf. Eq. (|1.33|) ), 



9nqq f d*p 



2 J 

, ^9nqq 


(2vr) 4 
r d 4 p 


4 „ 


' (2tt) 4 


^9irqq 


f d 4 p 


2 


1 (2tt) 4 


^9irqq 


f d 4 p 


4 o 


' (2vr) 4 



Tr [ i^^T 0, iS(p + q) 275T 6 iS (p) } 
r ai , Ml ,M2(<l,P) D Mi(p) 

MiM 2 

r 7r ,M 1 ,M 2 (-9) -p) D M2 {~p - q) 



J2 Ta u M lt M u *{q,P, ~P) D Ml {P) , 

N'h 

Mi 



Here the modified pion-quark coupling constant is defined as 

_ 2 GflX(g) 



9n, 



<n 



dq 



2 \q'=m% 



;i.97) 



1.981 



analogously to Eq. ( |1.19|) . Now we take the divergence of the axial current and then use 
the axial Ward-Takahashi identity (Eq. (|1.80| )) to simplify the expressions p9| . One finds 



35 



Sabfn 



9irqq ( 

mH n (q) — mll n (0) 
+m5rf a \q) 1 1 r d ' P 



o J " D n (p+ q ))T a n b ^(- q , -p) 

{D^p)-D a (p + q))T%^-q,- P )} 



+m5H {b \q) + A 

1 r d 4 p 



+ » 



2/ (2^ { ^(p)r?^(p + 



-3D w (p)r^_(p + g,-g)} 



+mOT( c )(g) + I y -^{Z^r^fr + g, -p) 

+ ^(p)r^> + g,-g)}) . 



(1.99) 

To further proceed, let us first consider the terms which still contain an integration over 
p. If one takes into account that the 7T7rcr-vertex should be symmetric with respect to 
exchange of the two pions and performs a shift of the integration variable p such that all 
meson propagators depend on p and not on p + q, one immediately sees that these terms 
cancel. Using Eq. ( |1.93| ) we then arrive at the following expression: 



fn = 9 



m 



U^q) -n^(o) 



;i.ioo) 



In the chiral limit, — > 0, the difference ratio on the r.h.s. of Eq. (|1.100|) can be replaced 
by the pion-quark coupling constant (Eq. ( |1.98| ). This leads to the Goldberger-Treiman 



relation 



fn9 



nqq 



m 



(1.101) 

The same remark of caution as in the Hartree+RPA case concerning the choice of the 
constituent quark mass m should be considered here. Although Eq. ( |1.100|) could formally 
be derived without any reference to the extended gap equation (Eq. ( |1.70| )), for further 
proceeding we should restrict the discussion to the selfconsistently determined quark mass 
ml because otherwise the pion mass would not vanish in the chiral limit. 

For the pion mass we start from Eqs. ( |1.24j) and (|1.25 ) and expand the inverse pion 
propagator around q 2 = 



2g s fL*(0) 



2 9s (^n*(<z) 



1.102) 



As explained in section |1.4.1| , to find m 2 in lowest non- vanishing order in m , we have to 
expand 1 — 2g s U n (0) up to linear order in mo, while the derivative has to be calculated in 
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the chiral limit, where it can be identified with the inverse squared pion-quark coupling 
constant, Eq. (|1.111|) . The result can be written in the form 



<^ s k + ° (<) - (L103) 

Multiplying this by f% as given by Eq. ( |1.101| ) we get to linear order in mo 



m' 



mlf* = mo (1.104) 
2# s 



Obviously this is consistent with the GOR relation (Eq. ( |1.78| )) if the quark condensate is 



given by Eq. ( |1.73| ), but not if it is given by Eq. (|1.72| ). In section |1.3.3| we saw that within 
the effective action formalism the quark condensate is given by Eq. (|1.73p . Therefore at 
this point we clearly see that the interpretation of m! as a constituent quark mass, which 
would mean that we have to calculate the quark condensate according to Eq. (|1.72|) , 



leads to a contradiction with the GOR relation. Therefore, in the numerical part, we will 
calculate the quark condensate according to Eq. ( |1.73|) . 



1.4.3 1/iVc-expansion 



Again we have to show the validity of Eq. ( |1.85| ). In the l/iV c -expansion scheme the 



function 11^ is given by Eq. ( |1.23|) , i.e. it differs from the corresponding function in the 
ML A (Eq. (|1 . 75|) ) by the fact, that it contains the additional diagram <JIlj$ (As we 



already discussed in the MLA it is implicitly contained in the RPA diagram.). The other 
diagrams have the same structure as before and we can largely use the results of the 
previous subsection. 

Since the gap equation is not changed in the perturbative 1/N C expansion, Eq. (|1.76| ) 
remains true, i.e. Eq. ( 1.85|) is already fulfilled by the RPA polarization loop alone. There- 



fore we have to show that the contributions of the correction terms add up to zero, i. e. 

]T 5iTj fc) (0) = for m = 0. (1.105) 



k=a,b,cA 



In the previous subsection we calculated the contribution from diagrams (a) — (c). The 
sum of these correction diagrams is given by Eq. fll.93j ), evaluated at the Hartree quark 
mass. Thus, we are left with the calculation of 511^ (0). According to Eq. fli~33?) , it can 
be written in the form 

<5n?° 6 (0) = -2^(0,0)^(0) A . (1.106) 

The next step is to realize that r* n a (0, 0) can be written with the help of tt- and a-meson 
polarization functions in the RPA, 

ri i<T (o,o) = ^(n.(o)-n CT (o)) . (1.107) 

Till now the evaluation of the various correction terms has been independent of the choice 
of the quark mass. Primarily the next step requires the use of the Hartree quark mass 
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In the chiral limit 11^(0, m#) is equal to l/2g s (see Eq. ( [L.76|) ), and combining Eq. ( 1.1071) 



with the definition of the RPA propagator of the cr-meson in Eq. (|1.12|) one finds that the 



product of the first two factors in Eq. ( |1.106| ) is simply 8 a b/m H , i.e. one gets 

8Tl^ ab (0) = 8 ab — . (1.108) 
m H 



With these results it can easily be checked, together with Eq. (|1.94j) that Eq. (|1.105| ) 
indeed holds in this scheme. 

The discussion concerning the regularization procedure can be repeated here. The 
structure of the proof again leads to the conclusion that we have to regularize all quark 
loops in the same way, whereas we have the freedom to choose the regularization for the 
meson loops independently. 

Another important observation is that we, in both schemes, do not need the explicit 
form of the RPA propagators. D a {p) and D n (p) only need to fulfill Eq. ( |1.89| ). Thus, 
approximations to the RPA propagators can be made as long as Eq. (|1.89|) remains valid. 

In section 1.4. 1| we demonstrated that the Hartree + RPA quantities, i.e. the leading- 



order quantities in 1/N C , exactly fulfill the GOR relation (Eq. (|1.78 )) which determines 



the behavior of the pion mass for small current quark masses mo- This is also true for the 
MLA. In contrary, in the l/iV c -expansion scheme we cannot expect that the GOR relation 
holds exactly. We already discussed in section |1.3.1| that the pion mass, in spite of the fact 
that we perform a strict l/A^-expansion of the mesonic polarization functions, contains 
arbitrary orders in 1/N C because of its implicit definition (see Eq. ( |1.25| )). A second point 
is that we, even if we carefully expand m\ and in orders of 1/N C , generate higher orders 
simply by multiplying both quantities on the left hand side of Eq. (|1.78| ). Thus, because 
we certainly expect the GOR relation to hold in each order in 1/N C separately, we have 
to expand both sides of this relation up to a definite order. In section |1.3.1| we calculated 
the quark condensate in leading order and next-to- leading order in 1/N C . Hence, to be 
consistent, we should also expand the l.h.s. of the GOR relation up to next-to- leading 
order in 1/N C 

ml i0) f 2 J 0) + mT ] 8fl + 8ml f? ] = ~m (<^> (0) + <W>) . (1.109) 

Here, similar to the notations we already introduced for the quark condensate, and 
/^°) denote the leading-order and 8m\ and 8f% the next-to-leading order contributions 
to the squared pion mass and the squared pion decay constant, respectively. In contrast 
to the Hartree + RPA scheme and the MLA Eq. ( |1.109| ) corresponds here to a double 
expansion: Besides the usual expansion up to linear order in m both sides are expanded 
in orders of 1/N C . 

The leading-order and next-to-leading-order expressions for the quark condensate are 
given in Eqs. ( |1.6|) and ( |1.32|) . For the pion decay constant f n we follow the same steps 



as in the previous section. We have emphasized that the derivation of Eq. ( |1.100| ) is 
independent of the choice of the quark mass. Thus this equation also can be employed in 
the present case to describe the contribution of the RPA loop and the sum of correction 
terms (a) to (c). We therefore have to add only the contribution from diagram (d) in 
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Fig. |1.4j, which is given by 



i W£V = IWfa, -g) ,0,(0) A 



Here the l/iV c -corrected pion-quark coupling constant is defined as 



1.110) 



-2 
9ivqq 



,-2(0) + ^-2 = ^Mg) 



Jirqq 



-qq 



dq 



2 I9 2 =m? 



1.111) 



With the help of the axial Ward-Takahashi identity (Eq. ( |1.80| )) and Eq. ( |1.107| ) we find 



/W = ^ ( + m ff ffl (d) (?) - m H mW(Q) + AL> (T (0)(n 7r (g) - 11.(0))) . (1.112) 



g 



Adding up Eqs. ([TTT00|) and ( [1TTT2D we finally arrive at the following expression for the 
pion decay constant in the l/iV c -expansion scheme: 



/ um - r\(o) um - n.(o) 

Jtt = gnqq [ g m# H -5 Ar(0) A 



1.113) 



g^ g^ 

v 0, Eqs. ( p..lS| ) and ( |1.111| ) can be employed to replace the 



In the chiral limit, q 2 = m 
difference ratios on the r.h.s. by pion-quark coupling constants. If we square this result 
and only keep the leading order and the next-to-leading order in 1/N C we finally obtain 



/. 2(0) + *fl = m 2 + (m 2 6g£ + 2mD o (0)Ag^ ) . (1.114) 
Following the analogous steps which led us to Eq. ( |1.103| ) we find for the pion mass 

m gl m f , D a (0)A 



m H 2g s 



(i 



+ 0{mi) . 



1.115) 



Finally one has to expand this equation in powers of 1/N C . This amounts to expanding 
glqq-i which is the only term in Eq. (|1.115 ) which is not of a definite order in 1/N C . One 
gets 



2(0) I c 2 ilirqq 

n * °2g s m 2 



m 



2g s m v m ' 



Combining Eqs. Qg) , ( PH , and ( IlTTBD one finds that the GOR relation in next-to- 
leading order, Eq. (|1.109|) , holds in this scheme. 

However, one should emphasize that this result is obtained by a strict l/iV c -expansion 
of the various properties which enter into the GOR relation and of the GOR relation itself. 
If one takes /„. and m n as they result from Eqs. ( |1.113j ) and ( |1.115| ) and inserts them into 
the l.h.s. of Eq. ( |1.78|) one will in general find deviations from the r.h.s. which are due to 
higher-order terms in 1/N C . As such one can take the violation of the GOR relation as a 



measure for the importance of these higher-order terms [35 . 
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1.5 Transversality of the Rho-Meson 



In this section we want to show that the polarization function in the vector channel fulfills 
the transversality condition, Eq. ( |1.15| ), in the Hartree + RPA scheme, and in the MLA 
as well as in the l/iV c -expansion scheme. Our proof will mainly rely on the Ward identity 

^ = S-\ V + q/2)-S-\ V -q/2) , (1.117) 

where S(p) stands for a quark propagator. We will nowhere need the explicit form of 
this propagator, S 1-1 ^) = i> — ran or S~ 1 (p) = j) — m! respectively. Thus, concerning 
this proof, the only difference between the approximation schemes will be the number of 
diagrams we consider for the entire polarization function. This statement implies that 
the RPA polarization loop and 6TL^ are tranverse themselves, whereas only the sum of 
the three other contributions needs to be transverse. We will explain the procedure in 
detail by proving the transversality of the RPA polarization loop, which was defined in 
Eq. ( |1 . 1 1| ) . For the other terms we will be more brief since the procedure contains no 
essential difference, only the expressions become more lengthy. This is also the reason 
why we will neglect here again the contributions from intermediate vector and axial vector 
states. 

Contracting the RPA loop with and applying the Ward identity, Eq. (|1.117|) , we 
obtain 

q^' ab (q) = -i J ^ Tr [(^(p + |) - S~H P - \))r«iS{p + |) 7 VzS(p - |) ] . 

(1.118) 

This can be further simplified using the invariance of the trace under cyclic permutations, 

q,W^\q) = -i J ^{ T r[zrVr^(p-|)]-Tr[^z5(p + |) 7 V]} . (1.119) 

Performing a shift of the integration variable in the two addends, p' = p— | and p 1 = p+ 1 
respectively, we immediately see that Eq. ( p.,15 ) holds for the RPA polarization loop. As 



already pointed out in section |1.4j in connection with the proof of the Goldstone theorem 
and the GOR relation within the various approximation schemes, we have to keep in mind 
that, most of the integrals have to be regularized since they are divergent. Certainly the 
crucial point is here the shift in the integration variable, whose realization is accompanied 
by a shift in the integration boundary -in some regularization schemes- which would 
destroy transversality. 

Let us proceed with Inspecting the definition in Eq. ( |1.33| ) closely, we conclude 

that we only have to show the transversality of the ppu-triangle vertex. This will later 
on be helpful in connection with diagram (a) in Fig. Oj. Again, contracting the triangle 



vertex for one external p-meson and external mesons M2 and M3, see Eq. ( |1.26| ), with q^ 
and applying the Ward identity we obtain 



,M 2 ,M 3 (q,p) 
d 4 k 

+Tr [r a iS{k - q)T M3 iS(k + p)T M2 ] - Ti[r a iS{k)T Ah iS(k + p)T M2 ]}(l.l2ty 



I (2^{ Tr [ TaiS ( k ) T M 2 iS(k - p)T M3 } - Tr[r a ^(£; + q)T M2 iS(k - p)T M3 ] 
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Performing a shift of the integration variable, k! = k + p in the first two terms, k! = k + q 
in the last two terms, and using the invariance of the trace under cyclic permutations one 
immediately realizes that, if at least one of the remaining mesons, i.e. M 2 or M 3 , is an 
isoscalar state, this triangle vertex is tranverse. This is applicable especially for the case 
of the ppa-vertex. Hence, we can conclude that the correction term 5U^ is tranverse by 
itself. 

Let us now proceed with SH^. As mentioned above, we neglect here the p and a\ 
subspace for intermediate mesons and therefore have to deal only with a two pion inter- 
mediate state. The contraction of the p7T7r-vertex with can be written as a difference 
of RPA pion polarization loops, which in turn can be written as a difference of inverse 
RPA propagators 

q»T^(q,p) = 2ie abc (D-\p)-D-\p + q)) . (1.121) 

This equation represents the Ward identity for the RPA pions. Inserting this result into 
the definition of 5Ii^ we obtain 

% m^ a \q) = -2e acd J t^(At(p + ?) - Ar(p))r#5(-g, - P ) . (1.122) 

Following the same strategy we arrive at a similar result for the sum of the two remaining 
diagrams, 5U^ and 5U^\ 

q,(5Uf^ b (q) + 5U^ b (q)) = -2e acd J ^D n (p)(r^(-q, - p )-r%£(-q,p+q)) . 

(1.123) 

In principle there exists also a contribution to these two diagrams from the cx-meson. A 
straightforward calculation shows that their longitudinal part vanishes. We only have to 
suppose that we can substitute the integration variable p' = —p. If it is allowed to put 
p' = —p — q in the last addend of Eq. ( |1.123| ), we see that the sum of Eqs. ( |1.123| ) and 
( |1.122| ) vanishes, i.e. also the pionic contribution to the p-meson polarization function is 
tranverse. If all the integrals were convergent, this step would certainly be allowed and 
the p-meson polarization function would therefore be transverse. Since this is not the 
case and we have to apply a regularization scheme, we have to pay attention whether 
Eq. (|1 . 1 5|) still holds. In section |2.1| we will come back to this question. 



1.6 Relation to hadronic models 

This section discusses an approximation to our approximation schemes which points out 
the relation to hadronic models. For instance, a possibility to extract a hadronic de- 
scription from the NJL model is to construct an effective Lagrangian from its bosonized 
form, see e.g. Ref. [E7[] and references therein. We will follow here another strategy: In 
order to suppress the quark effects in the present model, it is suggestive to assume that 
the constituent quark mass is very large as compared with the relevant meson momenta. 
Certainly, in the various diagrams we have to evaluate we integrate over meson momenta, 
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i.e. they can in principle be arbitrarily large. Thus we assume, that the main contribu- 
tion to the integrals comes from a region where the momenta are small, i.e. we perform 
an effective low-momentum approximation for the quark-loop vertices. This amounts to 
expanding these vertices up to the first non-vanishing order in the incoming momentum 
("static limit"). In fact, in most cases this corresponds to a zero-momentum approxima- 
tion. In order to preserve chiral symmetry we then have to approximate the RPA-meson 
propagators consistently. It turns out that a straightforward generalization of the pre- 
scription we applied for the vertices to the propagators would lead to a contradiction with 
the Goldstone theorem. To make this obvious let us have a closer look at the way in which 
the latter comes out within this approximation. Since the essential step which determines 
the form of the RPA-meson propagators within this approximation is the same in both 
schemes, the static limit of the ML A and the 1/A^-expansion scheme, we will restrict the 
discussion here to the static limit of the MLA. 

Our aim is now to find a form of the RPA propagators which ensures the validity of 
Eq. Q1.85 ) within the static limit of the MLA, i.e. that the pion polarization function 



vanishes at q = within this approximation, too. For clarity of the expressions we will, 



as in section |1.4j , neglect vector and axial vector intermediate states for the proof. Before 
we look at the pion polarization function we have to investigate the form of the extended 
gap equation Eq. ( |1 . 70| ) in the static limit in order to determine the constituent quark 
mass. The completely momentum independent Hartree self energy remains unchanged 
but the constant A, see Eq. ( |1.29| ), occurring in the correction term to the self-energy ST, 
is approximated by 

A = AN c N f m J D a (p) (3 7(0) + Am 2 K(0)) 

+DM 37(0)} . (1.124) 

For the pion polarization function let us begin by diagram 5U[ a \ The 7T7r<7-vertex is given 
by 

K b ^(0,p = 0) = -5 ab AN c Nf 2m 7(0) , (1.125) 

with a and b again being isospin indices. The proof of the Goldstone theorem in the full 
momentum dependent case relied mainly on the relation between the it- and cr-propagator 
in the RPA, given in Eq. Ql.89|) . By a comparison of Eq. ( |1.86|) and Eq. (|1.89 ), we can 
write this relation in the form 

DMDM = - t mrWM , (1-126) 

where we have used the notation T® b n a = 5 a bT nt7Tta . A straightforward generalization of 
Eq. Ql.126 ) to the static limit would give, explicitly, 



DMDM=i D ^~ D /P\. (1.127) 
nKFJ AN c Nf 2m 2 7(0) V ; 

Exploiting this relation we arrive at the following expression for SU^ in the static limit, 
similar to Eq. ( |1.90 ): 
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m (a )ab{Q) = _ 5ab AKNf J ^P_2I(0){D a (p) - D n (p)} . (1.128) 

Performing the zero-momentum approximation for the vertices in the next two diagrams, 
we find analogously to Eq. ( |1.91| ), 



m (b)a b{Q) = _ 5ab ANcNf J ^{D a (p) (21(0) + Am 2 K(0)) 

+DM 6/(0)} , 

5n W^(0) = -5 ab AN c N f J^I(0){-D a (p)-D n (p)}. (1.129) 

Adding up Eqs. ( |1.129| ) and (|1.128|) and comparing the result with Eq. ( |1.124|) we conclude 
that Eq. (|1.93| ) and thereby also Eq. (|1.85| ) holds within the static limit provided that the 
RPA propagators ascertain the validity of Eq. ( |1.127| ). Thus, in order to be consistent 
with chiral symmetry the RPA propagators have to fulfill Eq. (|1.127|) . Including vector 
and axial vector intermediate states one can show that the propagators of the RPA-p- 
and ai-meson have to fulfill a similar relationship. For all mesons this can be achieved 
by simply replacing the function I(p) in the RPA polarization functions by 1(0). For 
instance, for the inverse a-propagator —2g s D~ l this corresponds to replacing 

l-2s s 2iN c N f (2I 1 -(p 2 -4m 2 )I(p)) by l-2g s 2tN c N f (2I 1 -(p 2 -Am 2 ) 1(0)) . (1.130) 

1(0) can be related to the RPA pion decay constant in this approximation. The 
RPA-meson propagators then take the form of free boson propagators, 

2(0) 

D M (q) = - 9Mqq m , (1.131) 
q 2 - my 

with 

(1.132) 



and 



2(0) = 2(0) = 2 2( o) 2 2 (0) rr^_ 

y-nqq iicjqq ^ ' pqq ^a iq q 2(0) ' ^ J- . ±00 J 

We have to point out that, in a strict sense, this replacement does not correspond to an 
expansion of the RPA polarization functions up to order momentum squared. Consider, 
for instance, the inverse a-propagator. An expansion up to order p 2 would give 

- 2g s D~ 1 = l-2g s 2iN c Nf(2h - p 2 (1(0) + 4m 2 ^M )) (1.134) 

whereas the above prescription leads to 

- 2g s D- x = l-2g 8 2iN c N f (2h - (p 2 - Am 2 ) 1(0)) . (1.135) 
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Figure 1.10: Contributions to the self-energy of the p-meson in standard hadronic models, 
pion loop (left) and pion tadpole (right). The dashed lines denote pions, the curly ones 
p-mesons. 



We proceed by expanding the remaining quark triangles and box diagrams to first 
non-vanishing order in the external momenta, as it was demonstrated above for the 7T7rcr- 
triangle diagram (Eq. ( |1.125p ). In this way one obtains a hadronic model with effective 



meson-meson coupling constants. Let us, for example, look at the p-meson self-energy 
which is generated as an approximation to the NJL model in MLA. Neglecting the p 
and a% subspace for intermediate mesons, we are left with a pion loop diagram, which 
is generated from the diagram shown in Fig. |1.4j (a) and a pion and a sigma tadpole 
diagram coupling directly to the p-meson, which arise from the sum of the diagrams 
in Fig. [0](k) an d (c). In the static limit of the 1 /^-expansion scheme we would, from 
diagram (d) in Fig. [O], in addition get a pion and a sigma tadpole, coupling to the p-meson 
via an intermediate sigma-meson. Assuming the contributions from intermediate sigma- 
mesons to be negligible we exactly recover the diagrams which are calculated in standard 



hadronic descriptions of the p-meson in vacuum, which are displayed in Fig. |1.10| (see e.g. 
Refs. [[5], 0). In fact, a simple reasoning reveals that in a consistent static approximation 
the contributions from the cr-meson have to vanish exactly. As mentioned above, the 
cr-mesons occur in the p-meson polarization function within this approximation only in 
tadpole diagrams, i.e. their contribution is independent of the momentum of the p-meson. 
The pionic contributions, however, preserve the consequences of current conservation on 
their own, in particular H p (q = 0) = for q 2 = 0. Therefore the contributions from the 
u-meson have to vanish for vanishing momentum of the p-meson. Thus, since they are 
momentum independent, they have to vanish identically. 

The comparison of the standard hadronic calculations with the static limit provides 
us with a possibility to estimate the importance of dynamical quark effects. For that 
purpose we will look at the pion loop in the p-meson self energy E M ^. We adopt here the 
usual notation for the self-energy, i.e. the p-propagator is written as 

{q 2 - m p 0)2 - E) _1 T^ + longitudinal parts , 

with X = 1/3T^S M ^ '. Straightforward evaluation of the corresponding diagram gives 

g^{q) = ^{^K) - (q 2 - 4m^(?,m,)} , (1.136) 

with a p7T7r-coupling constant g p7rK . In the above expression we introduced an additional 
mass argument for the integrals Ji and I(q) in order to indicate that they have to be 
evaluated at the pion mass and not, as usual, at the constituent quark mass m. If we 
want to compare this with the static limit of an NJL model calculation we have to pay 
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attention on the fact that the polarization function we usually consider contains a factor 
g 2 pqq as compared with the usual definition of the self-energy, 

1 - 2g M Il M (q) = 1 - 2g M (U M (q) + 5U M (q)) = -2g M g1? q { q \q 2 - m^ 0) + g%° q ' q 5U M (q)) . 

(1.137) 

From the structure of that expression we conclude that £ is given by 

= -9%%SU M (q) . (1.138) 

For the contribution of the pion loop to the p-meson self-energy in the static limit we find 

g,^(q) = lBi^j/xK) - (q 2 - Aml)I(q,m n )} . (1.139) 

Comparing this result with Eq. ( |1.136| ) and replacing the p-quark coupling constant by the 
pion-quark coupling constant using Eq. ( |1.133| ), we obtain for the effective pirn coupling 



constant 

2 

4L = 64°) = 6-^. (1-140) 

If we take the commonly used value of about 6 for g pn7T in hadronic models we obtain for 
the constituent quark mass 

m « V§fi 0) . (1-141) 

We should remark here that there is some uncertainty in the value we have to take for the 
pion decay constant in RPA, f(°\ In the next chapter we will perform a refit such that 
the corrected quantities, i.e. /„. determined either in the MLA or in the l/A^-expansion 
scheme, are equal to the empirical value of f n . As will be seen in section |2.6.1| , is 
larger but not much larger than the corrected quantities. In any case, since we are only 
interested in a rough estimate here we will take a value of ~ 100-150 MeV. According 
to Eq. ( |1.141| ) this leads to a constituent quark mass of about 250-350 MeV. This is 
obviously in contradiction to our original assumption of very heavy quarks . Thus, if 
the main contribution to the momentum dependent quark vertices does not arise from 
momenta much smaller than 250-350 MeV, we would expect that this approximation does 
not describe the full model very well. This point will be discussed in more detail in the 
next section. 

In hadronic models one is also interested in the coupling of a bare p-meson to a photon, 
which is e.g. needed to calculate the pion electromagnetic form factor via vector-meson 
dominance. In the NJL model the bare p-meson corresponds to the RPA meson and its 
coupling to a photon is basically given by the RPA polarization loop in the vector channel, 
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(0) 

/"/</ 



<» = -i^n p (p)T^(p) . (1.142) 

If we now perform the same low-momentum approximations as for the polarization func- 
tions we find that the vertex is given by 

g%(p) = i{e/g p ^)p 2 T^{ P ) (1.143) 
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which exactly corresponds to the 7p vertex in the vector dominance model of Kroll, Lee 
and Zumino ||59| . 

Similar to the proceeding for the p-meson, approximations to the self-energies of the 
other mesons can be performed. For instance, this approximation to the l/iV c -corrected 
NJL model generates an effective one loop approximation to the linear sigma model in 
the 7r-<7 sector. The resulting effective meson-meson coupling constants depend on the 
quark mass and which in turn depend on temperature and density. This dependence 
is used in Ref. within a linear sigma model calculation for in-medium pion properties 
to simulate effects of (partial) chiral symmetry restoration at nonzero temperature and 
density. This undertaking should, of course, be interpreted with great care since all the 
effects which can be modulated by their dependence on the quark mass are closely related 
to unphysical quark effects present in the NJL model due to the lack of confinement. We 
will discuss this question in more detail in Chapters || and [37|. 

Another suggestive approximation to the full momentum dependent calculation is 
to evaluate the vertices for on-shell intermediate mesons instead of performing a low- 
momentum expansion. This approximation is discussed in Ref. [Q . However, at least for 
processes dominated by intermediate pions this gives very similar results to those obtained 
with the low-momentum expansion. 
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Chapter 2 

Numerical results at zero 
temperature 



In this chapter we will discuss our numerical results at zero temperature. We will begin 
by a brief description of the regularization scheme and then discuss peculiarities related 
to the solution of the gap equation in the MLA. After that we will study the influence 
of mesonic fluctuations on quantities in the pion sector, thereby focusing on possible 
instabilities. Finally we will perform a refit of these quantities within the l/A^-expansion 
scheme and the MLA and apply the model to observables in the p-meson sector. 



2.1 Regularization 

Before we begin the explicit calculation, we will have to fix our regularization scheme. As 
discussed in Sec. |1.4| , all quark loops, i.e. the RPA polarization diagrams, the quark trian- 
gles and the quark box diagrams must be regularized in the same way in order to preserve 
chiral symmetry. Besides Lorentz-covariance being a desirable feature for its own, we 
prefer for computational convenience a regularization scheme for the quark loops which 
does not destroy Lorentz-covariance. There are several possible covariant schemes, e.g. 
proper-time regularization, Pauli-Villars regularization, subtracted dispersion relations 
or an 0(4)-cutoff. Some of these have the great disadvantage of displaying the wrong 
analytic behavior of the RPA meson propagators. For instance an investigation of Bro- 
niowski et al. [ BT ] shows that the proper-time regularized RPA meson propagators have 



poles in the complex plane near the imaginary axis. Subtracted dispersion relations lead 
to poles for space-like momenta and an 0(4)-cutoff itself induces a non- analytic struc- 
ture in the corresponding meson propagator. Since, for the calculation of the mesonic 
fluctuations, we need the RPA meson propagators for arbritrary four-momenta, we would 
like to choose a regularization scheme which generates RPA propagators with the correct 
analytic behavior. At first sight, Pauli-Villars-regularization fulfills both requirements, 
the covariance and the analyticity of the RPA meson propagators. We will therefore use 
this regularization scheme. Soon we will become aware that Pauli-Villars regularization is 
not completely free of peculiarities, either. We will, however, be able to circumvent most 
of the occurring problems by slightly modifying the standard Pauli-Villars scheme. 
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Figure 2.1: Inverse o -propagator in RPA, D~ l (p 2 ), multiplied by —2g s as a function of 
momentum squared, entire polarization loop regularized (solid line) and only I\ and I(p) 
regularized (dashed line). The parameters used are listed in Table p?J| , for a meson cutoff 
A M = MeV. 



To be consistent, all divergent terms should be regularized with the same number of 
regulators, which is determined by the highest occurring degree of divergence, i.e. the 
quadratic divergence of the integral 

r d 4 k 1 
1 ~ J (2vr) 4 k 2 - m 2 + ie 

which we encounter, for instance, in the Hartree gap equation (Eq. (|1.4j )), 

mil = mo + 2ig s 4N c ]S jmifii(m#) . (2.2) 
A quadratic divergence requires two regulators, thus we replace 

Sm^t^-j^tcmiH), (2.3) 



j=0 



by 

1$ = m 2 +jA 2 q ; c = 1, c x = -2, c 2 = 1 . (2.4) 

Here A q is a cutoff parameter. Let us illustrate the standard way to apply this prescrip- 
tion within the NJL model with the RPA polarization loop in the u-channel, defined in 
Eq. ( p.. lip . A straightforward evaluation gives 

f afk / 2 p 2 - Am 2 \ 

^ P ' = 1 c f J (2^U 2 -m 2 + ie ~ (k 2 -m 2 + ie)((k+p) 2 -m 2 + ie)) ' 



The most obvious procedure would be to replace the entire integrand in Eq. (|2.5|) by the 
sum over the regulators, see Eq. (|2.3p . This would be more in the spirit of the original 
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Pauli-Villars regularization |62j] than the usual strategy [1?|, where the polarization loop 



is written in terms of the already defined elementary integrals I\ and I(p), 

IL a (p) = 2iN c N f (2h - (p 2 - Am 2 )I(p)) . (2.6) 

Here only the integrals I\ and I(p) are regularized, the factor m 2 in front of I(p) is not 
modified. The advantage of this strategy can be seen most easily from the inverse a- 



propagator which is displayed in Fig. |2.1| for both methods. This calculation has been 
performed with the Hartree quark mass. The naive application of Eq. ( |2.3| ) to the entire 
polarization loop obviously leads to a zero of the inverse propagator, i.e. a pole in the prop- 
agator, for space-like momenta. The same observation can be made for the ai-propagator. 
This is the reason why we will adopt the usual form of Pauli-Villars regularization in the 
NJL model and regularize only the divergent elementary integrals throughout this inves- 
tigation. 

We should add a comment on a rather general problem which is related to the regu- 
larization of the RPA polarization loops. A straight-forward evaluation of the vector and 
axial vector polarization diagrams gives 

n p (p) = -i^N c N f (-2I 1 + (p 2 + 2m 2 )I(p)), (2.7) 
n ai (p) = -i^N c N f (-21, + (p 2 - 4m 2 ) I{p)) . (2.8) 

Because of vector current conservation U p should vanish for p 2 = 0. This is only true if 

m 2 I(0) = h, (2.9) 

which is not the case if we regularize I(p) and I\ as described above. Alternatively one 
could perform the replacement Eq. ( |2.3| ) for the entire polarization loop. Then the factor 
m 2 in Eqs. Q2.7Q and (|2.9j ) should be replaced by a factor /i| inside the sum over regulators 
and one can easily show that Eq. ( |2.9| ) holds (see Eqs. ( |A.10| ) and ( |A.12| )). So it seems 



as if the strategy we rejected above since it generates poles in the RPA a- and ai-meson 
propagators for space-like momenta would be favorable to apply here. However, this 
scheme leads to another severe problem: From the gap equation (Eq. |2.2| ) we conclude 
that il\ should be positive. On the other hand the pion decay constant in the chiral limit 
and in leading order in 1/N C , i.e. RPA, is given by, see Eq. ( |1.81|) 



/2(o) = _ 2iNcNf m 2 1(0) , (2.10) 

which implies that im 2 I(0) should be negative. Thus, irrespective of the regularization 
scheme Eq. (|2.9| ) cannot be fulfilled if we want to get reasonable results for m# and 
at the same time. Therefore we will retain our choice of the standard form of Pauli-Villars 
regularization in the NJL-model and replace the term Ii in Eq. ( |2.7|) by hand by m 2 1(0). 
This leads to the following p-meson propagator 

Dp ^ P) = 1 + 2ig v ±N c N f (-2m 2 7(0) + (p 2 + 2m 2 ) I(p)) ' (2 ' H) 
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A priori it is not clear whether the symmetry properties of our approximation schemes 
can be preserved if we perform the above subtraction for the vector polarization function. 
We noted, however, that this is the case if we take care that the relationship between the 
RPA p- and ax-meson propagators, Eq. ( |C.11|) , is not affected. This can only be achieved 
if the ai-meson polarization function is treated analogously to the vector one. We then 
arrive at the following RPA ai-meson propagator 



1 + 2ig v lN c N f 



'-2m 2 1(0) + (p 2 - 4m 2 ) 7(p)) 



(2.12) 



After having explained the regularization scheme applied to the quark loops, we will 
come to discuss the remaining meson- loop integration (integration over d 4 p in Eq. ( p..33| )). 
This integration is not automatically rendered finite by regularizing the quark loops due 
to the non-renormalizability of the model, so that we have to regularize it separately. This 
regularization is not constrained by chiral symmetry and independent from the quark-loop 
regularization. For merely practical reasons we choose a three-dimensional cutoff Am in 
momentum space. In order to obtain a well-defined result we work in the rest frame of the 
"improved" meson. The same regularization scheme was already used in Refs. |39], |4(], [12j . 

From an esthetic point of view, it would certainly be desirable to employ a regulariza- 
tion scheme which allows to render both quark loops and meson loops finite without the 
introduction of an additional parameter (or as in our case, even another regularization 
scheme). This could be achieved in a NJL-type model with a separable non-local inter- 
This non-local interaction generates a form factor at the quark vertices. 



action 
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That this makes all integrals finite can most directly be seen in the version of Ref. 



The authors of Ref. |28] use a form factor which depends on the three-momentum of the 
quarks. It can be chosen in such a way, that the three-momentum of each quark is limited 
to absolute values below sum cutoff parameter A. Then simple kinematical considerations 
lead to the conclusion that the absolute value of the three-momentum of the mesons is 
automatically restricted to lie below 2A. The disadvantage of this model is that it is 
manifestly non-covariant. In Refs. |63| similar models with form factors depending on 

These are obviously covariant, but they have the 



In Refs. [ 

the four-momentum are considered, 
disadvantage of destroying the analytic properties of the RPA meson propagators. Re- 
cently mesonic fluctuations have been investigated [§4, [yj within the model presented in 
Ref. p§. 



2.2 Evaluation of the various two- loop diagrams 

Within this section we will present some details concerning the rather involved numerical 
evaluation of the two-loop diagrams [] contributing for instance to the mesonic polarization 
functions or the modified gap equation in the MLA, Eq. ( |1.70| ). 



A general two-loop diagram contains a four-dimensional (momentum) space integra- 
tion for each loop, i.e. altogether eight integrations. Some of these can usually be per- 



^ote that diagram (a) in Fig. 1.4 is in principle a three-loop diagram. But, since the two quark loop 
vertices decouple, we actually have to deal only with a two-loop diagram. A similar remark concerns the 
RPA mesons which in principle contain an arbitrary number of loops. 
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formed analytically while the remaining ones have to be performed numerically. Of course 
one is interested in performing as few as possible integrations numerically since comput- 
ing time rises and computing precision is reduced steeply the more integrations one has. 
In our case all the quark loop integrals, i.e. those for the RPA meson propagators and 
the effective meson-meson vertices, can be written in terms of five "elementary" integrals 
defined in App. [D] which in turn, with one exception, can be evaluated analytically. Only 
for the function 

u , _ f d A k 1 

\Pl,P2,P3) J (2?r)4 ^ 2 _ m 2 + ie ^},2 _ m 2 + te ^ k 2 _ m 2 + ^2 _ m 2 + k ) ' 

(2.13) 

with ki = k + pi, an analytic expression exists exclusively in kinematical regions where 
at least one of the momenta lies above the corresponding (in our case quark-antiquark) 
threshold, i.e. p\,p\,p\, (P1+P2) 2 , (P1+P3) 2 or (P2+P3) 2 is larger than 4m 2 . This function 
appears in the four-meson vertex T Ml)M2tM , i ^ Mi (p ll p2 1 pz) ) Eq. Ql.271 ). Thus, we are at 
most left with a one- dimensional integration for the quark loops. 

Since we work in the rest frame of the mesons the angular part of the meson-loop 
integration merely contributes a factor of 4tt and therefore only an integration over energy 
and the absolute value of the three-momentum remains. The former is most difficult to 
handle because we have to take care of the various singularities and discontinuities arising 
from poles and thresholds of the RPA meson propagators and the quark loop vertices. 
Let us consider for instance the pionic part of the energy integrand for SU^ in the scalar 
channel 



3z 

8^ 



/•Am roc 

/ d\p\ dpop^r ann (q,p-q/2)D^(j) + q/2)r atnn (-q,-p + q/2)D w (j)-q/2) 

JO J -co 

3 Mm roc 

= g^3 I d \Pl /_ oo d Po f(Po, \p\,q) , (2.14) 

where we have included a factor 3 due to the sum over the isospin indices of the intermedi- 
ate pions. We have performed a shift of the integration variable p — > p + q/2 as compared 
with the definition of 5U^ in Eq. ( |1.33| ) in order to obtain a function f(po, \p\,Q) which 



is symmetric in p , enabling us to restrict the energy integration to the interval [0, oof by 
replacing f(po, \p\, q) with 2f(p , \p\, q). This symmetry property becomes obvious if we 
take into account that the RPA pion propagator D v (p) is only a function of p 2 instead 
of p and that the ainr- vertex function T a ^ j7T (q,p) only depends on the momentum of the 
incoming mesons squared, i.e. q 2 , (p — q/2) 2 and (p + q/2) 2 . The real and imaginary parts 
of the symmetrized integrand are shown in Fig. [2.2| for \p\ = 100 MeV and two values of qo, 
100 MeV and 500 MeV. We have used the parameters listed in Table JH\ for A M = MeV . 
The most pronounced peaks we observe are the pion poles at p = \ ± q /2 + J m 2 + p 2 \ . 
Note that we have introduced by hand a small but nonzero imaginary part by replacing 
the quark mass m in the quark loop integrals with m — ie, e ~ 1 MeV. One consequence 
of this nonzero imaginary part is the finite width of the pion poles. More or less steep 
quark-antiquark thresholds appear if one of the momenta, i.e. q 2 , (p + q/2) 2 or (p — q/2) 2 , 
equals 4m 2 . These are invisible in Fig. |2^ because of the much larger scale of the singu- 



larities arising from the pion poles. Additional thresholds exist due to the Pauli-Villars 
regulator masses m 2 + A 2 and m 2 + 2A 2 , see Eq. (j2.4|). 
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Figure 2.2: Integrand for the po-integration, pionic part of the contribution 6TL^ to the a- 
meson polarization function with \p\ = 100 MeV, real (solid) and imaginary part (dashed), 
for q = 100 MeV (left) and q = 500 MeV (right). 



The values of go = 100 MeV and go — 500 MeV have been chosen for the following rea- 
son: The former value for g lies below any threshold and we expect no imaginary part to 
be present at the end, the latter above the two-pion threshold at 280 MeV. The results we 
obtain for the imaginary part below any thresholds are very sensitive to the details of the 
numerical evaluation. Calculating this imaginary part, which we know to vanish, serves as 
a sensitive measure for the precision of our delicate numerical integration procedure. An 
imaginary part is always related to a decay process, e.g. a — > nn. The three-dimensional 
cutoff we apply for the |p| -integration can at most restrict the phase space for that decay. 

If for example go > l^jm^ + A^, then no decay is possible, thus the cutoff can suppress 
a decay process but in no way generate one, and consequently no imaginary part. From 
this we conclude that the energy integration alone has to be responsible for the imagi- 
nary part vanishing below any threshold. Thus the integration over the imaginary part 
of f(po, \p\, q) has to vanish for g < 280 MeV, e.g. for g = 100 MeV. Looking at Fig. ^]2| 
one realizes that this is a completely non-trivial undertaking. Particular attention should 
be paid to three crucial points: to proper integration over singularities and thresholds, to 
an extrapolation to vanishing e, and to the inclusion of high-energy contributions. The 
first and the second point are connected, since the quantity e has been introduced by 
hand to soften the otherwise numerically not tractable singularities and thresholds. This 
artificial width enables us, by suitably adapting the location and number of mesh points 
in the vicinity of singularities and thresholds, to perform the integration with sufficient 
precision. The limit e — > is approximated thereafter by linear extrapolation. Including 
also contributions from very high energies turns out to be difficult because the evaluation 
of the integrand becomes numerically unstable for very large energies, i.e. po > 15 GeV. 
However, at least for the imaginary part these contributions cannot be neglected since 
they are necessary for cancellations below any threshold. An integration up to po — > oo 
can be achieved by first applying one of the standard substitutions of the integration 
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variable: For large energies, i.e. above any threshold, we substitute po = tan(x). The 
resulting function of x can very well be approximated by a quadratic function and thus 
can easily be extrapolated up to the upper integration boundary of tt/2. The real part 
remains almost unaltered by the operations described above. 




500 



Figure 2.3: Integrand for the \p\-integration, pionic part of the contribution SH^ to the 
a -meson polarization function, real (solid) and imaginary part (dashed), with (left) go 
I ()() MeV and (right) g = 500 MeV. 

The finally rather satisfying result can be seen in Fig. |2.3| where we have plotted the 
corresponding integrand for the ^-integration. The left panel represents the result for 
go = 100 MeV, the right panel that for g = 500 MeV. For g = 100 MeV the imaginary 
part is suppressed by 5 orders of magnitude as compared with the real part which reveals 
the efficiency of our method of integration. The |p|-integration is comparatively easy to 
handle. For go = 100 MeV the integrand is completely smooth, for go = 500 MeV only 
in the vicinity of \p\ = \J ql/A — m^, i.e. \p\ = 207 MeV, we have to be a little careful 
with placing the mesh points for the integration. In most other cases, more precisely in 
all cases where we do not have to deal with a two-pion intermediate state, we also end up 
with a completely smooth function for the |p|-integration. 



2.3 Solution of the gap equation in the ML A 



In the l/iV c -expansion scheme all diagrams are calculated with the constituent quark mass 
in Hartree approximation, i.e. as it emerges from the Hartree gap equation, Eq. (|1.2|) . In 
the MLA, on the contrary, the quark mass is determined by the solution of the extended 
gap equation, Eq. ( |1.70|) , which was derived in section |I . 3 . 3| from a one- meson- loop approx- 
imation to the effective action. The structure of this extended gap equation was indispens- 
able in order to prove the validity of various symmetry relations, namely the Goldstone 



theorem, the Goldberger Treiman relation and the GOR relation, in section |1.4.2| . This is 
not astonishing since the entire scheme is based on the extended gap equation. However, 
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only the structure of this equation has been used but not the explicit solution. This sec- 
tion is now devoted to a discussion of the explicit solution of the extended gap equation. 
Unforeseen problems occur which will force us to slightly modify the scheme. 

First indications for these problems can already be inferred from the behavior of the 
RPA mesons as a function of a trial constituent mass m. As already pointed out, the 
RPA mesons, which enter the extended gap equation (see Fig. |1.9|) , have to be evaluated 
with the selfconsistently determined quark mass m', which is in general different from 
the Hartree quark mass m#. Hence, the propagators of these mesons differ from those 
in the Hartree + RPA scheme. To illustrate this let us look at the masses of the pion 
and u-meson in the RPA as a function of a trial quark mass m displayed on the l.h.s. of 
Fig. 2.4 . We have used the parameters listed in Table 2T with Am = MeV. The effects 
of 7r-ai-mixing can be governed by manipulating the vector coupling constant g v : With a 
vanishing g v , 7r-ai-mixing is turned off, with an increasing g v it becomes more and more 
important. For comparison we have computed the pion mass with g v — (solid line) and 
g v = 2g s (dotted line). The mass of the a-meson (dashed line) does not depend on g v . 
The principal observation is that the pion becomes tachyonic, i.e. m^ 2 becomes negative, 
for quark masses smaller than the Hartree quark mass. Strictly speaking this is only true 
in the chiral limit, for nonvanishing current quark masses we find negative m^ 2 at first 
for quark masses slightly below m#. 

Qualitatively the results do not change if 7r-ai-mixing is included. If we keep all 
parameters fixed, it can be seen from the explicit form of the pion propagator in the RPA 
(Eqs. QC.4D and ( B.5 )), that in the chiral limit the results even do not change at all. Away 
from the chiral limit the pion mass is slightly increased and hence the trial quark mass, 
where the pion mass vanishes, is slightly decreased. In principle, the parameters should 
be refitted with the pion mass fixed to its empirical value at m = m#, which diminishes 
the influence of 7r-ai-mixing on the value of the quark mass m where the pion becomes 
tachyonic. 

A similar observation can be made for m^ 2 , but only for m much smaller than the 
Hartree mass. This finding of tachyonic RPA mesons for certain quark masses causes 
the meson- loop term in the effective action (second term in Eq. ( |1.53|) ) to be no longer 
positive definite in that range of quark masses. As discussed in section |1.3.3| this reveals 
that the effective action is ill-defined in that range. 

The latter point becomes obvious if we look at the quark self-energy which becomes 
complex if it is calculated with tachyonic RPA mesons. Therefore we have to hope that 
a solution of the extended gap equation can be found which is larger than the Hartree 
mass. Any other solution necessarily must be complex or negative. To illustrate this point 
we consider the difference between the l.h.s. and the r.h.s. of the extended gap equation, 
Eq. (|1 . 70| ) , as function of a (real) trial mass m. This difference is plotted in the right panel 
of Fig. ^4| . Let us begin the discussion by the results obtained with g v = 0, i.e. we turn 
off effects of 7r-ai-mixing and the contributions of vector and axial intermediate states in 
<5£. The solid line denotes the real part for this case, the dashed line the imaginary part 
of m — m — E(m). Evidently the self-energy indeed becomes complex for quark masses 
below some mass m c . A comparison with the pion mass reveals that m c coincides with the 
quark mass below which the pion becomes tachyonic. A kink in the real as well as in the 
imaginary part indicates the point where also the sigma becomes tachyonic. Moreover, we 
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see that there is no solution of the gap equation for real constituent quark masses besides 
the trivial solution at m — 0. Hence, solutions, if there exists any, lie somewhere in the 
complex plane off the positive real axis. Consequently the RPA mesons would be built of 
quarks with either complex or negative masses leading to peculiar properties. In any case 
this renders a reasonable description of the mesons in MLA completely impossible. This is 
most clearly seen for the p-meson, whose properties are mainly determined by intermediate 
RPA pions. Including vector and axial vector intermediate states as well as 7r-ai-mixing 
does not alter this troubling result, since, qualitatively, the function m — m — E(m) shows 
the same behavior in the case of g v = 2g s . It is displayed on the r.h.s. of Fig. |2.4j , the 
dashed-dotted line corresponds to the real part, the dotted line to the imaginary part. 

A possible way out of this dilemma is to perform a modification, which has been 
introduced in Ref. |3(|. Of course, any approximation we perform has to preserve the 
symmetry properties of the MLA. As discussed in section |1.4.2| this can be achieved if 
the validity of Eq. ( |1.89|) for the RPA meson propagators is not affected. This remains 
true if we include vector and axial vector intermediate states, see App. |U.2| . In that case 
in addition to Eq. (|1.89| ) for the pion and a-meson propagators a similar relationship for 
the RPA p- and ai-meson propagators has to hold. Thus the approximation suggested in 
Ref. 0] for the purely scalar and pseudoscalar case can straightforwardly be generalized 
to the case with vector and axial vector interaction. This will be shown in more detail 
in App. |(J.2| , here we will only illustrate the principle of this approximation which can be 
done most easily if we set g v = 0. 

The authors of Ref. [30] simply replace the RPA pion propagator in the extended gap 
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Figure 2.4: (Left) Squared pole masses of the pion with g v = (solid), g v = 2g s (dotted) 
and the a -meson (dashed) in RPA as functions of a trial constituent quark mass m in 
units of the Hartree quark mass. (Right) Difference m — m — between the l.h.s. and 
the r.h.s. of the MLA gap equation, Eq. ( \1.70j ), as a function of the trial constituent quark 
mass m. The real part is denoted by the solid line for g v = and by the dashed-dotted 
line for g v = 2g s , the imaginary part by the dashed line for g v = and by the dotted line 
for g v = 2g s . 
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equation, 



D v (p) = -2g s 



by 



2ig s AN c N f 



DM = -2g s 



(2tt) 4 k 2 - m 2 + is 



2tg s (2N c N f )p 2 I(p) (2.15) 



-i 



m 



m 



+ 2ig s 2N c N f p 2 I{p) 



(2.16) 



and analogously for the a-propagator. It is clear from the proof of the Goldstone theorem 
in section 11.4.21 that the same replacements have to be performed for the RPA meson prop- 

c- (k) 

agators in the correction terms 5H M to the mesonic polarization diagrams. In addition 
chiral symmetry requires that the RPA contribution to the mesonic polarization function, 
Hm, itself is not changed. By performing this approximation we not only obtain a real 
(positive) solution of the extended gap equation but also the masses of the intermediate 
mesons remain real. Moreover, we obtain massless intermediate pions in the chiral limit, 
which becomes obvious by looking at Eq. ( |2.16|) . Besides, the various symmetry relations 
proved in section |1.4.2| still hold. 

We can certainly ask whether there is any other motivation for this approximation in 
addition to the practical reasons discussed above. In the Hartree approximation, i.e. if 
the RPA mesons are calculated with quarks of mass m#, the above replacements would 
be exact. Following the arguments put forward by the authors of Ref. |HJ, the correction 
terms we neglect by performing the above replacements would be suppressed, because 
they are of higher order in 1/N C . This is of course a questionable argument in the MLA 
since already the quark mass mf, the selfconsistent solution of the extended gap equation, 
contains arbitrary orders in 1/N C . In fact, we have seen that this approximation drastically 
changes the results. 

We will keep the name "one-meson-loop approximation" for this scheme, including the 
above replacements, although this is strictly speaking only an approximation to the MLA 
as it was originally introduced in section |1.3.3| . 



2.4 Meson- loop effects on quantities in the pion sec- 
tor 

The additional cutoff Am we introduced in section |2TT1 , enables us to study the influence of 
mesonic fluctuations on several quantities, e.g. the quark condensate or the pion mass m n 
and the pion decay constant f n . This cutoff restricts the momentum of the intermediate 
meson states. We can therefore, by slowly turning on this cutoff, control the importance 
of the mesonic fluctuations. For Am = we of course recover the usual RPA result. 
It is advantageous to keep all other parameters constant. We use mo = 6.13 MeV, 
A q = 800MeV and g s A 2 = 2.90. These parameter values are chosen in such a way that 
mf£\ and (ipip)^ reproduce the empirical values. We obtain, for g v = 0, = 140.0 
MeV, ffl = 93.6 MeV and ($V) (0) = -2(241.1 MeV) 3 . For the constituent quark mass in 
the Hartree approximation we get m# = 260 MeV. These values are slightly changed if we 
include vector mesons, i.e. g v ^ 0. In the following section we will discuss our results for 
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the quark condensate, in section 2.4.3| we will consider the properties of the pion, among 



other things m w and f n which are related to the quark condensate via the GOR relation. 

Of course, an application of the approximation schemes, the l/iV c -expansion scheme 
and the MLA, to physical processes will only be meaningful if a refit of the parameters is 
performed to reproduce the empirical values of m^,^ and (ipip)- This will be discussed 
in section [2.6. 



2.4.1 Quark condensate 

The behavior of the quark condensate as a function of the mesonic cutoff Am is displayed 
in Fig. |2.5|. The upper panel shows the results in the chiral limit and the lower panel 



those with mo = 6.13 MeV, i.e. with a realistic pion mass. The solid lines correspond to 
g v = 0, i.e. excluding vector interactions, and the dashed lines to g v = 2g s . We will begin 
the discussion by the results obtained without vector and axial vector intermediate states 



which were already investigated in Ref. [Q. We observe that in the chiral limit as well 



as with a nonvanishing current quark mass the absolute value of the quark condensate 
decreases first with an increasing cutoff and then, if the cutoff exceeds a certain value, 
goes up again. This feature is common to the l/iV c -expansion scheme (left panel) and the 
MLA (right panel). This effect is related to an unphysical pole in the pion propagator 
we detect for large values of the cutoff, see section P~5] . Studying the effects of vector 
and axial vector fluctuations is an interesting issue in the context of a suggestion made 
by Plant and Birse [fUl , saying that these instabilities do not survive if one takes into 



account fluctuations of vector and axial vector mesons and especially the 7r-ai-mixing. 
They corroborate this suggestion by a calculation of the quark self-energy in a generalized, 
non-local, version of the NJL model (cf. section |2?1| ). Without vector meson intermediate 
states they encounter in next-to- leading order in 1/N C a minimum in the self-energy for 
large momenta, whereas the vector and axial vector states seem to suppress this effect. 

Considering the results with g v = 2g s we should primarily state that the influence 
of the vector and axial vector fluctuations is remarkably big. In the 1 /A^-expansion 
scheme the minimum does not vanish, neither in the chiral limit nor with a nonzero 
current quark mass. It is only relocated at higher values of the cutoff. It becomes even 
more steep compared with the case without vector and axial mesons. In the MLA the 
minimum has entirely disappeared. Looking at the results in the chiral limit we instead 
observe a behavior typically for a phase transition: For small values of the cutoff the 
quark condensate rapidly approaches zero, at Am ^,450 MeV it vanishes. The slope of 
the condensate at the "critical cutoff" diverges, indicating indeed a phase transition of 



second or first order. In section |0| we will establish a conclusive answer to the nature of 
this phase transition from the effective potential. In that section we will also explain how 
the calculations in the chiral limit, numerically a difficult task, are actually performed. 
For mo = 6.13 MeV the condensate shows the typical cross-over behavior: For small 
values of Am it rapidly approaches zero, at Am ~ 550 MeV the slope rather abruptly 
changes such that the quark condensate seems to reach zero asymptotically. 

The authors of Ref. |41| state that mainly the longitudinal parts of the vector mesons 



are responsible for the disappearing of the minimum. In our case, as shown in section [L5 



the ppu-vertex is transverse and hence the longitudinal part of the p-meson gives no 
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contribution at all, whereas the longitudinal part of the axial polarization, absorbed in 
the pion propagator via 7r-ai-mixing, can very well be important. 

To investigate this question we will look at the constant A which is the main element 
of the corrections to the quark condensate, both in the l/A c -expansion scheme and in 
the MLA. In the 1 /^-expansion scheme this can be seen directly (cf. Eqs. ( |1.6| ) and 
( |1.32| )): As a function of Am the Hartree contribution and the RPA cx-propagator remain 
constant. Thus the mesonic fluctuations can only have an effect on A. In the MLA this 
connection is not as obvious since we have to solve Eq. ( |1.70| ) selfconsistently. Nevertheless 
the qualitative behavior of the quark condensate can be deduced from that of A. In order 
to understand the peculiar "turn-around" as a function of the meson cutoff Am, it is 
convenient to look at the integrand instead of A itself. The function displayed 




Figure 2.5: The ratio (ipip) / (ip^) ^ as a function of the meson cutoff Am- The solid line 
corresponds to g v = 0, i.e. the result without vector and axial intermediate states, the 
dashed line to g v = 2g s , (left) 1/N C - expansion scheme, (right) MLA. The curves in the 
upper panel are calculated with a vanishing current quark mass, in the lower panel the 
corresponding results with a realistic value of m = 6.13 MeV are shown. 
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in Fig. 2J3, has to be integrated over \p\ to obtain A. This representation allows us to 
estimate the influence of the cutoff value on the result for A. We show the contributions 
of the different mesonic intermediate states and the corresponding sum, on the left hand 
side we neglected the p- and ai-meson subspace, whereas on the right hand side it is 
included (together with 7r-ai-mixing). 

Let us begin our discussion by the purely scalar-pseudoscalar case. The sum of both 
contributions first rises, then, at |p| ~ 445 MeV, it starts to decrease. At \p\ ~ 900 MeV 
it even crosses the axis and its slope indicates that it will become more and more negative 
for larger values of \p\ . We see that the pionic contribution is much more important than 
the contribution from cr-meson intermediate states. It also mainly produces the negative 
contributions which are responsible for the "turn-around" of the quark condensate as a 
function of the cutoff. The contribution of the cr-meson of course remains completely 
unchanged if we include vector and axial intermediate states. The pionic one, however, 
is affected by vr-ax-mixing. As can be seen on the r.h.s. of Fig. [2.6| the decrease of the 
pionic contribution is weakened and shifted to higher momenta. The (transverse) vector 
contribution is mainly responsible for the decrease of the integrand, whereas the tranverse 
part of the axial one, i.e. the ai, tends to increase it. The maximum and the zero of 
the overall sum is relocated at higher absolute values of the three-momentum compared 
with the purely scalar-pseudoscalar case, but qualitatively the behavior is not altered. 
This observation is directly reflected in the behavior of the quark condensate in the 
l/Af c -expansion scheme, where vector interactions only lead to a shift of the observed 
minimum as a function of the mesonic cutoff Am to higher energies. In the MLA, however, 
the relation between A and the quark condensate is less direct mainly because of the 
selfconsistent determination of quark condensate. In any case also in the MLA the general 
behavior remains unaltered. 




Figure 2.6: Integrands for A as a function of \p\. (Left:) g v = 0, contribution from the 
pion (dashed line), a-meson (dotted), and the sum of both contributions (solid). (Right:) 
g v = 2g s , contribution from the pion (dashed), a-meson (dotted), p-meson (short- dashed), 
a\-meson (dashed- dotted) , and the sum of all contributions (solid). 
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2.4.2 Problems with vector intermediate states 



Before we begin our discussion of pionic properties let us explain the difficulties we en- 
counter if we try to include vector and axial vector meson intermediate states for the 
calculation of the corrections to the mesonic polarization functions. The problem is re- 
lated to the rather general problem discussed in section |2.1| of preserving the consequences 
of current conservation for the vector propagator in RPA, i.e. transversality and n p (0) = 0, 
and simultaneously not inducing wrong analytical properties for a- and ai-mesons. As an 
example we will consider the diagram (b) in Fig. |1.4|. It contains among other contribu- 



tions one four-meson vertex where two external pions with momentum q couple to an axial 
state with momentum p. Contracting this vertex with p^Pv/p 2 we receive contributions 
to the integral for SH^(q), defined in Eq. (|1.33|) , from an ai-meson as well as from the 
axial part of the pion. In fact, the resulting scalar vertex function, 



■PuPt 
— i- 



P 2 



-12JV C JV>( lip + q)(p + qf - I(q)(q 2 + p • q) - l(0)p • q 

+ K(q)q 2 (p -q-p 2 )+ Am 2 p 2 K(p) - Am 2 p 2 q 2 L(p, -q, 0))(2.17) 



is multiplied by the difference of the RPA ai-meson propagator and the axial part of the 
pion propagator, D ai {p) - D Vv (p). 

The problem arises from the combination of two types of terms: The first type are 
addends proportional to l/p 2 , e.g. the term q 2 /p 2 I(q) in Eq. ( |2.17[ ), the second type 
are terms depending only on the momentum q, i.e. being constant with respect to the 



integration variable p, e.g. the term q 2 K(q) in Eq. ( 2.17 ). Thus, on the one hand the 
difference of the RPA ai-propagator and the axial part of the pion propagator should 
vanish as p 2 for p 2 — > in order not to produce unphysical poles at p 2 = 0, and on the 
other hand it should vanish at least as 1/po for large energies to guarantee the convergence 
of the energy integration. Using the explicit expressions for the ai-propagator and the 



axial part of the pion propagator given in Eqs. ( 2.12|) and ( |U.4j ), respectively, we see 



that the difference vanishes at p 2 = 0. Hence we can conclude that the first condition 
is fulfilled for the RPA propagators in their present form. Since the three-dimensional 
cutoff Am restricts the absolute value of the three-momentum \p\, we can gain insight into 
the behavior of the difference for large energies by looking at it for large p 2 . Using the 
properties of the Pauli-Villars regulators introduced in Eq. ( |2.3| ) we infer that the integral 
I(p) in leading order behaves like 

I(p)-» 2^5 for p 2 >m 2 . (2.18) 

This enables us to conclude that the difference of the RPA ai-propagator and the axial 
part of the pion propagator contains besides terms vanishing at least with l/p 2 a constant 
term proportional to 

-m 2 I(0) + h . 

This term is reminiscent of the a\ polarization function. Note that the subtraction we 
performed in section |2.1| in order to preserve gauge invariance for the vector polarization 
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Figure 2.7: Squared pion mass as a function of the current quark mass m for differ- 
ent meson-loop cutoffs; (left) 1/N C - expansion scheme: Am = MeV (solid), 500 MeV 
(dashed), 900 MeV (dashed- dotted) and 1300 MeV (dotted); (right) ML A: A M = MeV 
(solid), 300 MeV (dashed), 500 MeV (dashed- dotted) and 700 MeV (dotted). The calcu- 
lated points are explicitly marked. 



function here induces this problem. On the other hand a closer examination shows that 
the difference would not vanish at p 2 = if we had not performed this subtraction. Thus 
the only solution seems to be providing for Eq. ( [2.9|) to hold. This is a dilemma since we 



saw in section |2J] that this leads to even more severe problems. We therefore decided to 
put the vector and axial vector meson intermediate states aside and will only treat scalar 
and pseudoscalar states throughout the remaining part of this paper. 



2.4.3 Pion Properties 

In this subsection we want to study the influence of mesonic fluctuations on the pion 
propagator, which will enable us to study especially the effects on the pion mass and 
the pion decay constant, both within the l/A c -expansion scheme and within the MLA. 
Major parts of this analysis can be found in Refs. |39|, We proceed in the same way 



as in the previous subsection: Since the strength of the fluctuations is controlled by the 
meson cutoff A^, we first keep all other parameters fixed and investigate how the above 
quantities change, when A M is varied. 

Fortunately, the pion mass provides us with a decisive test for the stability of the rather 



involved numerics. The Goldstone theorem, which was proved analytically in section [L4 
for the i/iV c -expansion scheme and the MLA, clearly states that the pion mass has to 
vanish for a vanishing current quark mass rriQ. For numerical reasons we cannot compute 
the pion mass exactly in the chiral limit, but its value can be extrapolated from the 
behavior of the squared pion mass as a function of the current quark mass m . This 
is displayed in Fig. ^L7\ for different values of Am, on the l.h.s. evaluated within the 
l/A c -expansion scheme, on the r.h.s. in the MLA. Obviously the dependence in both 
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Figure 2.8: The ratios m^/m^ ^ (solid), f%/f% (dashed), and the combination 
—mo(ipi(})/m^f^ (dotted) as a function of the meson loop cutoff Am in the 1/N C - expansion 
scheme (left) and the ML A (right). 



schemes is almost linear, i.e. the analytically calculated point at m = 0, = and the 
numerically computed points, corresponding to the same cutoff, lie almost on a straight 
line. This result is suited for convincing us of the consistency of our calculations with 
chiral symmetry and of the stability of the numerics. 

Now we fix the current quark mass again to mo = 6.13 MeV and slowly turn on the 
mesonic cutoff Aj^ and with it the strength of the mesonic fluctuations. The resulting 
behavior of and f% as a function of Am is displayed in Fig. |2~lj| . The left panel shows the 
results for the l/iV c -expansion scheme, the right panel for the MLA. It is clearly visible, 
that the qualitative behavior of f w and m n turns out to be the same in both schemes: 
The first (dashed lines) is reduced whereas the latter (solid) is increased. If only a and 
7r intermediate states are considered it was observed in the previous section that both 
schemes show similar results for the quark condensate. 

In section |1.4.1| we showed that in the Hartree approximation + RPA the quantities 
m^°\ fl^ and (■?/>?/>) < -°\ should be in good agreement with the GOR relation, Eq. Q1.78 



This should also hold for the corresponding quantities m n ,f n , and {ipip) in the MLA, 
see section pL . 4 . 2| . The l/A c -expansion scheme, on the contrary, is consistent with the 
GOR relation only up to next-to-leading order and, as discussed in section |1.4.3| , the 
"corrected" quantities violate the relation by higher-order terms. For this reason we 
expect a less perfect agreement in this scheme, becoming worse with increasing values of 
Am- 

These expectations are more or less confirmed by the results. In the Hartree + RPA 
scheme we find only deviations -due to corrections in higher orders in m - of about 0.1%. 
In the MLA the deviations are slightly larger, due to the numerics being much more 
difficult to handle. One reaches up to about 3% for a cutoff of Am = 900 MeV. These 
deviations can be inferred from the ratio of the r.h.s. and the l.h.s. of the GOR relation, 
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Eq. (|1.78|) , displayed in Fig. ^8] by the dotted lines. As can be seen, in the MLA (right 
panel) it almost coincides with the desired value of 1, whereas in the l/iV c -expansion 
scheme (left panel) we indeed observe considerable discrepancies. Within this scheme 
the GOR relation holds at least within 30% for Am < 900 MeV, tolerable for higher- 
order corrections in this perturbative scheme. However, when the meson cutoff is further 
increased the agreement with the GOR relation rapidly deteriorates, indicating that in 
this regime higher-order corrections become important. Our perturbative scheme should 
not be trusted there any more. 

One may wonder why we do not plot the various curves for cutoffs larger than Am = 
1250 MeV in the l/iV c -expansion scheme and Am = 900 MeV in the MLA. The reason 
is that upon further increasing Am we detect a second, unphysical, pole in the pion 
propagator |39|, This phenomenon will be discussed in more detail in the next section. 



2.5 Instabilities due to strong mesonic fluctuations? 



As already mentioned it has recently been claimed by Kleinert and van den Bossche P5| 
that in the NJL model chiral symmetry, which is spontaneously broken in the mean field 
ground state, gets restored due to strong mesonic fluctuations. One principal objection 
one can raise against this supposition is the non-renormalizability of the NJL model. 
Including fluctuations beyond mean-field we cannot avoid encountering additional diver- 
gencies. Thus a priori an additional cutoff parameter has to be introduced as discussed 
previously. This additional cutoff parameter then controls the strength of mesonic fluc- 
tuations [Pm . Nevertheless it might be possible to observe some kind of "chiral symmetry 
restoration" at a certain value of the cutoff parameter. 

In the previous section we stated that a second, unphysical, pole emerges in the pion 
propagator if the cutoff exceeds a certain value. This pole has a residue of the "wrong" 
sign, leading to an imaginary pion-quark coupling constant and a negative value for the 
squared pion decay constant f 2 . This becomes clear from Fig. |2.9| , where the inverse pion 
propagator D^ 1 is plotted as a function of the momentum squared q 2 for different values 
of the cutoff Am- The left panel shows the results in the l/iV c -expansion scheme, the 
right panel those in MLA. In the chiral limit all curves would go through —2g s D~ l = 
at q 2 = with a negative slope. In the Hartree approximation + RPA the effect of a 
non-zero current quark mass is to shift this value up to mo/m^, i.e. the zero of the inverse 
propagator, corresponding to the pole of the propagator, is then located at some positive 
(time-like) value of q 2 . In the 1 /^-expansion scheme the main effect of the meson-loop 
contributions for Am < 900 MeV is to further increase the value of D^ 1 at q 2 = and 
with it the value of q 2 at the zero, i.e. the pion mass. Upon further increasing Am the 
value of D~ l (0) is again reduced. From Eqs. ( |1.93|) and ( |1.108|) we conclude that the 



l/A^-correction terms to D n 1 (0) are proportional to A. Thus, the behavior of D 7r 1 (0) 



reflects that of the constant A, discussed in section |2.4.1| . At the same time the curve 
flattens considerably, leading to a strongly rising pion mass. For Am ?t 1250 MeV we 
observe that it turns around below the quark- ant iquark threshold, causing the second 
pole in the pion propagator. Obviously the slope of D^ 1 at this second zero is inverse 
as compared with the first zero. This reveals the above mentioned fact of a residue with 
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Figure 2.9: Inverse pion propagator, D^ 1 , multiplied by —2g s as a function of the 
^-momentum squared for various meson cutoffs; in the 1/N C - expansion scheme (left): 
A M = (solid), 900 MeV (dashed), 1300 MeV (dashed- dotted) and 1500 MeV (dotted) 
and the ML A (right): A M = MeV (solid), 500 MeV (dashed), 700 MeV (dashed- dotted) 
and 1100 MeV (dotted). 



the "wrong" sign. When the cutoff is further increased the two zeros move towards each 
other, merging at Am ~ 1350 MeV and disappearing thereafter from the real positive 
g 2 -axis. Qualitatively the same behavior is observed within the MLA. The only difference 
is that the second unphysical zero emerges already for A M > 910 MeV and the two poles 
merge already for Am ~ 1095 MeV. 

Within both schemes this unphysical behavior can be traced back mainly to the con- 
tributions from 5Tl£\ This term has a peak in the imaginary part above the gg-threshold 
with the "wrong" sign. This can be seen from Fig. |2.10| , where the various contributions 
to the pion polarization function are displayed. For low values of the cutoff the sum 
of all other contributions cancels this effect such that the imaginary part of the entire 
polarization function behaves reasonably. In the l/A^-expansion scheme especially ^II^, 
proportional to A, adds to this cancellation. Here the growing of the peak in STlffl to- 
gether with the diminishing value of A for Am > 900 MeV lead to the imaginary part 
of the polarization function changing sign if the cutoff exceeds a certain value. This 
imaginary part with the "wrong" sign can then be related via dispersion relations to the 
"turn-around" in the real part below the gg-threshold. 

The appearance of the unphysical second pole already indicates that the pion becomes 
unstable if the strength of the mesonic fluctuations is increased too much. In Ref. 
we suggested that these instabilities of the pion propagator might indicate an instability 
of the underlying ground state. This was motivated by referring to similar situations in 
many-body systems, where e.g. complex energy eigenvalues in the excitation spectrum 
are a hint at an instability of the ground state ||65|| . In that reference we exclusively 
dealt with the l/A^-expansion scheme, where it is evidently impossible to investigate the 
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Figure 2.10: Real (upper panel) and imaginary part (middle) of the different 1/N C - 
correction terms to the pion polarization function, 6U.fi (long-dashed), 6Ufi (dotted), 
6Ufi (short-dashed), and 5Ufi (dashed- dotted), for A M = 700 MeV (left), and 1100 MeV 
(right). Different corrections term to the pion polarization function in the MLA for 
A M = 1000 MeV (lower panel), 6U^ (long- dashed), 6U^ (dotted), and 6U^ (short- 
dashed), real part (left) and imaginary part (right). 
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question whether the ground state develops an instability against mesonic fluctuations, 
since these are built perturbatively on the Hartree ground state. However, the MLA 
offers a possibility to examine the influence of mesonic fluctuations on the structure of 
the ground state more closely f42fl . Since we encounter the same type of instabilities within 



this scheme as in the l/iV c -expansion scheme, we can probably gain deeper insight into 
the nature of these instabilities. The principal tool to study the structure of the ground 
state is certainly the effective action, Eq. ( |1.53| ). In our case it can be written in terms 



of the effective potential, defined in Eq. ( |1.51| ), which describes the energy density of the 
system. It is explicitly given by 



rr/ x _ «AT at f a P i , m P \ (m — m ) 
V(m) = 2iN c N f / — Vrlnf^ — £-) + — 

-\ J 4^u{Hl - 2g s U a (p)) + 31n(l - 2g s Il 7r (p))} + const. . (2.19) 



The irrelevant constant can be chosen in such a way that V(0) = 0. Here we only 
considered 7T- and cr-meson as intermediate states. Of course this can be generalized to 
include also p- and ax-meson states. The resulting explicit expression can be found in 
App. |03. 



Since we plan to study a possible "chiral symmetry restoration" the only sensible 
procedure would be to perform the calculations in the chiral limit. We have already 
mentioned in connection with the pion mass that it is numerically an extremely difficult 
task to evaluate the expressions in vacuum exactly in the chiral limit. Thus we are forced 
to build this limit either by extrapolating from nonzero current quark masses or from 
nonzero temperatures. We have checked that both methods lead to almost the same 
results. To be precise we should point out that we will proceed in the following way 
to build the limit: Starting from the parameter set given above, we keep the Hartree 
constituent quark mass, win = 260 MeV, fixed, while mo is reduced from 6.1 MeV to 
zero. The consequence is a slightly enhanced coupling constant, g s A 2 = 2.96 as compared 
with g s A 2 q = 2.90. The same procedure was also applied to the calculations of meson loop 
effects on the quark condensate described in section |2.4.1 . 



In which way can we now identify a possible "chiral symmetry restoration" in the 
behavior of the effective potential? First let us remind ourselves that the positions of the 
extrema of V{m) correspond to the solutions of the gap equation ( |1.70|) . As the system 
always tends to minimize its energy, the vacuum expectation value m! is therefore given 
by the value of m at the absolute minimum of V. For a vanishing current quark mass m' 
is, according to Eq. ( |1.73|) , proportional to the quark condensate and consequently to the 



order parameter of chiral symmetry. Hence, for a given value of Am, a non-zero value of 
m! indicates a spontaneously broken ground state whereas chiral symmetry is restored if 
the absolute minimum of V is located at m = 0. 

V(m) as a function of mjmu for different values of Am is displayed in Fig. [2.11[ 
Taking into account only fluctuations of n- and a-mesons we obtain the results shown 
on the l.h.s. in the other case those on the r.h.s. For Am = we of course encounter 
the usual picture: A local maximum is located at m = and we find a minimum at 
m = niu — 260 MeV. This is the well established result in the Hartree approximation. 
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Figure 2.11: Effective potential V(m) as a function of m/mn in the MLA for different 
values of the meson cutoff Km, on the l.h.s. without vector meson intermediate states, 
i.e. g v = 0, on the r.h.s. with g v = 2g s . (left) MeV (solid), 300 MeV (long-dashed), 
500 MeV (dotted), 900 MeV (dashed- dotted) and 1200 MeV (short- dashed), (right) MeV 
(solid), 300 MeV (dashed), 400 MeV (dotted) and 500 MeV (dashed- dotted). 



We have already explained that we would expect the maximum at m = to convert into 
a minimum upon increasing Am if there was indeed some kind of "phase transition" . It is 
obvious from Fig. 2.11 that in the purely scalar and pseudoscalar case no phase transition 
occurs, although the results for Am < 900 MeV seem to point in this direction. In this 
regime not only the vacuum expectation value ml is considerably reduced -to about 30% 
of the Hartree mass for Am = 900 MeV- but also the "bag constant", which is given as 
the difference in energy between the maximum at m = and the minimum at m = rnf, 
i.e. B = V(0) — V(jn'), moves towards zero. It decreases from 48.7 MeV/fm 3 at Am = 
to 0.8 MeV/fm 3 at Am = 900 MeV. However, when Am is further increased we also 
observe the traces of the "turn-around" of A (see section |2.4.1|) : Both, ml and B, go up 
again. This means in particular that V(0) always remains a local maximum leading to 
the conclusion that a "phase transition" is ruled out as an explanation for the instabilities 
we found. 

Note that the results do not completely exclude the possibility of a "phase transition" 
due to strong mesonic fluctuations. The small bag constant at A M = 900 MeV indicates 
that probably a careful parameter study brings to light a "phase transition" in a certain 
parameter range. This supposition is corroborated by the results including the p- and 
ai-meson subspace for intermediate states. As we already suspected from the results for 
the quark condensate in section |2.4.1| , we observe a second order phase transition. The 
existence of a phase transition is clearly discernible from the behavior of the effective 
potential for different cutoffs displayed on the r.h.s. of Fig. |2.11| . For Am = MeV 
the curve certainly does not differ from that in the purely scalar and pseudoscalar case. 
However, as already stated in the previous sections, the vector and axial intermediate 
states amplify the effects of mesonic fluctuations, leading to an accelerated decrease of 
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ml and B with an increasing cutoff. For Am = 500 MeV the maximum has disappeared 
and the only extremum we find is a minimum at m — 0. A closer examination of the 
development strongly suggests that we deal here with a second order phase transition with 
a "critical" A^ iC = 461.82 MeV. The exact value of this "critical" meson cutoff should 
not be taken too seriously since we have chosen the value of g v = 2g s rather deliberately. 
For instance, if we take g v = g s the critical cutoff is augmented by almost 100 MeV. 

We have to remark here that our investigations are not suitable to give a conclusive 
answer as to whether for a certain choice of parameters "chiral symmetry restoration" 
according to the conjecture by Kleinert and van den Bossche [I3| is observed. The phe- 
nomenon the authors of Ref. f|3|| describe is closely related to one from strong-coupling 
superconductors. In addition to the superconducting phase in strong-coupling supercon- 
ductors there exists a so-called pseudo-gap phase above T c |66|, [67|, [68| , where Cooper pairs 
are still formed but do not condense. Analogously in such a phase the quarks would still 
acquire a non-vanishing constituent mass if the mass was identified with the modulus of 
the field $ introduced in section |1.3.3| . Since the phase of the <3>-field is strongly fluctuating 
the expectation value of (ipip) would nevertheless vanish, indicating that chiral symmetry 
would not be broken within this phase. Our approach implicitly assumes a uniform phase 
factor and we can therefore not exclude a phase transition into a pseudo-gap phase of 
that type. 

In any case we would like to understand the reason for the instabilities we found. 
A hint is probably the discussion of another type of vacuum instability which has been 
performed recently by Ripka |64| . In section |2J] we pointed out that several regularization 
schemes produce unphysical poles of the RPA propagators, i.e. poles in the complex plane 
where they are forbidden by microcausality. Ripka states that the instabilities he finds 
by examinating the effective potential are caused by these unphysical poles of the RPA 
propagators, which are induced by the regularization procedure. His analysis is performed 
with a 4-momentum cutoff and a Gaussian form factor. As discussed in section [2.1| at 
the first sight the Pauli-Villars scheme does not generate such difficulties. However, the 
regulators lead to an "overshooting" of the imaginary part of the RPA propagators, i.e. 
in some kinematical regions a gg-continuum is created with the wrong sign. Although 
the above discussion, in connection with the observation that the instabilities in the pion 
propagator can be traced back mainly to an imaginary part with the wrong sign in 5H%\ 
strongly corroborates the supposition that the instabilities are related to the analytical 
properties of the RPA propagators. Certainly further investigations are necessary in order 
to give a conclusive answer. From that point of view it is also not very astonishing that 
recently a second (unphysical) pole in the pion propagator has also been found in a non- 
local generalization of the NJL model jy] using a Gaussian form factor similar to that 
discussed by Ripka |64j] . The above reflections render it implausible that these instabilities 
could principally be removed by including vector and axial vector intermediate states, as 
suggested by Plant and Birse JO], since the analytical properties of the RPA propagators 
are thereby not altered. Besides, we observed in section [2.4.1| that the behavior of the 
quark condensate, strongly related to the instabilities we found in the pion propagator, 
is not changed qualitatively if we include vector interactions. Probably vector and axial 
vector intermediate states could enlarge the region in parameter space, where we are far 
off the region where instabilities occur and at the same time obtain a reasonable fit of 
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observables. We will enter into this question in the next section. 



2.6 Parameter fit 

Up to now we have investigated the influence of mesonic fluctuations on various quantities 
in the pion sector. To that end all parameters have been kept fixed, which have been 
determined by fitting f@\m£\ and ($ijj)( Q \ except the meson cutoff Am- It is clear that, 
if the model should be applicable to describe physical processes, a refit of the parameters 
to reproduce the empirical values of /„., m n , and {4>ip) should be performed. For the 1/N C - 



expansion scheme this has been done in Ref. [40]. In the MLA, however, we did not achieve 



such a fit [B2J for reasons explained below. Of course, the three observables in the pion 
sector do not suffice to conclusively determine the five model parameters m , g s , g v , A q , 
and Am- The fact that we are compelled to exclude vector and axial vector intermediate 
states, see section |2.4.3| , considerably simplifies the fitting procedure because in this case 



pionic properties do not depend on the vector coupling constant g v . We will therefore 
follow the strategy established in Refs. [fH], For various values of Am we first fix the 
current quark mass m , the quark- loop cutoff A g , and the scalar coupling constant to fit the 
quantities in the pion sector, i.e. m n , f n , and {ipip). With the two remaining parameters, 
g v and Am, we will then try to reproduce the data for the pion electromagnetic form factor 
in the time-like region which is, assuming vector-meson dominance, primarily affected by 
the properties of the p-meson. The latter is very well suited for this purpose since it, 
besides being a vector state, cannot be described reasonably without including mesonic 
fluctuations. Since the phenomenologically important two pion intermediate state in our 
model consists of RPA pions, we are forced to fit m^' and not m n to the empirical pion 
mass if we desire to get the correct threshold behavior for the p-meson. This is certainly 
a slight inconsistency and will restrict the applicability of our model to regions where the 
deviations are not too large. 

Another point additionally restricts the regime of possible parameter sets: The unphys- 
ical gg-threshold must lie well above the peak in the p-meson spectral function in order to 
obtain a realistic description of the p-meson. Consequently the constituent quark mass, 
tuh in the l/A c -expansion scheme and m! in the MLA, has to be larger than about 400 
MeV. We have to emphasize here that the constituent quark mass is entirely determined 
by the particular parameter set fixed by fitting observables. Nevertheless we have some 
freedom: The empirical value of the quark condensate is not known very precisely because 
it is not a directly measurable quantity. We have to rely on sum rule analyses or results 
from the lattice. From the former an upper limit for the absolute value of (ipip) of about 
2(260 MeV) 3 at a renormalization scale of 1 GeV can be extracted |6!|. The value obtained 
on the lattice lies in the same range, recent results give (ipip) = -2((231±4±8±6) MeV) 3 
[70]. Thus, as long as the value of the quark condensate stays within the above-given 



boundaries, we can try to increase the constituent quark mass as much as possible. In 
the MLA it will be much more difficult to obtain a quark mass which is large enough 
than in the l/A c -expansion scheme. The reason is obvious: The meson-loop corrections, 
which in the MLA affect the constituent quark mass via the solution of the extended gap 
equation, Eq. ( 1.70Q , tend to diminish it. In the l/A c -expansion scheme, on the contrary, 



69 



the determining quantity is the Hartree quark mass tuh, which is not influenced by any 
mesonic fluctuations. In addition, in the l/iV c -expansion scheme, the mesonic fluctuations 
lower the absolute value of the quark condensate, whereas in the MLA the direct relation 
between the quark condensate and the constituent quark mass (Eq. ( |1.73|) ) at the same 
time leads to a decreasing quark mass. 



2.6.1 Pion sector 



Our results jH], Q in the pion sector are listed in Table [2TT| for the l/A^-expansion scheme 
and in Table |2.2j for the MLA. We show five parameter sets, corresponding to meson-loop 
cutoffs of MeV (RPA), 300 MeV, 500 MeV, 600 MeV, and 700 MeV, together with the 
respective constituent quark mass and the resulting values of m n , f T , and (ipip). Since the 
properties of the intermediate pion states, in particular the pion mass , are important 
in connection with the p-meson, we also display the RPA quantities m^, and (^ip)^ . 
In the MLA, however, these are strictly speaking no RPA quantities since within that 
scheme the intermediate "RPA" states consist of quarks with the constituent quark mass 
m! . This, among other consequences, results in the quark condensate (ipip) coinciding 
with the "RPA" value cf. Eqs. (|L~6l ) and (|1.73|) . We therefore renounced to list 



Am / MeV 


0. 


300. 


500. 


600. 


700. 


A, / MeV 


800. 


800. 


800. 


820. 


852. 


m / MeV 


6.13 


6.40 


6.77 


6.70 


6.54 




2.90 


3.07 


3.49 


3.70 


4.16 


m H / MeV 


260. 


304. 


396. 


446. 


550. 


/ MeV 


140.0 


140.0 


140.0 


140.0 


140.0 


m T 1 MeV 


140.0 


143.8 


149.6 


153.2 


158.1 


/ MeV 


93.6 


100.6 


111.1 


117.0 


126.0 


U 1 MeV 


93.6 


93.1 


93.0 


93.1 


93.4 


/ MeV3 


-2(241. I) 3 


-2(249.3) 3 


-2(261. 2) 3 


-2(271.3) 3 


-2(287.2) 3 


(V>i) / MeV 3 


-2(241. I) 3 


-2(241. 7) 3 


-2(244.1) 3 


-2(249. 5) 3 


-2(261.4) 3 




1.001 


1.007 


1.018 


1.023 


1.072 



Table 2.1: The model parameters (Am, A q , m and g s ) and the resulting values ofm n , f n 
and {^tp) (together with the corresponding leading-order quantities) , the constituent quark 
mass mn in the 1/N C - expansion scheme. The ratio -mo(^)/m^/^, is also given. 



both values in Table 2.2. For comparison we also give the corresponding Hartree quark 



mass in Table |2.2| . The ratio of the right hand side of the GOR relation, Eq. (|1.78|) , 



divided by the left hand side, cf. the last line in Tables EO and|2.2j, provides us on the one 



hand with a measure for higher order (beyond next-to-leading order) corrections in the 
l/Af c -expansion scheme, and on the other hand with the possibility to estimate numerical 
uncertainties. The latter can be deduced from the discrepancies between the exact value 
of unity and the obtained value in the MLA. The almost perfect agreement of the RPA 
quantities with the GOR relation suggest that we deal here indeed mainly with numerical 
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/ MeV 


140.0 


145.1 


156.3 


164.5 
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/ MeV 


93.6 


96.7 


103.6 


108.4 


120.0 


U 1 MeV 


93.6 


93.2 


92.9 


92.9 


92.8 


(^) / MeV 3 


-2(241.1) 3 


-2(241.7) 3 


-2(246.2)3 


- 2(250.8)3 


-2(260.9)3 


-m (ip^)/mlf% 


1.001 


1.001 


1.006 


1.01 


1.02 



Table 2.2: The same as in Table \2 
with {ipip) , see text. 



for the ML A. (iptp)( ' is not shown since it agrees 



uncertainties and not with higher-order corrections in mo- We obtain deviations of at 
most 2%, assuring our confidence in the numerics. With the above parameter sets we 
also stay in a region where higher order corrections in the l/iV c -expansion scheme remain 
small: the deviations in the l/iV c -expansion scheme are less than 10%, for Am < 600 MeV 
even less than 3%. 

Although we find in both schemes that the constituent quark mass, m# in the 1/N C - 
expansion scheme and m' in the ML A, increases with an increasing cutoff, the values in the 
ML A are considerably smaller than those in the l/A^-expansion scheme. Within the latter 
we find a region of meson cutoffs, 500 MeV < Am < 700 MeV, where on the one hand 
the constituent quark mass is large enough to shift the gg-threshold above the peak in the 
p-meson spectral function, and where on the other hand we can reproduce the empirical 
values of fn,m n and (ipip). The upper boundary for Am is caused by the impossibility to 
simultaneously stay below the limit for (ipip) and reproduce the correct value for /„- upon 
further increasing Am- If one remembers the above discussion concerning the constituent 
quark mass, it is not astonishing that the region of meson cutoffs where we obtain sensible 
results is much more narrow in the MLA. From the values listed in Table [2.2| it can be 
seen that for Am = 600 MeV m! is still too small whereas for = 700 MeV the value 
of (ipip) already slightly exceeds the limit. The main reason is that simultaneously with 
the constituent quark mass the absolute value of the quark condensate rises, rendering it 
much more difficult to stay below the limit for (ipip). 



2.6.2 Description of the p-meson 

In the previous section we started the parameter fit by determining parameter sets for 
various values of the meson cutoff A M , which reproduce the empirical values of m^, f n , 
and {^i/j). In both schemes we established a lower and an upper boundary for the meson 
cutoff Am- The lower boundary was chosen such that the quark threshold lies above the 
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p-meson mass, whereas the upper boundary mainly resulted from the impossibility to 
reproduce the empirical values of {^i/j) and f n at the same time. In this section we will 
attempt to decisively fix Am by looking at quantities related to the p-meson. This will 
allow us as well to fix the vector coupling constant g v . In the l/iV c -expansion scheme 



we will succeed [[IIJ, whereas in the MLA we will enter into a region where instabilities 
occur [42], similar to those observed in the pion propagator (see section |2.5| ). This will 
become more clear from the following discussion. 

Including meson-loop effects the polarization function of the p-meson reads (cf. 



Eqs. (pi) and (|L75D ) 

tl^ a \q) = Il^ ab (q) + £ fllf) ^ a \q) . (2.20) 

k 

Here k runs over {a, b, c, d} in the l/iV c -expansion scheme and over {a, b, c} in the MLA. 
It is understood that the functions have to be evaluated at the Hartree quark mass ran 
in the l/iV c -expansion scheme and at m — ml in the MLA. In section |1.5| we proved, 
with the help of Ward identities, that this polarization function fulfills the transversality 
condition, Eq. ( |1.15[ ), provided we apply a regularization scheme which preserves this 
property. The delicate point was the shift in the integration variable needed to establish 



for instance the cancellations in Eq. (|1.119|), which is accompanied by a shift in the 



integration boundary in some regularization schemes, destroying transversality. This is 
obviously not the case for the Pauli-Villars scheme we employed to regularize the RPA 
polarization loop. With the three-dimensional cutoff Am we use to regularize the meson 
loops we have to be more careful: In principle not only Lorentz covariance is violated but 
also the transverse structure of the correction terms 511 to the polarization function is 
in general not retained. However, we are working in the rest frame of the p-meson, i.e. 
q = 0, so at least transversality can be preserved. Thus, if we assume Lorentz invariance 
to be fulfilled, the ansatz for the tensor structure of the polarization function in Eq. (|1.16| ) 
is also applicable, enabling us to evaluate a scalar function tl p (q) = —\/?>g^ v Ii^ u {q) (cf. 
Eq. ( |1.17| )) instead of all tensor components. 

A second consequence of vector current conservation is, that the polarization function 
should vanish for q 2 = 0. In section |2J] we saw that already in the Hartree approximation 
+ RPA this can only be achieved by a subtraction. The reason should not be searched 
for in a particular regularization procedure we have applied, in this case Pauli-Villars, 



but to a greater extent in the rather general problem we noted in section |2J] concerning 
the inconsistency of several regularized RPA quantities. The three-dimensional cutoff we 
apply to render the meson-loop integrals finite violates current conservation directly, i.e. 
the correction terms to the p-meson polarization function do not vanish a priori. We cure 
this problem again by performing a subtraction, 

E 8np\ q ) — E K%) - *nf( )) • (2.21) 

k k 

We are now in a position to attack the evaluation of the pion electromagnetic form 
factor, enabling us to determine the two remaining parameters, g v and Am. This form 
factor is (in the time-like region) defined as the ratio of the cross section for e + e~ — > 7r + 7r~ 
divided by the cross section for the production of point-like pions. The amplitude for this 
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Figure 2.12: Contributions to the pion electromagnetic form factor in the 1/N C - expansion 
scheme. The propagator denoted by the curly line corresponds to the 1/N C - corrected rho- 
meson, while the double lines indicate RPA pions and sigmas. 



process is determined by the exchange of a (virtual) photon. Assuming vector meson 
dominance the coupling of the exchanged photon to the pions is mediated by a p-meson. 
This ansatz exhibits great success in describing the experimental data. It can also be 



simulated within the NJL model. The standard Hartree + RPA calculation [[71] for the 
corresponding amplitude, however, suffers from the completely unphysical nature of the 
intermediate p-meson, containing only quark-antiquark loops. We consider now two types 
of improvements. The first one is obvious: The RPA p-meson is replaced by a full one, 
containing mesonic, in particular pionic, fluctuations. The second one is more or less a 
question of consistency. In the standard scheme the photon couples to the p-meson via 
a quark-antiquark loop similar to the RPA polarization loop. Therefore we modify the 
coupling of the photon to the p-meson in the same way as we did for the polarization func- 
tions. The resulting diagrams for the l/iV c -expansion scheme are displayed in Fig. |2.12f 
In the MLA the diagram corresponding to 5Uy^ (fourth diagram in the lower line) is 
certainly missing. Note that the diagrams in the upper line would correspond to those in 
the standard scheme if we replaced the "improved" p-meson propagator (curly line) by 
an RPA one. 



Another type of possible improvement has been discussed in Ref. |[31|| : modifications of 



the quark triangle responsible for the coupling of the p-meson or the photon, respectively, 
to the pions. These would be suggestive if we considered full pions, i.e. including meson- 
loop effects. We decided, however, to take the external pions to be RPA ones. This is 
motivated by the following fact: it is more consistent with dressing the p-meson by RPA 
pions and in addition, the mass has been fitted to the empirical pion mass. 

Evaluating the diagrams of Fig. [2.12| for on-shell pions we notice that the resulting 
expression for the amplitude displays the same tensorial and isospin structure as that for 
point-like particles, i.e. it is proportional to (p M — p ,AI )e a fc c times the electric charge e. Here 
p and p' are the four-momenta of the external pions, a, b, and c are isospin indices of the 
pions and the photon. Calculating the ratio of the cross sections we become aware that 
this common factor and the common phase space factor cancel. We finally obtain for the 
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form factor, q = —p — p'^\ 



(2.22) 



where f(p,p') is a scalar function appearing in the p7r7r-vertex function (see App. [D]). Of 
course, the external pions are to be taken on-shell, i.e. p 2 = p' 2 = m 2 . an&p-p' = q 2 /2—m 2 r . 
Since a real photon (q 2 = 0) should "see" the actual charge of the pion, the form factor 
should equal unity at q 2 = 0. The subtraction performed in Eq. Q2.21Q guarantees that 
the p-meson self energy Yl p vanishes at this point. Therefore we have to check whether 
\9nqqf(PiP')\ 2 = 4 holds for q 2 = 0. This can most easily be shown with the help of the 
Ward-Takahashi identity derived in section |1.5| for the pirn- vertex (see Eq. ( |1.121| )). This 
identity is not directly applicable to our problem since the vertex is not contracted with 
q^ in our case. However, this can be handled by using the following relation, 



Using the Ward-Takahashi identity, Eq. ( |1.121| ), we can write the first term on the r.h.s. of 
this expression as the derivative of a difference of pion polarization functions in the RPA. 
Performing the limit q — > and comparing the result for the vertex with the definition of 
the pion quark coupling constant g^j we deduce 



<V(P,P')U = 2, (2-24) 

leading to the desired value for the form factor at q 2 = 0. 

Our results for \F V \ 2 in the l/iV c -expansion scheme as a function of the center-of-mass 
energy squared are displayed in the left panel of Fig. |2.13|, together with the experimental 



data f72|| . The three theoretical curves are calculated with different values of the meson 
cutoff, A M = 500 MeV (dashed line), A M = 600 MeV (solid), and A M = 700 MeV 
(dotted). We have already emphasized that these values for the meson cutoff lie within the 
small window of cutoff values where we can hope to produce reasonable results. The vector 
coupling constant g v is chosen such that the maximum of the form factor lies at the right 



position. The other parameters are taken from Table [2.1| . For = 500 MeV we observe 
a kink in the form factor slightly above the maximum ascribable to the quark- ant iquark 
threshold at s = 0.63 GeV 2 . Below the kink the form factor drops very steeply due to 
the sub-threshold attraction in the vector channel. Thus the poor description of the data 
above the maximum can be attributed to quark-antiquark threshold effects. Therefore 
we would conclude that our lower boundary for the meson cutoff, determined by the 
prescription that the threshold lies above the maximum, has not been chosen restrictively 
enough. However, the behavior of the form factor below the maximum, where threshold 
effects are less important, indicates that a higher value of the cutoff presumably produces 
better agreement with the data. Since increasing the cutoff reinforces the mesonic, in this 
case pionic, fluctuations we expect that with a larger Am the width of the form factor 
will be less underestimated. In fact, with a cutoff of Am = 700 MeV the width is already 

2 Note that in the corresponding formula in Ref. Efl] a factor of 1/2 is missing. 
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Figure 2.13: Pion electromagnetic form factor (left panel) and the nir-phase shifts in 
the vector-isovector channel (right panel) for Am = 500 MeV and g v = 1.0g s (dashed), 
Am = 600 MeV and g v = l.6g s (solid) as well as Am = 700 MeV and g v = 2Ag s (dotted) 
in the 1/N C - expansion scheme. The other parameter values are taken from Table The 
data points are taken from refs. /|7|/ and j7|/ ; respectively. 



noticeably overestimated coming along with a diminished height of the maximum. Below 
the maximum the description of the data is almost perfect with a cutoff of Am = 600 MeV. 
Above the maximum we certainly miss the detailed structure of the form factor around 
s = 0.61 GeV 2 which is due to p — cu-mixing, not included in our model. Approaching the 
threshold at s = 0.80 GeV , effects induced by quarks become more and more important 
resulting in an underestimation of the data in this region. Probably the overall agreement 
with the data can be somewhat improved by taking a slightly larger value for Am but 
since we are not interested in fine-tuning here, we keep the fit with Am = 600 MeV. 

Two quantities closely related to the pion electromagnetic form factor are well suited 
for confirming the consistency of our calculations with experimental data. The first one 
is the charge radius of the pion, which is defined as 

<^#U- (2 - 25) 

Obviously at q 2 = the details of the p-meson peak are less important than in the 
higher energy part. A simple pole ansatz for the form factor leads to (?" 2 ) 1//2 = 
^/6/m p = 0.63 fm ||74j| , a value which is already close to the experimental value of 
( r 2)i/2 = (o.663 ± 0.006) fm 0. Our results, (r 2 ) 1 / 2 = 0.59 fm for A M = 500 MeV, 



( r 2)i/2 = o.61 fm for A M = 600 MeV and (r 2 ) 1 / 2 = 0.66 fm for A M = 700 MeV, also 
satisfactorily agree with the experimental value. Since this quantity is not very well suited 
to fix the model parameters, which is clear from the above reasoning, we shall not worry 
about the fact that for A M = 700 MeV, which has been found to be too large to reproduce 
the data in the region of the maximum, the agreement with the data is almost perfect 
whereas for lower values of the cutoff the charge radius turns out to be slightly too small. 
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Figure 2.14: Diagrams contributing to the -scattering amplitude: Quark box diagram 
(left), t- and u-channel p- or a-meson exchange (middle) and s-channel p-meson exchange 
(right). Only the latter is taken into account in our calculations. 



Other observables are the 7T7r-phase shifts in the vector-isovector channel. We assume 
here that the dominant contribution to these phase-shifts arise from the s-channel p-meson 
exchange, visualized on the r.h.s. of Fig. |2.14j . The other processes shown in Fig. [2.14 



the direct scattering via quark box or the t- and u-channel exchange of a p- or a a- 
meson, respectively, are neglected. The result for different values of Am, together with 
the empirical data [73], is displayed in the right panel of Fig. 2.13| . Again, below the p- 
meson peak, which is to be understood as the location where the phase shifts equal 90 deg, 
the best description of the data is obtained with a cutoff of Am = 600 MeV, whereas for 
higher energies, where threshold effects begin to influence the results for the latter cutoff 
substantially, A M — 700 MeV produces better agreement with the data. 

The choice of contributions to the phase shifts has been made on purely phenomeno- 
logical grounds, motivated by the decisive role of the s-channel p-meson exchange. In 
addition to the s-channel p-meson exchange the direct scattering process via a quark box 
seems to be a suggestive contribution since it is of leading order in 1/N C , i.e. Hartree + 



RPA |f76| . Obviously it destroys unitarity of the ^-matrix if it is not iterated. Because 
of the non-trivial momentum dependence of the quark box diagram, this would be a very 
difficult task. On the other hand phase shifts are only properly defined if the S'-matrix 
fulfills unitarity. Our evaluation of the phase shifts in the vector-isovector benefits from 
the dominant role of the s-channel p-meson exchange which, when taken apart, is not in 
contradiction to unitarity. In any case, it can be checked that including the quark box 
diagram leaves the results for the scattering amplitude almost unchanged, thus justify- 
ing our assumption. In the scalar-isoscalar channel, on the contrary, we do not expect 
this pragmatic approach to be as successful. In that case the t- and -u-channel cr-meson 
exchange is essential in order to fulfill some of the low-energy theorems. Anyway, the 
a-meson within our model suffers from severe shortcomings which render it senseless to 



attack the problem of calculating phase shifts. This will be discussed in section 2.8 



The preceding presentation of the results in the l/A^-expansion scheme offered the 
possibility to decisively fix a consistent parameter set within this scheme. However, the 
much more narrow range for variations of the meson cutoff we found for the MLA in 
the previous section, renders it unlikely that such a fit is possible also in the MLA. In 
the l/A^-expansion scheme we found threshold effects to cause a steep drop in the pion 
electromagnetic form factor just below this threshold, leading to the conclusion that a 
constituent quark mass of about 400 MeV, shifting the threshold only slightly above 
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Figure 2.15: Real part of the p-meson polarization function il p as a function of the energy 
y/s in the rest frame of the meson. The solid line corresponds to the 1/N C - expansion 
scheme with Am = 600 MeV, the dashed line to the ML A with Am = 700 MeV. The other 
parameters are given in Table \2.1\ and Table [O, respectively. 



the p-meson peak in the form factor, is not sufficiently large. This further restricts the 
possible values for the meson cutoff to a small range just below Am = 700 MeV. Hence 
the only variable parameter to fit the form factor is the vector coupling constant g v . In 
fact, if we take A M = 700 MeV, it turns out that we unfortunately run into instabilities 
in the p-meson propagator similar to those found in the propagator of the pion (see 
section |2~5f) |||]. This observation is of course independent of the value of g v since it can 
be detected already in the polarization function Ii p . 

To this end the real part of the of the p-meson polarization function U p is plotted 
in Fig. |2.15| as a function of the energy \fs in the rest frame of the meson. The MLA 
result corresponds to the dashed line. For comparison we also show this function in the 
l/iV c -expansion scheme, using the 'best-fit parameters' given above (solid line). 

The "proper" behavior of the polarization function can be derived from what we 
expect of the inverse propagator. It has to develop a zero for \fs ~ m p to produce 
the well pronounced peak in the form factor. Thus 2g v Retl p has to become equal to 
unity for an energy of about 770 MeV. In the l/iV c -expansion scheme this can easily 
be achieved by an appropriate choice of the coupling constant since the polarization 
function is (at least below the quark- ant iquark threshold) a steadily rising function. In 
the MLA, on the contrary, this condition can by no means be fulfilled due to the maximum 
at \fs ~ 740 MeV and the subsequent steep drop the polarization function displays. 
Consequently for a weak coupling, i.e. small g v , we observe no zero at all in the inverse 
propagator whereas for stronger couplings we find two solutions, one at lower energies 
and one above the maximum. Obviously within those scenarios we would obtain only a 
very unrealistic description of the p-meson. Thus we conclude that we did not succeed in 
performing a fit of the observables in the pion sector, m^,/^, and (ipip), simultaneously 
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Figure 2.16: Real parts of the RPA contribution U p (solid) and the various correction 
terms toti p : 8H p a > (long -dashed), STly (dotted), 5Hy) (short- dashed), and8H p d ' (dashed- 
dotted). For all contributions we performed a subtraction, such that they vanish at ^/s = 0. 
The left panel corresponds to the 1/N C - expansion scheme, the right panel to the MLA. The 
model parameters are the same as in Fig. \2.15 . 



The cusp in the p-meson polarization function responsible for the peculiar behavior of 
the inverse p-propagator in the MLA looks very similar to that emerging in the inverse 



pion propagator for large values of the mesonic cutoff, see the discussion in section 2.5 



It has also the same origin: the strong momentum dependence of the only negative con- 
tribution, diagram 5Uf. This can be seen from the right panel of Fig. |2.16| , where we 
have separately plotted the various contributions to the polarization function of the p- 



meson. On the l.h.s. of Fig. [2.16| the same is shown for the 1/iVc-expansion scheme. 
One might wonder, why the results in the l/iV c -expansion scheme permit a reasonable 
description of the p-meson, although qualitatively the behavior of diagram 5Uy , respon- 
sible for the unphysical behavior in the MLA, agrees with that in the MLA. However, in 
the l/iV c -expansion scheme its contribution is almost cancelled by that of diagram 
which is not present in the MLA where the corresponding quark self-energy insertion is 
selfconsistently absorbed in the modified gap equation (Eq. ( |1.7U| )). Thus probably the 
purely perturbative treatment of both next-to-leading order quark self-energy insertions, 
<ffl (fe) and within the 1/N C -expansion scheme helps to suppress unphysical effects. 

However, we should emphasize that any unphysical behavior we detect here, although 
observed below the quark- ant iquark threshold, is mainly induced by quark effects. In sec- 
tion [2.5| we already noted that in particular the momentum dependence of the real part 
of 5U^ below the threshold is related via dispersion relations to an imaginary part with 
the "wrong" sign above the threshold. In the same way the real part of all other diagrams 
below the threshold is influenced by the quark-antiquark continuum. Only digram SU^ 
is not exclusively governed by effects which can be traced back to the quark-antiquark 
continuum: Here the imaginary part due to the two-meson intermediate state in a large 
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range of momenta dominates. In the following section we will carry on this discussion in 
connection with a comparison of these results with those obtained in the static limit. 

In any case the results suggest that the unphysical behavior of the polarization function 
can be at least weakened in energy regions relevant for the description of hadrons by 
pushing the constituent quark mass further up. This presumption is mainly supported by 
the steep drop in SU^ occurring only just below the threshold. One might suggest that 
the constituent quark mass is increased if further intermediate states are included in the 
model. However, we have seen that it is an extremely difficult task to implement p- and a±- 
mesons within the present model (see section |2.4.3|) . In addition the observations made in 
section |275| indicate that we have to cope with an even more severe problem when including 
vector and axial intermediate states: We observed "restoration" of chiral symmetry if the 
mesonic fluctuations became too strong. Thus, it seems as if the inclusion of the p- and 
ai-subspace for the intermediate states would not help without further modifications (or 
extensions) of the model. 



2.7 Comparison with the static limit 

This section is devoted to a comparison of the results obtained in the full model with 
those in the static limit, i.e. the approximation introduced in section [O]. This approx- 
imation is (in vacuum) equivalent to a purely hadronic description, i.e. all quark effects 
are suppressed. The only reminiscence of the full model can be found in the quark mass 
dependence of the effective coupling constants. Since the approximations to the MLA 
and the l/iV c -expansion scheme do not differ qualitatively from each other, it is com- 
pletely sufficient for our purpose to restrict the discussion within this section to one of 
the schemes. We have decided to take the 1/iVc-expansion scheme because only within 
that scheme we are able to consistently describe the data in the pion sector as well as for 
the p-meson. 

From the very beginning, we are for two reasons very distrustful of the static limit 
being a good approximation to the full calculation. First, the estimate of a realistic value 
of the constituent quark mass within this approximation (see Eq. ( |1.141|) ) seems to be in 
contradiction with the motivation of large quark masses as compared with the relevant 
momenta, with which we started in section |1.6| . In addition indications were found in the 
previous sections that quark effects are important for the mesonic polarization functions. 



The estimate performed in section |1.6| can also be applied vice versa: If we insert the 
parameters of our best fit to the pion electromagnetic form factor in the l/iV c -expansion 
scheme (those listed in Table ^J] for Am = 600 MeV) into Eq. ( |1.140| ), we obtain a value 



of 

9 pirn — 9.4 for the effective pirn coupling constant, which is considerably larger than 
the value of g pn7T ~ 6 usually needed in hadronic models to describe the data. Already 
with Am = 500 MeV, which in the full calculation noticeably underestimates the width 
of the p- meson, the value of g p7T7T ~ 6 is exceeded by about 40%. 

To compare the static limit with an exact treatment of the i/iV c -expansion scheme we 



will look at the imaginary part of the p-meson propagator which is displayed in Fig. [2.17 , 
The solid line represents a full calculation with Am = 600 MeV, the dashed line the 
corresponding static limit. Since we have fitted the parameters used in the full calculation 
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Figure 2.17: The imaginary part of the 1/N C - corrected p-meson propagator. The solid line 
represents the full calculation while the dashed line denotes the static limit. 



to the pion electromagnetic form factor, the maximum of the corresponding imaginary 
part is close to the p-meson mass of 770 MeV whereas this does not necessarily need 
to be the case in the static limit. In fact, the peak in the static limit is much broader 
and shifted to higher energies. Another obvious difference is the existence of the quark- 
antiquark threshold at \/s = 892 MeV in the full model which is not present in the static 
limit since within that approximation all quark effects are suppressed. 

For a better understanding of this intriguingly different behavior already below the 
unphysical threshold we also compare the real and imaginary parts of the corresponding 
self-energies (Fig. 2.18Q . Note that the imaginary part in this region is exclusively gener- 
ated by the decay process p — > tttt. Differences in the imaginary part therefore directly 
reflect the non-trivial momentum dependence of the quark loops in the meson-meson ver- 
tices. They become important only near the threshold. For \fs < 800 MeV the two curves 
almost coincide. In the real part the differences are much larger. Coincidence of the two 
curves can hardly be acknowledged up to y/s ~ 400 MeV. The real part in the static 
limit is obviously much less attractive than that in the full calculation. Besides, the slope 
of the curve is much steeper in the full calculation, which explains the broadening and 
the shift to higher energies of the peak in the p-meson propagator as shown in Fig. [2.17[ 
As already emphasized several times, the real part below the threshold is, via dispersion 
relations, related to the entire imaginary part including that above the threshold. Hence 
the discrepancies between the static limit and the full model here not only reflect the mo- 
mentum dependence of the quark loops contained in the vertices but additionally indirect 
effects arising from quark decay channels. 

The only conclusion we can draw from the preceeding discussion is that quark ef- 
fects are indeed important. At this point one might argue that one has to distinguish 
between "physical quark effects", for instance those modifying the point-like structure 
of the meson-meson vertices, and unphysical quark effects leading to unconfmed quarks. 
However, these two types of effects cannot be considered separately since they are linked 
to each other via dispersion relations. This is related to a rather fundamental question, 
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Figure 2.18: The real part (left) and the imaginary part (right) of the 1/N C - corrected p- 
meson self-energy, multiplied by —2g v . The solid line represents the full calculation, while 
the dashed line indicates the static limit. 



namely whether "physical quark effects" exist at all. Of course, at sufficiently high en- 
ergies quarks should be the relevant degrees of freedom which in one way or another 
influence hadronic properties. An approach like QCD sum rules [77] in principle offers 
the possibility to map at least different energy regions. But in any case an explanation of 
this influence requires first the understanding of quark confinement itself. 

Our scheme can only exemplify how a hadronic model emerges naturally from the 
underlying quark structure provided that the necessary degrees of freedom are included. 
This can be seen nicely from the low energy region where the static limit almost coincides 
with the exact treatment of our model. 



2.8 Results for the sigma- meson 

The sigma-meson is not a very well established particle in nature. In the "Review of 
Particle Physics" [50| we read concerning the a: "The interpretation of this entry as a 
particle is controversial." The mass of this "particle" lies between 400 and 1200 MeV with 
a width of the same magnitude, 600-1000 MeV. However, in chiral models where chiral 
symmetry is linearly realized, like the NJL model, it inevitably accompanies the pion as its 
chiral partner. Most of the models deal with a sharp a-particle on tree-level. However, due 
to its strong coupling to a two-pion state the cx-meson is significantly broadened. Within a 
linear sigma model approach this notion has been successfully applied to predict not only 
a broad resonance in vacuum but also its sharpening due to partial restoration of chiral 
symmetry in medium (see Schuck et al. |78| and references therein). The description of 



the u-meson in the NJL model in the standard Hartree + RPA scheme is certainly not 
suitable to obviate this problem. Within that approximation the cr-meson is an (almost) 
sharp particle, located at the gg-threshold or just above. Thus, the inclusion of mesonic 
fluctuations, in particular pionic ones, in the l/iV c -expansion scheme or the MLA seems 
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Figure 2.19: The imaginary part of the "improved " a-meson propagator, D a in the 1/N C - 
expansion scheme (left) and in the ML A (right) for Km — 600 MeV. 



to be a promising starting point to improve the phenomenology of the cr-meson within 
the NJL model. 

Actually, in both schemes, the l/7V c -expansion scheme and the ML A, the cr-meson 
is substantially broadened above the two-pion threshold as can be inferred from the left 
and right panel of Fig. [2.19| , respectively, where the imaginary part of the cr-propagator 
D a (s) as a function of the energy y/s in the rest frame of the meson is displayed. For 
the 1/iVc-expansion scheme the calculation has been performed with Am = 600 MeV and 
the corresponding "best fit" parameter set listed in Table |2.1| . In the MLA, where we did 
not succeed in performing a consistent fit, we also took A^ = 600 MeV together with 
the parameters shown for that value of the cutoff in Table |2.2| . In the MLA we find, in 
addition to the broadening above the two-pion threshold, a sharp peak just below that 
threshold which would be a ^-function if we had not introduced by hand some small but 
nonzero imaginary part. The origin of this peak will become clear from the behavior of 
the polarization function itself. 

The left panels of Fig. |2.20| show the real parts of the different contributions to the 
cr-meson polarization function as a function of s. The momentum dependence of the real 
part is mainly caused by the contribution arising from SU^ (long-dashed line), i.e. the 
correction term containing the two-pion intermediate state. Below the gg-threshold this 
term is in fact the only one producing a nonvanishing imaginary part. This imaginary 
part is shown on the r.h.s. of Fig. [2.20| . These features are common for the l/A^-expansion 
scheme and the MLA. 

This rather encouraging result, at least for the l/A^-expansion scheme, is spoiled, 
however, if we investigate the analytic properties of the propagator more closely. Recalling 
that the propagator acquires a pole if 2g s Yl a (s) becomes equal to unity, one understands 
the existence of the sharp peak in the cr-meson propagator in MLA. In the l/A^-expansion 
scheme, on the other hand, the existence of a pole is not as obvious. In the time-like region 
we find no pole at all, explaining the smaller scale of the imaginary part of the propagator 
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Figure 2.20: The real part (left) and the imaginary part (right) of the various correc- 
tions terms to the a-meson polarization function, multiplied by 2g s , in the 1/N C - expansion 
scheme (upper panel) and in the MLA (lower panel) : total polarization function (solid 
line), 5U^ (long dashed line), STL^ (dotted line), 511$ (short dashed line), and SU^ 
(dashed- dotted). Note that the contribution of 5U^ is missing in the MLA. 



as compared with that in the MLA (Fig. |2.19|) . We only guess from the slope of the 
polarization function at s = that this might very well occur for negative, i.e. space-like, 
s. This would of course be a very disturbing unphysical feature of the cx-propagator. 
Unfortunately we cannot support this guess by explicit results: Since we work in the rest 
frame of the mesons, we are not able to perform calculations for space-like momenta. 

However, another fact which can already be deduced from the value of the real part of 
the polarization function at s = corroborates the doubts concerning our results for the 
ex-meson: The sign of the inverse propagator at s = differs from what we except from a 
comparison with a free propagator l/(s — m 2 a ). On the l.h.s. of Fig. [2.21| we have displayed 
the inverse cr-propagator as it emerges from the l/iV c -expansion scheme (solid line), in the 
static limit (short-dashed line), the free one with an arbitrarily chosen mass oim a = 2m# 
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(dotted line) and for comparison the RPA one (long-dashed line) and the result obtained 
by Aouissat et al. [[fjj (dashed-dotted line). Note that we have normalized the results in 
RPA and in the l/^-expansion scheme to the RPA a-quark coupling constant g^jq. in 
the vicinity of s=0 the slope of the various inverse propagators agrees at least concerning 
the sign, while the function itself has positive values in the l/iV c -expansion scheme and 
negative ones in the two other cases. This behavior is basically not altered if we vary 
the mesonic cutoff. This can be inferred from the right panel of Fig. |2.21| where we have 
plotted the inverse cx-propagator in the 1/iVc-expansion scheme using different values of 
Am while all other parameters are kept constant. The effect of varying the cutoff is 
mainly to shift the entire curve up or down. The shape remains almost the same. For 
small values of Am, cf. the dashed line with Am = 200 MeV, we encounter the even more 
peculiar situation, which has already been observed in the MLA (see Fig. |2.20| ), where 
we have two zeros in the real part of the inverse cr-propagator in the time-like region, one 
below the two-pion threshold at about 280 MeV and one above. With increased cutoff 
the curve is first shifted upwards, then, if the cutoff exceeds Am > 900 MeV it begins 
to move downwards again. Nevertheless it seems to be hopeless that we can obtain a 
more realistic description of the a-propagator for very large values of the cutoff. On the 
one hand, if the shape of the curve is not altered, we cannot avoid finding finally two 
poles in the time-like region instead of one. On the other hand the investigations in the 
previous sections, concerning the pion and the p-meson, render it plausible that we run 
into instabilities when applying such large values for the cutoff, besides completely failing 
to describe the data in the pion and p-meson sector. 

Till now each unphysical phenomenon we detected could be ascribed to strong quark 




Figure 2.21: Real part of the inverse a-meson propagator. (Left): g^D~ 1 (s) in the 1/N C - 
expansion scheme with A M = 600 MeV (solid line) and in the static limit (short-dashed 
line), g^ ( j q 2 D^ 1 (s) in RPA (long-dashed line), s — Am 2 H (dotted line), and the corresponding 
linear sigma-model result by Aouissat et al. fiTQI (dashed-dotted line). (Right): 1/N C - 
expansion scheme with Am = 600 MeV (dashed), Am = 600 MeV (solid), and Am = 900 
MeV (dotted). 
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Figure 2.22: Example for contributions to the a-meson polarization function similar to 
those included in Ref. ^7!% (see text). 



effects. This is not true for the present case. The left panel of Fig. |2.21| also displays the 
result in the static limit (short-dashed line), where all quark effects are suppressed, which 
shows essentially the same behavior, at least in the low-energy region, as the calculation 
in the l/iV c -expansion scheme. This leads us to the conclusion that the strong attraction 
of the two-pion state in the scalar-isoscalar channel is responsible for this behavior. It 
is therefore suggestive that it can be mended by imitating the terms which produce 
reasonable results in a linear sigma-model calculation by Aouissat et al. |7{| (see l.h.s. 

consists of 



of Fig. [2.21] , dashed-dotted line). The RPA used by the authors of Ref. |T5| 
iterating the entire pion-pion scattering amplitude in the scalar-isoscalar channel for the 
self-energy of the cx-meson. In our case this would mainly correspond to terms like those 
presented in Fig. |2.22|. Terms of that type would have been included in the "improved" 



polarization function if we had applied the ^-derivable method discussed in section 1.3.3[ 



The scattering kernel of the quark-antiquark Bethe-Salpeter equation visualized in Fig. 1.8 



precisely contains those non-local contributions necessary to generate the desired cx-meson 
polarization function. In any case, we resigned to pursue that scheme further because it 
is very difficult to handle due to the non-local structure, which already shows up in the 
modified gap equation (cf. Fig. 01). 
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Chapter 3 

Results at nonzero temperature 



It is commonly believed that chiral symmetry, which is spontaneously broken in vacuum, 
gets restored at high temperatures as well as at high densities. For a recent review on 



the present view of the QCD phase diagram see Ref . . In wide ranges of densities and 



temperatures this notion is, however, based on model calculations for lack of fundamental 
knowledge. Only at low temperatures and low densities model independent results can 
be obtained by considering a gas of pions or a gas of pions and nucleons, respectively. 
Pions are the dominant degrees of freedom in that range because of their small mass 
as compared with other possible degrees of freedom. At zero chemical potential, i.e. 
zero density, lattice calculations provide us with results also in the vicinity of the phase 
transition, but at nonzero chemical potential no realistic lattice results exist. Thus in 
that region we have to rely entirely on model calculations. 

One of the models which has been applied by many authors to study the effects of 
(partial) chiral symmetry restoration at nonzero temperature and density on various quan- 
tities, in particular the quark condensate as the order parameter of chiral symmetry and 
properties of mesons, is the NJL model. Most of these investigations have been performed 



in the standard approximation scheme, i.e. Hartree + RPA (see e.g. [IT], [18], [19], ^0|). The 
main drawback of these examinations is certainly the unphysical nature of the degrees of 
freedom, deconfmed quarks, taken into account to describe the thermodynamical proper- 
ties of the system. The presence of unphysical degrees of freedom is not the only weak 
point of these studies. In addition, the relevant degrees of freedom, at low temperatures 
and densities mainly pions, are completely missing such that, for instance, the model in- 
dependent results within that range cannot be reproduced by these calculations. Hence, 
although the fundamental problem of lack of confinement cannot be overcome, we can 
hope to improve the results by introducing mesonic degrees of freedom within an approx- 
imation beyond Hartree + RPA, i.e. in the l/iV c -expansion scheme or the MLA. Since 
both schemes do not contain any nucleons, essential ingredients at nonzero densities, we 
will restrict our examinations in the l/iV c -expansion scheme and the MLA to nonzero 
temperatures but zero density. 
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3.1 Results in Hartree approximation + RPA 



Before we present first results at nonzero temperature in the 1 /^-expansion scheme 
and the MLA, we will briefly review the results in the standard approximation scheme, 
namely Hartree + RPA. To this end we first have to generalize the formalism presented in 
Chapter [l] to nonzero temperature and chemical potential. We will keep the latter since it 
enters into the expressions in almost the same way as nonzero temperature. In principle 
we can choose between two different methods to deal with thermodynamic properties 
of the system: the imaginary time or Matsubara formalism (see e.g. |j8~0|), and the real 



time formalism (see e.g. |8l|). Within this paper we will adopt the Matsubara formalism 
for simplicity. We should remark that it is not applicable to determine the temperature 
dependence of the mesonic properties in the l/iV c -expansion scheme or the MLA. This 
requires the analytic continuation of functions only known numerically for a finite number 
of points, which cannot be done uniquely. In any case we will restrict the later discussion 
on the quark condensate which is a static quantity and can therefore be treated within 
the Matsubara formalism. 

The basics of the Matsubara formalism can be found in textbooks, e.g. in Ref. [ jB0|| . 
In App. [E] we summarize the Feynman rules we obtain starting from the Lagrangian in 
Eq. fll.lj) . The resulting expressions at nonzero temperature and chemical potential p, can 
be written formally in almost the same way as in vacuum if we replace the integration 
over energy in the vacuum expressions by a sum over fermionic or bosonic Matsubara 
frequencies iu n , i.e. 

If we have a nonvanishing chemical potential for the quarks, p should be added to the 
fermionic Matsubara frequency in all expressions, i.e. iu n — > iu n + p. The replacement 
prescription in Eq. ( p.l|) can be applied to define the analogue to the elementary integrals 
in vacuum, e.g. 



I{p) 



d 4 k 1 



I(iuji,p) = iTj2 



(2tt) 4 (k 2 -m 2 + ie)((k + p) 2 - m 2 + ie) 
d 3 k 1 



(2vr) 3 ((^ n + Ai )2 _ k 2 _ m 2)^ iuJn + fl + i Ul )2 _ { J f+ £)2 _ m 2) 

(3-2) 



Note that the expressions at nonzero T and p depend separately on energy and three- 
momentum because a preferred frame of reference exists in the heat bath. This simplifies 
our notation: the second argument enables us to clearly distinguish vacuum quantities 
from their medium analogues. Functions with one argument, e.g. p, denote vacuum 
quantities, whereas functions with two arguments, e.g. iuii,p, are the corresponding in- 
medium quantities, such that we can use the same symbols to mark analogous quantities. 
The only exception are static quantities which do not depend on momentum at all. There 
we will add a suffix T, e.g. the quark condensate in medium will be denoted by (ip^T- 
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Figure 3.1: (solid line) and mf^ (dashed) as a function of temperature in the chiral 
limit (left panel) and with m = 6.13 MeV (right panel). For the remaining parameters 
see Table fj| for A M = 0. 



The proportionality of the constituent quark mass to the quark condensate in the 
Hartree approximation, see Eq. ( |1.6| ), suggests that this quark mass is strongly influenced 
by medium effects in the vicinity of the phase transition. In fact, the medium analogue 
of the gap-equation, Eq. (|1.2[), 



rriHT 



m 



2g s AN c N f m HT Tj2 



d 3 k 



1 



(27r)3(^ n + /x) 2 -£ 2 



m + 2g s AN c N f m HT I 1T 



(3.3) 



with E = y k 2 + m HT and u n = (2n+l)irT being fermionic Matsubara frequencies, reveals 
the T as well as the [i dependence of the constituent quark mass tuht- The self-consistent 
solutions of this equation correspond to extrema of the thermodynamic potential 

^( % ,r,,) = /^-ln( E2 Z^ + ^ ) + ^- m ^ +C onst. , 

(3.4) 

which can be employed to study the phase structure of the system. In the above formula 
we set g v — for simplicity. For details how to include vector interactions see e.g. 

In this approximation the NJL model exhibits a second order phase transition 

200 MeV and, depending strongly 

0,/i ^ 



Ref . § 

at /i = 0, T ^ with a critical temperature of T c « 150- 

on the choice of parameters fl8~3| , a second or a first order phase transition at T 
with a critical chemical potential of the same order of magnitude as the constituent quark 
mass in vacuum, /i c rs 300 — 400 MeV. 

In addition, the properties of the mesons in the RPA are modified by medium effects. 
The polarization functions for the RPA pion and the cx-meson read 



U^i^p) = 4iN c N f I 1T - 2iN c N f {(iu l f - f - Am 2 T )I(iuo u p) 
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Figure 3.2: Imaginary part of the p-meson (solid line) and a\-meson (dashed line) prop- 
agator in RPA for different temperatures. 



Yi^uP) = AiN c N f I 1T -2iN c N f {{iuo l f-f)I{iuo h p), (3.5) 

with uoi = 2lnT being bosonic Matsubara frequencies. Below the phase transition the in- 
tegral Iit can again be replaced with the help of the gap equation Eq. ( |3.3| ) (cf. Eqs. ( [B.2| ) 
and ( B.3|) ). The p-meson and the ai-meson are a little more difficult to handle since we 
have to keep in mind, that, although the four-dimensional tranverse structure of the cor- 
responding polarization functions is preserved in medium, we have to distinguish between 
a three-dimensional longitudinal and a three-dimensional transverse part. This will be 
explained in more detail in App. [C|. 

Although the thermodynamics within this approximation is exclusively driven by un- 
physical degrees of freedom, deconfined quarks, some features of chiral symmetry restora- 
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tion can nicely be illustrated. One example is the degeneracy of the spectral functions of 
the chiral partners, i.e. ir and a as well as p and a\, above the phase transition can be 
mentioned here. For the 7r-<r-system this can be inferred from Fig. |3.1| where the masses 
of pion and cr-meson [] as a function of temperature for p = are displayed. For param- 
eters we chose our standard RPA parameter set. This leads to a critical temperature of 
T c = 164.4 MeV. In the chiral limit (left panel) the mass of the pion, denoted by the 
solid line, vanishes below the phase transition. Away from the chiral limit the pion mass, 
although non-vanishing, remains almost constant up to T w T c . In both cases the mass 
of the a-meson (dashed line) decreases with temperature, approaching zero at T c in the 
chiral limit, and increasing again above the phase transition. Away from the chiral limit 
pion and a mass do not exactly coincide at (and above) T c . For the p-ai-system it is not 
sufficient to consider masses alone since both, the p- and the ai-meson, lie well above the 
gg-threshold already for temperatures T <C T c . In Fig. |T2]we have therefore displayed the 
imaginary part of the propagators for T = 0, T = .61T C and T = .99T C . Obviously these 
imaginary parts coincide for T = T c in the chiral limit. We have to emphasize again that 
these results can only serve as an illustration. Among other points the observation that 
the p- and the ai-meson lie above the gg-threshold over wide temperature ranges reveals 
the unphysical nature of these results. 



3.2 Quark condensate at T ^ beyond mean field 

Within this section we will study the behavior of the quark condensate at nonzero tem- 
perature (and p = 0) in the l/iV c -expansion scheme and the MLA. At the beginning of 
chapter ^ we expressed the hope that at least the low-temperature behavior can be im- 
proved considerably as compared with the Hartree + RPA scheme by including mesonic 
fluctuations within the l/iV c -expansion scheme or the MLA. To support this we will first 
compare the low-temperature behavior we obtain with the expectations from model in- 
dependent considerations. Within this section we will again neglect fluctuations of p- or 
ai-mesons. 



3.2.1 Low-temperature behavior 

Model independent results for the quark condensate can be obtained by a strict low- 
temperature expansion in chiral perturbation theory. In the chiral limit and at vanishing 



baryon density this reads 44 



(^>T=<W(l-^-g^ + ...) • (3-6) 

Here (V>V) denotes the quark condensate at zero temperature. The term proportional to 
T 2 arises from a pure (massless) pion gas. It can be derived by considering the pressure 
of a free pion gas, 

P. = -3T/^ln(l-e-^), (3.7) 
1 We define masses here as the zero in the real part of the inverse propagator of the meson at rest. 
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where E w = \j k 2 + m 2 . The quark condensate is in general given as the negative of the 
derivative of the total pressure with respect to the current quark mass. Assuming that 
for small T the change of the quark condensate with temperature is determined by the 
pion gas and using the GOR relation, Eq. (|1.78|) , we can express the quark condensate as 
follows: 



(V>V)t = (W>) + S(W)t = (#) + 



(0-0) dp n 
fl dm* 



7 (00)3 r d 3 k 1 1 
= (#) -T2iW l-e-^ ' (3 ' 8) 

In the chiral limit the integral in the last line of the above expression can be evaluated 
analytically to give the T 2 -term in Eq. fl3.6| ). The higher-order terms in addition contain 
interactions between the pions. In the chiral limit it is not astonishing that the massless 
Goldstone bosons, the pions, are the dominant degrees of freedom at low temperatures. 
Thus we expect at least the T 2 -term to reflect only the underlying chiral symmetry and 



not any model specific contributions. However, it can be shown |44| that aside from the 
T 2 - also the T 4 -term of this expansion follows from chiral symmetry alone. This implies 
that every model which exhibits chiral symmetry, among others the NJL model, should 
in principle show the same behavior up to T 4 . In contrast an expansion of the quark 
condensate in Hartree approximation [pO, 



+ (3-9) 



{2m H Tfl 2 _^ 

7T 3 / 2 (#)(°) r 

obviously does not reproduce these terms. There exists a very simple argument |45| why 
the Hartree approximation is not suited to reproduce these terms: From the fact that 
/„- is of the order \/W c we gather that the T 2 -term is of the order 1/N C and the T 4 -term 
of the order 1/N 2 . Remembering that the Hartree approximation corresponds to the 
leading order in an l/iV c -expansion we immediately see that this approximation must fail 
in the desired low-temperature behavior. This reasoning also reveals that by extending 
the calculations to next-to-leading order in l/N c we will at least be able to reproduce 
the T 2 -term [fl2 |. Although the MLA does not correspond to an expansion in orders of 



1/N C the T 2 -term can also be reproduced within that approximation scheme because 
it contains all l/iV c -correction terms in next-to- leading order. We will now proceed with 
demonstrating this first for the 1 /Ac-expansion scheme, and then for the MLA. 

Using the notations introduced in the previous section for the temperature dependent 
quantities, we can write the quark condensate in next-to-leading order in 1/N C in almost 
the same way as in vacuum (see Eqs. ( p..6|) and (|1.32p . 



/7/\ /7 /\(0) ,jf/7 /\ m HT -m D a (0, 0) A T 
(#) r =(#)V+W>T = 2g ' ^ ' 

From Eq. ( |3.9| ) we conclude that to order T 2 thermal effects in the leading-order term 
(■0'0)( ) can be neglected. In the same way one can reason that the temperature depen- 
dence of the cr-propagator does not contribute. Thus we are left with an expansion of A T 
at low temperatures. is explicitly given by 
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Figure 3.3: Contours C , C and C" for the evaluation of At- The poles at the Matsubara 
frequencies and at the pion energy are explicitly marked as well as the quark- antiquark 
cuts. 

d 3 p 



A T = AiN^fmrTjjJL^ 



D a {ioo h p)(2I{ioo h ft) + 7(0,0) - {{luutf -P 2 ~ ^m 2 T )K{iuo u p)) 
+D^iu h p) (3/(0,0) - 3((^) 2 -/) K(iu h p) )} , (3.11) 

where lo\ are bosonic Matsubara frequencies. If standard techniques are used the sum 
over the Matsubara frequencies in Eq. ( |3.11| ) can be converted into a contour integral [SO 



where C is a contour encircling the imaginary axis in the positive sense, see Fig. ^3 



At = 4iN c N f ™t±J ^ j cl ^{D^ (3/(0,0) - 3(z 2 - f) K(z,p)) 

+D a (z,p) (2/( 2 ,^ + /(0,0)-(^ 2 -p 2 -4m|) K(z,p})} . (3.12) 

One proceeds by deforming the contour C, first to C, then to C", in such a way that only 
the poles of the addends in Eq. ( p.llj) contribute. In our case we encounter poles at 



z = ±\Jp i + mT ) = ±E n 

and cuts for 



z < — y Am 2 HT + P 2 an d z > y Am 2 HT + p 2 . 
The resulting contour C" is also shown in Fig. |3l| In the evaluation this contour integra- 
tion each pole at z = z p is weighted with a factor l/(e Zp / T — 1). Thus at low temperatures 
the lowest-lying pion pole at z = ±E m gives the main contribution. Subtracting the 
vacuum part it can be expressed in the following way: 

(A T — A)^ = -AiN c Nf m J ( ^ 3eE J_ 1 Res{^(^,^(3/(0,0)-3(/; 2 -/)K(^,p))} 

(3.13) 
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The exponential suppression of the other contributions is certainly the stronger the lower 
the pion mass is as compared with two times the constituent quark mass rriHT- in the chiral 
limit we can therefore approximate Ay for low temperatures by the pionic contribution 

A T -A = 4i Vo i V/m / I g 5iH7 |- T {^i^-}. (3.14) 

This integral can be evaluated analytically and we obtain 

A T -A = m- . (3.15) 

Up to now we have shown that the changes of the quark condensate for low tempera- 
tures are indeed proportional to T 2 in the l/iV c -expansion scheme. We are left with a 
comparison of the coefficients in Eqs. ( |3.6| ) and ( |3.10| ), respectively. In the chiral limit 
the vacuum cr-meson propagator can be expressed with the help of the leading-order pion 
decay constant as, (cf. Eqs. ( |B.4j ) and fl2.1U| )), 

Ar(O) = --^y • (3.16) 
4/vr 



Inserting this expression together with Eq. ( |1.6|) for the leading order quark condensate 
and Eq. ( 13.15 ), evaluated at the Hartree mass, into Eq. ( |3.10|) we finally arrive at the 



following expression for the quark condensate in next-to-leading order in 1/N C at low 
temperatures: 

(^) T = (^)-(^) (0) -iy • (3.17) 

8/tt 

A comparison with Eq. ( |3.6| ) shows that the coefficients in front of the T 2 -term are similar 
but not identical. In the l/iV c -expansion scheme the leading-order quantities, (ipip)^ ' and 
determine the coefficient in contrast to (V>V) an d fn i n Eq- ( |3.6|) . On the one hand 
this simply reflects the fact that we only get the lowest order contribution to this coef- 
ficient by expanding the quark condensate up to next-to-leading order in 1/N C . On the 
other hand it can be understood from physical arguments: The l/iV c -corrections to the 
quark condensate consist of fluctuating RPA mesons. Hence the low temperature behav- 
ior within the 1 /iV c -expansion scheme is driven by excited RPA pions and consequently 
quantities corresponding to the RPA pions. 



For the MLA, a similar result has been derived in Ref. |]4"5" |. With the preceding results 
this derivation is easily comprehensible. In the chiral limit the quark condensate in the 
MLA, cf. Eq. ( |1 . 73| ) , is proportional to the quark mass m' T , 

(W>) = , (3-18) 



where m' T is a solution of the modified gap equation, cf.Eq. ( |1.70| ), 



m' T = mo + J2ht(itT'' t ) — 2g s A.T{m' T ) . (3.19) 
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This equation implicitly defines the temperature dependence of m' T and therefore also that 
of the quark condensate. To obtain the Independence we will now take the derivative of 
Eq. ( |3.19| ) with respect to T 2 at T = and solve this for dm' T /dT 2 



1. e. 



dm'n 



dT 2 



<9(£ 



HT 



T 



T=0 



dT 2 

m! , 
-2^ T (1 - 



T=0 



(1 



d(Y, H -2gA),-x 



dm' 



d(E H -2g a A) y 1 
dm' 



(3.20) 



where in the last step we have employed Eq. (|3.15| ) and used the fact that the derivative 
of Tjht with respect to T 2 vanishes for T = 0. From the definition of the inverse meson 
propagators in the MLA in section |1.3.3j it is obvious that the derivative of the quark self 
energy with respect to the quark mass is related to the a-meson polarization function 
n <J (0). Thus the above expression can be written in terms of the a- meson propagator for 
vanishing momentum 

dm'rr ~ , , ml 

(3.21) 



dT 2 



T=0 



The authors of Ref. ^§1 argue that the leading order piece in 1/N C of the above expression 
is the RPA cr-meson propagator which in turn can be expressed with the pion decay 
constant to that order. Consequently the authors find 



(#>r = 



r 2 



8/. 



2(0) 



(3.22) 



We have to keep in mind that f 2 ^ has to be understood as the RPA-pion decay constant, 
Eq. ( 2.10| ), evaluated at the quark mass ml . This result is consistent with the fact that 
in the MLA the pionic excitations, which determine the thermal corrections to the quark 
condensate at low temperatures, consist of RPA pions built of quarks with the constituent 
quark mass m! . Since (?/>?/>) ^ evaluated at the new mass m', agrees with (^>ip), {^>i^)t to 
order T 2 is here, in contrast to the l/iV c -expansion scheme, proportional to the "improved" 
quark condensate (ipip). 

Of course, the l/iV c -expansion argument, leading from Eq. ( |3.21| ) to Eq. ( |3.22j ), is 
questionable because already the constituent quark mass contains arbitrary orders in 
l/N c . Nevertheless the result, Eq. ( |3.22| ), seems well justified on physical grounds, since 
the dominant degrees of freedom are in any case the excited RPA pions. Together with 
the discussion of our numerical results in section |3.2.3| we will come back to the question 
whether this approximation is well justified. 



3.2.2 Numerical results within the 1/N C - expansion scheme 

Our numerical results for the temperature behavior of the quark condensate within the 
1 /A^-expansion scheme are displayed in Fig. |3.4j . For the r.h.s. we have used the pa- 
rameters listed in Table |2.1| for Am = 600 MeV which allows a reasonable description 
of data in vacuum. In particular the RPA pion mass reproduces the empirical value of 
m (o) _ M e v. The l.h.s. of Fig. [3]4| corresponds to the chiral limit which has been built 



as described in section |2.5| starting from the above parameter set. 
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We will begin the discussion by the results in the chiral limit. Up to temperatures of 
T pa 100 MeV the next-to-leading order result almost perfectly agrees with the pure pion 
gas result whereas the leading order result remains almost constant. We have to admit, 
that the denotation "pure pion gas" result has to be considered carefully. Motivated by 
the low-temperature expansion in the previous subsection, see Eq. ( |3.17| ), we compare 
our result here with a gas of free leading-order pions. The applied parameter set, how- 
ever, has been fixed by fitting the l/iV c -corrected quantities (except the pion mass, cf. 
section |2.6.1| ), hence the value of f(°\ determining the behavior of the quark condensate 
at low temperatures, does not coincide with the empirical value. The deviation is of 
about 25%. Nevertheless we can state that the l/iV c -corrections lead to a considerable 
improvement as compared with the leading order result since our comparison reveals the 
dominant role of pionic excitations at least in the low temperature regime, where quark 
degrees of freedom are exponentially suppressed. 

However, we cannot by-pass the fundamental problem of lack of confinement in the 
NJL model. The larger degeneracy factor of the quarks (24 as compared with 3) will 
consequently favor thermally excited quarks at higher temperatures. To obtain a rough 
estimate of the temperature range where quark effects become important let us compare 
the pressure of an ideal quark gas, 



Pq = 4N c N fT J^-Hl + e' E ^) 



(27T) 1 



(3.23) 



E q = y k 2 + m 2 , with that of pions, see Eq. ( |3.7| ). For the case of massless pions, i.e. in 
the chiral limit, the latter can be evaluated analytically to obtain 



Pit 



T 



30 



(3.24) 
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Figure 3.4: Quark condensate as a function of temperature, normalized to the vacuum 
value, in the chiral limit (left) and with = 140 MeV (right). Leading order in l/N c 
(dashed line), next-to-leading order (solid line), and free pion gas (dotted line). 
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Figure 3.5: Pressure of a free gas of quarks with a mass of m = AA6MeV (solid line) 
compared with that of a gas of massless pions (dashed line) as a function of temperature. 



In Fig. |3.5| the solid line represents the pressure of a gas of free quarks with mass m = 
446 MeV, the dashed line that of a gas of massless pions. The temperature of roughly 
100 MeV where first deviations of the results in the l/iV c -expansion scheme from the pure 
pion gas result become visible lies in the same range where the quark pressure becomes 
compatible with the pion pressure in the free gas approximation. 

At this point we enter into the question about the physical meaning of quark effects 
at these temperatures. The excitation of quark degrees of freedom is not completely 
excluded by nature. At temperatures above the deconfinement phase transition quarks 
are, besides gluons, the principal degrees of freedom. In the NJL model, however, there 
is no confinement and consequently no deconfinement transition. On the other hand 
lattice calculations ]84]] indicate that the deconfinement transition coincides with the 
chiral phase transition which can be examined within the framework of the NJL model. 
Thus the critical temperature for the chiral phase transition determines the scale for the 
judgement on the relevance of quark effects. In principle one might take the position that 
quark effects acquire a physical meaning only above the phase transition. But as long 
as we are ignorant about the detailed dynamics of the deconfinement transition we could 
allow for quark degrees of freedom already in the vicinity of the phase transition. From 
a phenomenological point of view this is consistent with the finding that a resonance gas 
with many degrees of freedom can be, close to the phase transition, effectively modeled 
by a quark gas ( "quark-hadron duality" ) . 

Unfortunately, in the 1/A^-expansion scheme the perturbative treatment of the 
mesonic fluctuations prevents us from investigating a possible phase transition and con- 
sequently from estimating whether thermally excited quarks become important near the 
phase transition or much below. Anyway, probably the applicability of the perturbative 
expansion scheme has to be considered with great care already for much lower tempera- 
tures than T c . We will come back to the question about the relevance of thermally excited 
quarks in the next section in connection with the ML A. 
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Figure 3.6: Quark condensate in the chiral limit as a function of temperature, normalized 
to the vacuum value, Hartree approximation (dashed), meson loop approximation (solid) 
and free pion gas (dotted). 



This section will be completed by a discussion of the results away from the chiral limit, 
i.e. with mo 7^ 0, shown on the r.h.s. of Fig. 3.4 . The fundamental difference to the chiral 
limit is the presence of massive pions which are now, like all other particles, exponentially 
suppressed. This leads to the much more flat behavior of the quark condensate at low 
temperatures. Although the pions are massive they are still the lightest particles and 
can therefore be excited most easily. This is confirmed by comparing the pure pion gas 
result (dashed line) with the next-to-leading order result (solid): For temperatures below 
T < 100 MeV the two curves satisfactorily agree, though less perfectly than in the chiral 
limit. The former has been calculated according to Eq. ( p.8|) with a nonzero pion mass. 
Stronger deviations of the next-to-leading order result from the pure pion gas result above 
T fa 100 MeV indicate the increasing importance of quark effects, which become visible 
at almost the same temperature as in the chiral limit. 



3.2.3 Meson loop approximation 



Let us now compare the results of the previous subsection with the analogous calculations 
in the MLA. A study of the temperature dependence of the quark condensate within the 
ML A was first performed in Ref. |45j and can, for the chiral limit, also be found in 
Ref. 



12 



Here we likewise restrict ourselves to the chiral limit because our main aims 
are on the one hand to be able to judge the relevance of quark effects on more sensible 
grounds than in the l/iV c -expansion scheme and, on the other hand, to test the quality 
of the approximation which led us to Eq. (|3.22|) . 

Our results are shown in Fig. |3.6| . The calculations have been performed using the 

700 MeV, but m 



parameters of Table 2.2 for A 



0. As discussed in section 3.2.1 



at low temperatures the model behaves again like a free pion gas (dotted line) decreasing 
with T 2 for low temperatures. Had we assumed Eq. ( 3.22 ) to be exact we should take the 
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Figure 3.7: Thermodynamic potential per volume as a function of the constituent quark 
mass in the MLA with A M = 700 MeV (left panel) atT = 163.9 MeV (dotted), T = 164.5 
MeV (dashed), and 165.3 MeV (solid) and (right panel) with Am = 10 MeV (dotted), 50 
MeV (solid), and 100 MeV (dashed) at the corresponding critical temperature. 



coefficient determining the T 2 -behavior to be equal to l/8f^°\ where f^ corresponds 
to the pion decay constant of RPA pions evaluated at the selfconsistently determined 
quark mass m! . From Table [2.2| we read off a value of f^ = 120.0 MeV for the present 
parameter set with mo = 7.9 MeV. In the chiral limit, the value of f^ is slightly smaller, 
we obtain = 118.6 MeV. The best agreement of the MLA result with the pion gas 
result can be achieved by taking f(°) = 119.0 MeV supporting in hindsight Eq. (|3.22|) . At 
first sight the approximation which has been performed to arrive at Eq. ( |3.22| ) seemed a 
little deliberate although suggested by intuition. 

Deviations from the pure pion gas behavior become visible at T « 100 MeV, which 
is quite similar to our observations in the l/iV c -expansion scheme. These deviations arise 
from quark effects which could be tolerated only close to the phase transition. In contrast 
to the l/iV c -expansion scheme an examination of the phase transition is possible in the 
MLA. The present parameter set leads to a critical temperature of T c = 164.5 MeV. Hence 
quark effects cease to be negligible at a temperature of about .6T C , where one cannot avoid 
admitting that they are completely unphysical. 

Despite this obvious shortcoming of the model let us now have a closer look at the effect 
of the mesonic degrees of freedom on the phase transition itself. The critical temperature is 
considerably lowered as compared with the corresponding value in Hartree approximation, 
T c = 266.1 MeV. This reduction is caused on the one hand by the smaller constituent 
quark mass in vacuum, w! = 447.1 MeV whereas ma — 600 MeV, and on the other 
hand by the direct influence of mesonic fluctuations on the temperature dependence of 
the quark condensate. It is worth noting that the critical temperature almost agrees with 
the value we obtain in Hartree approximation with a parameter set fitted to the RPA 
quantities. Had we succeeded in determining a parameter set in the MLA which fits the 
MLA quantities this would of course be the only sensible comparison. But, since it is to 



98 



be expected that the properties of the p-meson do not influence the quark condensate and 
therefore the critical temperature considerably, T c presumably remains in the same range 
as long as the values of (t^ip) and f n are refitted. 

In accordance with the findings in Ref. j|5) the phase transition is of first order whereas 



it is of second order in Hartree approximation. This is already indicated by the jump in 
the order parameter at T = T c . This result can be assured by considering the behavior 
of the thermodynamic potential Q, the nonzero temperature analogue to the effective 
potential (cf. Eq. (pT9|) ), 

n(m T ,T) = fi (0) (m T ,T,/i = 0) 
1 r d?p 

+ 2 T 5I J 7^^{^( 1 - 2 ^n (T (ia; ri ,p)) + 31n(l-25' s n 7r (ia; n ,p))} + const. . 

(3.25) 

The irrelevant constant is chosen in such a way that Q(0,T) = 0. Though this choice 
hinders the direct comparison of the absolute values of Q for different temperatures, it 
facilitates the identification of the order of the phase transition. That we indeed observe 
a first order phase transition can be inferred from the l.h.s. of Fig. |3.7| , where we have 
displayed Q for different temperatures as a function of the constituent quark mass 
At T = T c = 164.5 MeV one can doubtlessly identify two degenerate minima at m = 
and m/0, infallible evidence for a first order phase transition at that temperature. One 
might ask whether this finding depends on the choice of the parameters, in particular on 
the choice of Am which controls the strength of the mesonic fluctuations. One indication 
that this is not the case is the fact that the authors of Ref. E| report the same although 



they use a different parameter set. In fact, varying the meson cutoff Am we find that 
the discontinuity decreases with decreasing A M but nevertheless the order of the phase 
transition remains the same. The results for Q at the corresponding critical temperature 



are shown on the r.h.s. of Fig. |37?] for different values of the cutoff. All other parameters 
are kept constant. We conclude that a first order phase transition seems to be a property 
of the present approximation. 

This is probably, besides the unphysical quark degrees of freedom, another shortcoming 
of the approximation, because universality arguments suggest that the finite-temperature 



chiral phase transition in QCD with two massless quarks is of second order ^(J. This 
conjecture is based on the assumption that the critical behavior is completely determined 
by the four (at T c massless) bosonic degrees of freedom and that QCD therefore lies in 
the same universality class as the 0(4) model which is known to exhibit a second order 
phase transition. The same arguments can be applied to the NJL model, which has 
the same underlying symmetry as QCD with two massless flavors. However, one of the 
objections one might raise against the above hypothesis is that a theory with composite 
boson fields not necessarily belongs to the same universality class as the 0(4) model [|57|| . 
This objection is among others valid for the NJL model. 

In the same context we should mention that in the linear sigma model including bosonic 
fluctuations in one-loop calculation similarly induces a first order phase transition 



This model is usually chosen as prototype of a model being in the same universality class 
as 0(4). Applying the RPA to this model also leads to a first order phase transition [R5 . 
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Here the application of renormalization group techniques led to a second order phase 
transition Thus it seems interesting to apply these techniques also to the NJL 

model. 
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Summary 



We investigated quark and meson properties within the Nambu-Jona-Lasinio model, in- 
cluding meson-loop corrections. Here an approximation scheme beyond the usual Hartree 
+ Random Phase approximation which is consistent with prescriptions following from 
the underlying chiral symmetry of the Lagrangian is required. One possibility to con- 
struct a consistent scheme beyond Hartree + RPA is to include the ring sum in the 
functional $ within the framework of the "^-derivable method". This scheme gener- 
ates a non-local contribution to the gap equation. We decided not to perform explicit 
calculations within this scheme because of the difficulties which arise due to the non- 
trivial momentum dependence of the non-local quark self-energy term. The two following 
schemes, an expansion of the self-energies in powers of 1/N C up to next-to- leading order, 
established in Refs. pp| , [H], [12] , and a one- meson- loop approximation (MLA) to the 



effective action [3C], were applied to examine quark and meson properties. The former is 
a perturbative scheme, where one starts from the usual Hartree quark self-energy. The 
latter generates a local correction term to the self-energy. We explicitly showed that both 
schemes are in accordance with the Goldstone theorem, i.e. in the chiral limit the pions 
emerge as massless particles. In the MLA in addition the validity of the Gell-Mann- 
Oakes-Renner relation, which describes the behavior of the pion mass for small current 
quark masses, can be proved. In the i/iV c -expansion scheme the GOR relation only holds 
if one carefully expands both sides of the relation up to next-to- leading order in 1/N C . 

The non-renormalizability of the NJL model leads to additional divergencies if mesonic 
fluctuations are included. We therefore decided to introduce an additional cutoff parame- 
ter Am, which a priori is not related to the regularization of the quark-loop contributions. 
This cutoff parameter controls the strength of the mesonic fluctuations. Thus, varying 
this cutoff, we are able to study the influence of mesonic fluctuations on various quanti- 
ties. This is interesting especially in the context of the question whether chiral symmetry 
gets restored due to strong mesonic fluctuations E|. In fact, for large values of the cutoff 



Am we found instabilities in the pion propagator as well in the l/A c -expansion scheme 
as in the MLA. However, due to the perturbative treatment of the mesonic fluctuations a 
closer examination of the question whether these instabilities are related to an unstable 
ground state, as suggested in Ref . , leading to a "chiral restoration phase transition" , 
was not possible within the l/A c -expansion scheme. In the MLA an examination of the 
effective potential as a function of Am enabled us to establish a more decisive answer. We 
observed a "phase transition" only if vector interactions are included. One consequence 
of this observation is that the instabilities we find in the pion propagator are not caused 
by a "phase transition" . A plausible explanation for these instabilities could be found in 
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some peculiarities in the analytic structure of the RPA meson propagators. 

For the quark condensate in MLA we as well found a "phase transition" or a "cross- 
over" as a function of the cutoff Am, depending on the value of the current quark mass 
only if we include vector and axial intermediate states. In the other cases, i.e. in the 
1/ A^-expansion scheme and in the MLA without vector interactions, the quark condensate 
exhibits a minimum as a function of Am in the same region of cutoffs where the instabilities 
in the pion propagator emerge. Thus this observation might be related to the instabilities 
in the pion propagator. In the l/A c -expansion scheme the effect of vector interactions 
was simply that the minimum is then relocated at higher values of the cutoff. 

For the investigations described above all parameters except the mesonic cutoff Am 
were kept constant. The corresponding values have been chosen in such a way that 
the empirical values of the RPA quantities, {'tpip)^, and /^°) could be reproduced. 
However, at the end the values of the parameters, including A M , should be determined 
by fitting physical observables. In addition to the quantities in the pion sector, the 
quark condensate, the pion mass, and the pion decay constant, the properties of the p- 
meson are very well suited for that purpose since in particular pionic fluctuations are 
absolutely crucial if one wants to include the dominant p — > 7T7r-decay channel. The 
standard Hartree + RPA scheme contains only unphysical gg-decay channels. A related 
problem, which additionally constrains the possible choice of parameters is the following: 
The presence of those unphysical decay modes can by no means be avoided because they 
are a consequence of the lack of confinement in the NJL model. One way in which this 
problem can be circumvented is to increase the constituent quark mass sufficiently, such 
that this decay channel opens far above the relevant region for the p-meson, i.e. far above 
its mass. 

In the l/A c -expansion scheme a set of parameter values which fulfills all these re- 
quirements could be found. With these parameter values a reasonable fit of the data for 
the pion electromagnetic form factor, corresponding roughly to a fit of the mass and the 
width of the p- meson, together with reasonable values for {ip^),m^ and f w can be ob- 
tained. With the corresponding value of the constituent quark mass, m H = 446 MeV, the 
gg-threshold lies at ~ 892 MeV, i.e. about 120 MeV above the maximum of the p-meson 
spectral function. Moreover, with these parameter values we are far away from the region 
where we encountered instabilities in the pion propagator. In fact, from the moderate 
changes of the quantities in the pion sector - f w and (ipip) are lowered by about 20 % 
and m n is increased by about 10% - we can conclude that higher-order terms in 1/N C are 
small and that the convergence of our perturbative l/A c -expansion scheme is satisfying. 

In the MLA, on the contrary, we did not succeed in finding a suitable parameter set. 
The main reason for this failure is the fact that in this scheme the meson-loop effects 
lower the constituent quark mass as compared with the Hartree mass mn- Therefore it is 
much more difficult to prevent the p-meson from decaying considerably into qq pairs. A 
relatively large cutoff of A^ m 700 MeV was necessary in order to shift the corresponding 
threshold above the p-meson peak and simultaneously to obtain reasonable values for 
other physical observables, in particular /„.. It turned out that, applying these parameter 
values, we encounter instabilities in the p-meson propagator similar to those which emerge 
in the pion propagator in certain regions of parameter space. Conclusively it can be stated 
that we are unable to describe p-meson and pion properties reasonably at the same time 
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within the MLA. At this point one might ask whether the inclusion of further intermediate 
states, e.g. vector and axial vector mesons, could improve the situation. We can give no 
conclusive answer to that question at this stage. We can only remark that the results 
for the quark condensate and the effective potential on the one hand reveal additional 
problems, in particular a possible "phase transition", and on the other hand probably 
enlarge the parameter space where we do not encounter any instabilities. Besides, a 
different way of regularizing the occurring divergencies is worth being investigated within 
this context. 

Not only for the description of the p-meson but also for that of the a-meson it is 
absolutely necessary to include the two-pion intermediate state which can be achieved in 
the l/iV c -expansion scheme or the MLA. Our results for the a- meson, however, suffered 
from another disease: the attraction in the two-pion channel turned out to be too strong, 
leading to the absence of one physical, time-like, pole in the propagator. This problem 
is not an artifact of the NJL model but is found in the same way in similar hadronic 
calculations. The similarity to difficulties arising in hadronic calculations provided us 
with a hint, in which way this problem probably could be cured: The a- meson self-energy 
as it would come out if we applied the "<I>-derivable-method" has a promising structure 
and should therefore be examined more closely. 

As well for the description of the p- as for that of the <r-meson the two-pion intermediate 
state gives the dominant contribution which is certainly an improvement compared with 
standard calculations in the NJL model, which incorporate only quark effects. As such it is 
interesting to see how a hadronic description emerges from an underlying quark structure. 
In this context we performed a comparison of our results with that obtained in the "static 
limit". The static limit is an approximation to the full schemes where all quark effects 
are suppressed by performing an expansion in the momenta of the incoming mesons up to 
the first non-vanishing order. This approximation was motivated by the assumption that 
the relevant momenta are much larger than the constituent quark mass. This comparison 
shows that quark effects in our calculations should be regarded with care: Although we 
were able to exclude the presence of unphysical gg-decay channels in the energy region 
relevant for the description of p- or cr-meson indirect quark effects considerably influence 
also the propagator of the p-meson near its mass. 

In the last part of the present paper we examined the temperature dependence of sev- 
eral quantities. First we reviewed known results for the temperature dependence of the 
RPA quantities, in particular the masses of pion and a-meson and the spectral functions 
of p- and ax-meson. In that way one of the principal consequences of chiral symmetry 
restoration at nonzero temperature could be illustrated: Above the phase transition the 
chiral partners, i.e. ir and a or p and a±, respectively, degenerate. Then we came to ad- 
dress the temperature dependence of the quark condensate including mesonic fluctuations. 
We find that as well in the l/iV c -expansion scheme as in the MLA the low-temperature 
behavior is almost exclusively driven by thermally excited pions. Furthermore, it could 
be shown that in the chiral limit it is consistent with the lowest-order chiral perturbation 
theory result, i.e. the contributions from a free pion gas. This is certainly a consider- 
able improvement as compared with the usual Hartree approximation which completely 
fails to describe this behavior because pionic degrees of freedom, which should be the 
dominant ones in this temperature region, are missing. Within that scheme thermally 
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excited quarks are the only possible degrees of freedom. At higher temperatures, how- 
ever, we cannot avoid having effects arising from thermally excited quarks also in the 
two extended schemes. One could argue that the existence of these unphysical degrees of 
freedom is not necessarily depreciable, because they could be tolerable in the vicinity of 
the deconfmement transition, which coincides, according to lattice results, with the chiral 
phase transition. 

Of course, in the l/iV c -expansion scheme we cannot judge whether quark effects become 
visible near the chiral phase transition or much below, because the perturbative treatment 
of the mesonic fluctuations does not allow an examination of the phase transition. In 
principle this is possible within the MLA. Here the only difficulty is to suitably choose 
the parameter values since we did not succeed in determining a conclusive parameter set. 
However, the estimate we obtained by, rather deliberately, taking the parameter set for 
Am = 700 MeV seems to be robust: We find a critical temperature of 164.5 MeV, and on 
the other hand, quark effects are visible already at a temperature of ~ 100 MeV, i.e. at 
about 0.6 T c . This is much too early to be realistic, i.e. at this point the model obviously 
fails. The strong influence of vector interactions on the results for the quark condensate 
at zero temperatures suggests that this point should be reexamined including additional 
intermediate states. 



We should mention that, in agreement with Ref. ||45|| , we find a first-order phase 
transition. This point should be investigated further since it is generally believed that 
the chiral phase transition at nonzero temperature in a model with two light flavors is of 
second order. In this context the application of renormalization group techniques seems 
promising. 
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Appendix A 

Definition of elementary integrals 



It is possible to reduce the expressions for the quark loops to some elementary inte- 
grals H9l| ) see App. [B| and [D[ In this section we give the definitions of these integrals. 



1 ~ J (2tt) 4 k 2 - m 2 + ie ' { ' ' 

d 4 k 1 



I<kP) J (2tt) 4 (k 2 -m 2 + ie)((k + p) 2 - m 2 + ie) ' ^'^ 

K ^ = I (2vr) 4 (k 2 -m 2 + ie) 2 ((k + p) 2 - m 2 + ie) ' ^'^ 

r d 4 k 1 
M( Pl ,p 2 ) = J ( 2 7t) 4 (k 2 -m 2 + ie){k\ - m 2 + ie){k 2 2 - m 2 + ie) ' (A ' 4) 

L( Pl ,p 2 ,P3) = J ( 27r )4 (fc2 _ m 2 + ^2 _ m 2 + ^2 _ m 2 + ^(£.2 _ m 2 + ^) >( A ' 5 ) 

rfMxfo,!*) +rfM 1 (p a ,p 1 ) =/ ( t)4 (fc 2 _ m 2 + , e)(fc 2 _ ^ 2 + iem _ m2 + , e) >(A.6) 
with kj = k + pi. The function Mi(px,p2) can be expressed in terms of the other integrals, 

A# f \ - P 1 'P 2 ^(Pl) ~^2 / (P2) + (j^-Pl'Pi) I(P1-P2)+P2 (Pl-Pl-Pi) M (Pl,P2) 
M 1 {Pi,p 2 ) — —77 2 — 2\ ' 

(A.7) 

All integrals in Eqs. (|A~T1) to (|A~6|) , are understood to be regularized. As described in 
section pTTJ we use Pauli-Villars regularization with two regulators, i.e. we replace 



" k [ f{k;m) — > / AA.^Cjfik;^) , (A.8) 



by 



(2vr) 4 ^ ' ' J (27r) 4 ^ 

=m 2 +jA 2 q ; c = 1, a = -2, c 2 = 1 . (A.9) 

One then gets the following relatively simple analytic expressions for the integrals I\, I(p) 
and K{p): 
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h = 
I(p) 



—I 



J2 c i $ Hrf) 



16tt 2 

J 

= ^ c 3 \ x n ln ( z ji) + x j2 ln(-x j2 ) + z,i ln(-p 2 x 3l ) + x j2 \n(p 2 x j2 



I(p = 0) 



167T 2 



£ c i ln (^ 2 ) 



K(p) 



—i 



1 (i ^ 



V(^i - ^2) 



(-ln( 



(A.10) 
(A.ll) 
(A.12) 

xji) - ln(— Xji) + ln(xj 2 ) + ln(— Xj 2 )) ,(A.13) 



with 



1 1 



1 



p2 



(A.14) 



An analytic expression for the three-point function (Eq. |A.4|) can be found in Refs. p2 



and [pH|] . In certain kinematical regions the four-point function (eq. |A.5| ) is also known 
analytically ||, |3[ . 
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Appendix B 
RPA propagators 



Using the definitions given in the previous section, we can write the gap equation 
(Eq. ( P-.2|) ) in the following the form 

m = m + 2ig s AN c N f m h . (B.l) 

Similarly one can evaluate the quark-antiquark polarization diagrams (Eq. fll.ll|) ) and 
calculate the RPA meson propagators. The results for cr-meson and pion read 

° M = 1 - 2ig s 2N c N f {2h- (p 2 - 4m 2 ) I(p)) ' (R2) 

DM = l~2zg s 2N c N f (2h-pi I(p)) " (R3) 

If we evaluate these propagators with the constituent quark mass in Hartree approximation 
we can simplify the above expressions with the help of the gap equation (Eq. |B.1|) to obtain 



D (v) = — (BA) 

°W + 2igs 2NcNf {p2 _ 4m 2) /(p ) > y ■ ) 

DAP) = f + 2ig s 2N c N f p 2 I(p) ■ (R5) 



As discussed in section [2.3| , this form is also used for the internal meson propagators in 
the MLA. 

A straight-forward evaluation of the vector and axial vector polarization diagrams 
gives 

n p (p) = -i^N c Nf (-2/i + (p 2 + 2m 2 ) J(p)) , (B.6) 
n oi (p) = -4iV c iv7 {-2h + (p 2 - 4m 2 ) I{p)) . (B.7) 



we 



Because of vector current conservation Il p should vanish for p = 0. in section [TT 
discussed that in order to achieve this we have to perform a subtraction. We then obtain 
for the p- and ax-meson propagator 
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-2g v 

1 + 2ig v \N c N f (-2m 2 1(0) + (p 2 + 2m 2 ) I(p)) 

= 2^ 

1 + 2ig v ±N c N f (-2m 2 1(0) + (p 2 - Am 2 ) I(p)) 
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Appendix C 



Details concerning vector and axial 
vector intermediate states 



C.l RPA meson propagators 

One direct consequence of the presence of possible axial vector excitations is 7r-ai-mixing 
which occurs already on the level of the Hartree approximation + RPA. As mentioned 



in section |1.2| the quark-antiquark T-matrix is no longer diagonal in the pion channel. 
Therefore we have to solve a matrix equation to obtain the pion propagator in RPA. A 
generalized ansatz for the T-matrix, 

-D^M^hK^M - D^M^K^ki , (ci) 



V¥ ai ' 



leads to the following equation 



/ D ns D nsv \ f 2g 8 \ 2g s \ 11^ ^ \ / £>„, D Wsv \ 
{ A*. At. ) ^ -2g v J + \ -2g v J \ n oi7r n aiai ) ^ D nvs D„ v ) ' 

(C.2) 

For simplicity we neglected the isospin indices and in the latter expression also the mo- 
mentum arguments. We have generalized here the definition for the polarization functions, 
cf. Eq. (imp , to 

n A fiv(g) = -i J yj^ Tz[T M iS{j)+ ^)T N iS(p- -)] . (C.3) 

At this stage we should mention that our notation for the axial coupling 7 M 7sr, 
introduced in Eq. ( |1.3j ) is here somewhat misleading because this vertex is not only related 
to the coupling of an ax-meson, but also to that of several parts of the pion propagator. 
Nevertheless we will keep that notation. Note that only the longitudinal part of the axial 
polarization function contributes to the pion channel. The solution of this matrix equation 
can be written in the following way: 
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Ar„(g) 
Ar„(?) 



-2g s {l-2g v 8N c N f m 2 tI{q)) 
D(q) 

-2g s 2g v AN c N f m^iI{q) 
D(q) 

2g v (l-2g s U^(q)) 
D(q) 



(C.4) 



with the determinant 



D(q) = (l-2g s 2tN c N f (2I 1 -p 2 I(p)))(l + 2g v 8tN c N f m 2 I(q))-4g s g v m 2 q 2 (AN c N f I(q)) 2 . 

(C.5) 

With the help of the Hartree gap equation, Eq. ( |1.4j ), the integral I\, occurring several 
times in the above expressions, can be replaced analogously to the case without ir-ai- 
mixing. 

The structure of the pion propagator, Eqs. flCTT ) and ( |C.4[ ), shows that in the present 
case, including 7r-ai-mixing, it is not sufficient to define a single pion-quark coupling 
constant analogously to Eq. ( |1.19| ). One has to define a pseudoscalar and an axial pion- 
quark coupling constant by taking the residue of the propagator at the pion pole. For 
more details see e.g. Ref. [ITH . 



C.2 Consistency with chiral symmetry 

In section O] we neglected possible vector and axial vector intermediate states while 



proving the consistency of our scheme with chiral symmetry. In this section we will make 
up leeway. However, since there is conceptually no difference to the MLA we restrict the 
discussion to the l/iV c -expansion scheme. 

One has to show that, in the chiral limit, the inverse pion propagator vanishes at zero 
momentum, 

2# s fL T (0) = l for m = 0. (C.6) 

As before we use the notation II" 6 = ^II,,-. Similarly to section 1.4.3 we only need to 
show that the contributions of the correction terms add up to zero, 

<fflt fc) (0) = for m = 0. (C.7) 

k=a,b,c,d 

Let us begin by diagram 5H^\ The external pion can now couple to a 7ra, a rrp, a aa%, 
as well as to a pa\ intermediate state. Evaluating the trace in Eq. ( |1.26|) for zero external 
momentum one gets for the corresponding triangle diagrams 

1^(0, p) = ±N c N f 2m I(j>) , 

r^(0,p) = 5 ab 4iN c N f p» J(p) , 

r^(0,p) = e abc 4iN c N f p»I(p), 

^Xi c (0,p) = e abc 4N c N f 2m gT I(p) . (C.8) 
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The vertices for the coupling of the mixed or axial part of the pion propagator, cf. sec- 



tion |C.1| , can be obtained by multiplying the corresponding ai-vertex with the momentum 
of the pion, see comment on the notation below Eq. ( |U.2|) . We can see that, as the 7r, tt, p- 
and the 7r, cr, ai-vertices are proportional to the momentum of the p- and ai-meson, the 
transverse parts of the propagators do not contribute. Inserting this into Eq. ( |1.33j ) we 
find 



d A p 



{T%„„(p,p) dm r^o, - P ) dm 



+ r £^(M ^FAJp) r^ i(7 (o, - P ) dm 
+ r%z,(o,p) (-ip„)D« sv ( P ) rt%, CT (o, -p) dm 

+ K%,Mp) (iPv)D« s M 1^(0, -p) dm 
+ ^Z cd p (0, p) T pK D ai (p) r£g(0, -p) T Xu D p {p) 

+ r^(o,p) D^(p) r^(o, -p) l, u d p (o) 

+ r^S(°'f) V^-Cp) r SS(°' -p) W?„(°) 
+ r^(o lP ) {-ip»)D« av {p) rJ£J(o, -p) w?,(o) 

+ I^(0,p)(<Pi,)i? w „(p)rS^(0,-p)L^(0)} . (C.9) 
Inserting the expressions for the vertices (Eq. |1.86|) one obtains 



l6ab I 7^i( 4iVciV / J W) 2 { DM(D w M4m 2 +p 2 D nv -4mp 2 D Wsv (p)) 

+ D ai {p)D p {p)2Am 2 

- D p (0)2 (p 2 D ns (p) + Am 2 D 7T M - Amp 2 D 7Tsv (p)fp^) 



The essential step in the proof in section 1.4.3 was the validity of Eq. ( 1.89| ) which enabled 
us to write the product of the RPA pion and cx-meson propagator as a difference. Of 
course this expression cannot be directly applied to the pion propagator in the present 
case since it has been modified by vr-ax-mixing. Nevertheless we can derive a similar 
helpful expression. Besides, for the transverse part of the RPA p-meson propagator and 
the ai-meson propagator Eq. ( |1.89| ) can be adopted almost unchanged. We find 
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D a (p) D ns (p) 

DaAv) D p {p) 
to finally obtain for SU^ 



DM - DM 4m 2 



4iV c iV/ 2m 2 J(p) p 2 AN c Nf 2m? I(p) 

DM -DM 

AN c N f 2m 2 I(p) 



+ DM)> 



(C.ll) 



5n (a)a fc(()) = _^ 4jVciV/ | J^J^j^) _ DiTs{p) + 6Dp(p) _ 6/}qi(p) 



(2tt)< 

-2G„ + 2D P (0)} . 

The next two diagrams can be evaluated straightforwardly. One finds 
<5n?>» = -<U 42V e JV) / 7^{Ar(p) (/(P) + HO) - (P 2 - 4m 2 ) K(p)) 



(C.12) 



+Ajp) (3/(p) + 3/(0) 3p 2 iT(p)) 

2 

+D ff „(p) 3^4m 2 K(p) 
m 

+D nv (p) (-3/(p) - 12m 2 if(p)) 

+£> p (p) ( - 3/(p) - 6/(0) + 6(p 2 + 2m 2 ) K(p)) 

+DM ( " 3J(p)) 

+D ai (p) ( - 3/(p) - 6/(0) + 6(p 2 - 4m 2 ) K(p)) } 



<5l4 c)a6 (0) 



-5 ab 4N c Nf 



j^M{3D ai (p)-3DM-DM 

-DM - d*m + d„m} ■ 



(C.13) 



Finally we have to calculate 511^(0). According to Eq. Q1.33 ), it can be written in the 
form 

fflf(0) = -iYf^MDM A , (C.14) 

with 



d A p 

WY 4 



^2(-iD M (p))(-iT M>Mj M -P)) 

M 

4N c N f mJ-0^{ DM (2 M + 1(0) ~ (P 2 ~ 4m 2 ) K(p)) 
+D« S (p) (3/(0) - 3 P 2 K(p) ) 



p 



+D*.M -(3/(p) + 12m 2 AT(p)) 



m 



+Ar„(p) ( -6/(p) - 12m 2 tf(p) ) 

+D p (p) ( 6/(p) - 6/(0) + 6(p 2 + 2m 2 ) K(p) ) 

+D ai (p) ( -12/(p) - 6/(0) + 6(p 2 - 4m 2 ) K(p) )\Clh) 
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Evaluating D a (0) in the chiral limit at the Hartree mass mn and comparing the result 
with Eq. ( |C.8f) one finds that the product of the first two factors in Eq. (|C.14|) is simply 
Sab/n^H, i-e. one gets 

611® ab (0) = 6 ab — . (C.16) 

With these results one can easily check that Eq. (|C.7|) indeed holds in our scheme. 

It now remains to check whether the replacement for the RPA meson propagators in the 
MLA, see section |2.3| , does preserve the symmetry properties of the schemes also if vector 
and axial vector intermediate states are included. The question is whether Eq. ( |C.11|) 
still holds for the RPA pion and cr-meson propagators. The analogue of the replacement 
in Eq. ( [2.16| ) would be to replace the determinant in the denominator of the RPA pion 
propagator, Eq. (|C.5|) , by 



D(q) = + 2g s 2iN c N f p 2 I(p)))(l + 2g v 8tN c N f m 2 I(q)) - Ag s g v m 2 q 2 (AN c N f I(q)) 2 . 

(C.17) 

Inserting this into Eq. ( |C.4j ) for the pion propagator one can easily verify that Eq. ( |C . 1 1| ) 



still holds for pion and sigma. Because of gauge invariance we had to perform a subtraction 
for the vector and axial vector polarization function in RPA anyway, such that no more 
replacement is needed here. This subtraction is consistent with Eq. (|C.11|) . 



C.3 Effective potential including Rho- and ai-mesons 

The explicit formula for the effective potential in vacuum including vector interactions 
reads 

d A p m 2 — p 2 (m - Wq) 2 



V(m) = -UN c N f [ T^hi^V^) + 
J (27T) 4 m„- v 2 

i f d A p 



(2vr) 4 V mg-P 2 *9. 
2 I J2W { H-D- l ip)/2g s ) + StvH-D- 1 ^) 

+ 3\n(-D; 1 (p)/2g v ) + 3 \n(-D~ i 1 (p)/2g v )} + const. ,(C18) 



where the pion propagator has to be taken from Eq. (|C2j ). The trace is here to be 
understood to sum only the components of the 2 x 2matrix in Eq. ( |C2[ ). 
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Appendix D 



Explicit expressions for the 
meson-meson vertices 



In this section we list the explicit formulae for the meson-meson vertices. We restrict 
ourselves to those combinations which are needed for the calculations presented in this 
paper. We introduce the convention that a vertex labelled with 'V is to be understood to 
contain a pseudoscalar coupling, whereas a vertex labelled with "ai" is to be understood 
to be calculated with an axial coupling although the latter, due to 7r-a!-mixing, couples 
also pionic states. 

We begin with the three-meson vertices T Ml ^ M2 ^ M3 (q,p) (see Fig. L4|) : 



h(q,p 

h(q,p 
f2(q,P 

-iK:t P (Q,P 

f(q,p 



i2mN[I{p') + I(q) + I(p) + (4m 2 - ~(p' 2 + p 2 + q 2 ))M(p, -q) 



i2mN6 ab (l(p')+p-qM(p, -q\ 
Oabh{q,p)[g* 



p 2 q 2 



(p-qf 



imN(l(q) + I(p) - 2I(p') + (4m 2 -2p-q- p' 2 )M(p, -<?)) 
p2 _|_ q2pfj.p\ _ (p' M q x + q^p x ) 

S * h &*W p^ - (p. g )2 

-2imN5 ab g^(j(p) + I(q) - (p' 2 - 4m 2 )M(p, -qf 

$ab(q> M fi{q,p)+P fl f2{q,p)) , 
-N(l(p') -p 2 M(p, -q) + 2p-p'M 1 (q, -p) 
-N(l(p') + J(p) - (q 2 - Am 2 )M(p, -q) - 2p-p'M 1 (-p, qfj , 
£abc(q'*f(q,p)-p tl f(p,qj) , 

N(- I{q) +p 2 M{p,-q) +2p-qM l (q,-pj) , (D.i; 



with p' = —p — q and N = 4N c Nf. To calculate the constant A we need the vertices 
coupling to a for the special case where the momentum of the a is vanishing. Since 
this limit is in several cases difficult to build, we will also list here the corresponding 
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expressions, 





-P,p) 


-jV ab ( 


~P,p) 




~P,p) 


ai,oi,tr\ 


-P,p) 



,2 2\ 



p) = 5 ab ^fT^(l(p)-I(0) + (p 2 + 2m 2 )K(p) 



+ pV(7(p) - 7(0) + (p 2 + 2m 2 )K{p))) , 
-zr^ )CT (-p,p) = -iV^(/(p)+4m 2 K(p)) . 

(D.2) 

For the four-meson vertices we only need to consider the special cases needed for the 
diagrams (b) and (c) in Fig. |TT3| : 

- iT am (q,p, -q) = ~N{l(p - q) + I(p + q) + 8m 2 (M(p, q) + M(p, -q)) 



+4(m 2 (4m 2 — p 2 — q 2 ) — — j-) £(p, — g, P — g) } 



-iT^^iq, p, -p) = -iV{/(p + g) + 7(0) + 4m 2 (7i» + K(q) + 2M(p, -g)) 

+2p-qM(p, -q) - q 2 K(q) - p 2 K(p) 

+m 2 (16m 2 - Ap 2 - Aq 2 + ?^L)L(p, -q, 0)} 
-iTl b m (q,p,-q) = 5 ab N{l(p + q) + I(p-q)+p 2 (Am 2 -q 2 )L(p,-q,p-q)} 
-iTl b m (q,p, -p) = S ab N{ - J(p + g) - J(0) - (4m 2 - q 2 )(K(q) - p 2 L(p, -q, 0)) 

+p 2 K(p) - 2p-g M(p, -g)| 
-ir^ !7r (g,p,-g) = -7V« obcd {/(p + g)+J(p-g)-^VL(p,-g,p-g)} 
-iT^(q,p,-p) = -NK abcd {l(p + q) + I(0)-p 2 K(p) 

-q 2 K{q) + 2p-qM(p, -q)+p 2 q 2 L(p, -g, 0)} 
-ir^ )CT (g,p, -q) = -25 ab N{l(p + q) + J(p - g) + 27(g) - p-g(M(p, -g) - M(p, g)) 

+ (4m 2 - 2p 2 )(M(p, q) + M(p, -g)) 

+m 2 (8m 2 - 6p 2 + 4g 2 + P ~ ^'^ )L(p, -g,p - g)} 

-ir^ )CT)P (g,p,-g) = -25 afc iV{-7(p + g)-7(0) + (p 2 -4m 2 )ir(p) 

+ (g 2 + 2m 2 )K(q) + (4m 2 - 2p-g)M(p, -g) 

+m 2 (8m 2 - 2p 2 + 4g 2 - ^-)L(p, -g, 0)} 



m z 
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-iTf^{q,p,-q) = 2NK abcd {-I(p + q)-I(p-q)-2I(q) 

+2p 2 (M(p, q) + M(p, -q)) + p-q(M(p, -q) - M(p, q)) 
+ {2m 2 p 2 - p 4 + (p-q) 2 )L(p, -q,p - g)| 

-ir^% p (q,p,-q) = 2NK abcd {l(p + q) + I(0)-p 2 K(p)-(q 2 + 2m 2 )K(q) 

+2p-qM(p, -q)+p 2 (2m 2 + q 2 )L(p, -q,0)} , (D.3) 

with T PyMMyP (q,p,-q) = g flu T^ MM ^q,p,-q),T pM ^ M (q,p,-p) = ^r^ )P)Af (g,p, ~p) 
and K abcd = 5 ab 5 cd + S ad 5 bc — S ac S bd . 
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Appendix E 



Expressions at nonzero temperature 



E.l Feyman rules at nonzero temperature and chem- 
ical potential 

• Quark propagator: 

S{tu! n ,p) = — — — , (E.l) 



with cj„ = (2n+ \)ixT fermionic Matsubara frequency, \i the chemical potential and 
E = y/fP + m 2 . 

• Vertices: 

-2g M T M ®r M (E.2) 



X 



A factor (-1) and a trace for closed fermion loops 

The four-dimensional momentum space integration is converted into a three- 
dimensional one and a sum over the corresponding Matsubara frequencies, i. e. 

^ rW^. (E.3) 



(2tt) 4 (2k) 

E.2 Rho- and ai-meson propagators in the medium 

In this section we will explain the structure of the p- and ai-meson propagators in medium. 
Denoting the four-vector describing the motion of the heat bath with -u M , we have in 
contrast to the vacuum here two independent four-vectors, and the momentum g M of 
the meson, from which we can construct the tensor structure of the polarization functions. 
Generally we can write 

n"" = <r /i + ^J-/2 + «"«Vs + q q h , (E.4) 

q z q ■ u 
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where we have already used the symmetry of the polarization function with respect to an 
exchange of ji and v. The functions fi to / 4 are, a priori independent, scalar functions, 
depending on the Lorentz scalars q 2 , and q ■ u. The most simple choice for the four- 
vector is = (1,0,0,0), i.e. we assume to be in the rest frame of the heat bath. 
With that choice the general polarization function can be splited into a four-dimensional 
transverse and three-dimensional transverse, a four-dimensional transverse and three- 
dimensional longitudinal as well as two four-dimensional longitudinal components. Here 
three-dimensional transverse or longitudinal, respectively, means that the polarization 
function is transverse (longitudinal) with respect to q. Since the vector polarization 
function is in our case four-dimensional transverse and the a± is constructed exclusively 
from the four-dimensional transverse part of the axial polarization function we can restrict 
the discussion to that part. The two projectors on the three-dimensional transverse, P^ T , 
and longitudinal, P^li P ar t are > 



if i,j = 1,2,3 
otherwise 

pp L = pp-Pg. (E.5) 




In the rest frame of the meson, i.e. for q = 0, the two parts of the polarization function 
are equal and can be extracted in the same way as in vacuum with the help of the P^ ■ 
For instance, for the propagators of p- and ai-meson in RPA we then arrive at 

Dp{qo > 0) = 1 + 2ig v lN c N f (-2m2 1(0,0) + (g 2 + 2m*) I(q , 0)) ' (K6) 
jDai(g °'° ) = 1 + 2ig v lN c N f (-2m2 1(0, 0) + (q 2 - 4m 2 ) I(q , 0)) ' (K7) 
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